Engineering Vibration 
Analysis with Application to 
Control Systems 



C. F. Beards BSc, PhD, C Eng, MRAeS, MIOA 

Consultant in Dynamics, Noise and Vibration 
Formerly of the Department of Mechanical Engineering 
Imperial College of Science, Technology and Medicine 
University of London 



Edward Arnold 

A member of the Hodder Headline Group 
LONDON SYDNEY AUCKLAND 





First published in Great Britain 1995 by 

Edward Arnold, a division of Hodder Headline PLC, 

338 Euston Road, London NW1 3BH 

© 1995 C. F. Beards 

All rights reserved. No part of this publication may be reproduced or 
transmitted in any form or by any means, electronically or mechanically, 
including photocopying, recording or any information storage or retrieval 
system, without either prior permission in writing from the publisher or a 
licence permitting restricted copying. In the United Kingdom such licences 
are issued by the Copyright Licensing Agency: 90 Tottenham Court Road, 
London W1P 9HE. 

Whilst the advice and information in this book is believed to be true and 
accurate at the date of going to press, neither the author nor the publisher 
can accept any legal responsibility or liability for any errors or omissions 
that may be made. 

British Library Cataloguing in Publication Data 
A catalogue record for this book is available from the British Library 

ISBN 0 340 63183 X 

12345 95 96 97 98 99 

Typeset in 10 on 12pt Times by 

PPS Limited, London Road, Amesbury, Wilts. 

Printed and bound in Great Britain by 
J W Arrowsmith Ltd, Bristol 




Acknowledgements 



Some of the problems first appeared in University of London B.Sc. (Eng) Degree 
Examinations, set for students of Imperial College, London. The section on random 
vibration has been reproduced with permission from the Mechanical Engineers 
Reference Book, 12th edn, Butterworth — Heinemann, 1993. 




Preface 



The high cost and questionable supply of many materials, land and other resources, 
together with the sophisticated analysis and manufacturing methods now available, 
have resulted in the construction of many highly stressed and lightweight machines 
and structures, frequently with high energy sources, which have severe vibration 
problems. Often, these dynamic systems also operate under hostile environmental 
conditions and with minimum maintenance. It is to be expected that even higher 
performance levels will be demanded of all dynamic systems in the future, together 
with increasingly stringent performance requirement parameters such as low noise 
and vibration levels, ideal control system responses and low costs. In addition it is 
widely accepted that low vibration levels are necessary for the smooth and quiet 
running of machines, structures and all dynamic systems. This is a highly desirable 
and sought after feature which enhances any system and increases its perceived quality 
and value, so it is essential that the causes, effects and control of the vibration of 
engineering systems are clearly understood in order that effective analysis, design and 
modification may be carried out. That is, the demands made on many present day 
systems are so severe, that the analysis and assessment of the dynamic performance 
is now an essential and very important part of the design. Dynamic analysis is 
performed so that the system response to the expected excitation can be predicted 
and modifications made as required. This is necessary to control the dynamic response 
parameters such as vibration levels, stresses, fatigue, noise and resonance. It is also 
necessary to be able to analyse existing systems when considering the effects of 
modifications and searching for performance improvement. 

There is therefore a great need for all practising designers, engineers and scientists, 
as well as students, to have a good understanding of the analysis methods used for 
predicting the vibration response of a system, and methods for determining control 
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system performance. It is also essential to be able to understand, and contribute to, 
published and quoted data in this field including the use of, and understanding of, 
computer programs. 

There is great benefit to be gained by studying the analysis of vibrating systems 
and control system dynamics together, and in having this information in a single 
text, since the analyses of the vibration of elastic systems and the dynamics of control 
systems are closely linked. This is because in many cases the same equations of motion 
occur in the analysis of vibrating systems as in control systems, and thus the techniques 
and results developed in the analysis of one system may be applied to the other. It 
is therefore a very efficient way of studying vibration and control. 

This has been successfully demonstrated in my previous books Vibration Analysis 
and Control System Dynamics (1981) and Vibrations and Control Systems (1988). 
Favourable reaction to these books and friendly encouragement from fellow 
academics, co-workers, students and my publisher has led me to write Engineering 
Vibration Analysis with Application to Control Systems. 

Whilst I have adopted a similar approach in this book to that which I used 
previously, I have taken the opportunity to revise, modify, update and expand the 
material and the title reflects this. This new book discusses very comprehensively the 
analysis of the vibration of dynamic systems and then shows how the techniques and 
results obtained in vibration analysis may be applied to the study of control system 
dynamics. There are now 75 worked examples included, which amplify and demon- 
strate the analytical principles and techniques so that the text is at the same time 
more comprehensive and even easier to follow and understand than the earlier books. 
Furthermore, worked solutions and answers to most of the 130 or so problems set 
are included. (I trust that readers will try the problems before looking up the worked 
solutions in order to gain the greatest benefit from this.) 

Excellent advanced specialised texts on engineering vibration analysis and control 
systems are available, and some are referred to in the text and in the bibliography, 
but they require advanced mathematical knowledge and understanding of dynamics, 
and often refer to idealised systems rather than to mathematical models of real systems. 
This book links basic dynamic analysis with these advanced texts, paying particular 
attention to the mathematical modelling and analysis of real systems and the 
interpretation of the results. It therefore gives an introduction to advanced and 
specialised analysis methods, and also describes how system parameters can be 
changed to achieve a desired dynamic performance. 

The book is intended to give practising engineers, and scientists as well as students 
of engineering and science to first degree level, a thorough understanding of the 
principles and techniques involved in the analysis of vibrations and how they can 
also be applied to the analysis of control system dynamics. In addition it provides a 
sound theoretical basis for further study. 



Chris Beards 
January 1995 
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General notation 



a damping factor, 

dimension, 
displacement. 

b circular frequency (rad/s), 

dimension, 
port coefficient. 

c coefficient of viscous damping, 

velocity of propagation of stress wave. 
c c coefficient of critical viscous damping = 2^J(mk). 

c d equivalent viscous damping coefficient for dry friction 

damping = AFJnaiX. 

c H equivalent viscous damping coefficient for hysteretic damping = rjk/a>. 

d diameter. 

/ frequency (Hz), 

exciting force. 

/ s Strouhal frequency (Hz). 

g acceleration constant. 

h height, 

thickness. 

j V- 1- 

k linear spring stiffness, 

beam shear constant, 
gain factor. 

torsional spring stiffness, 
complex stiffness = k( 1 + jr/). 
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xvi General notation 



/ length. 

m mass. 

q generalized coordinate, 

r radius. 

s Laplace operator = a + }b. 

t time. 

u displacement. 

v velocity, 

deflection. 

x displacement. 

y displacement, 

z displacement. 

A amplitude, 

constant, 

cross-sectional area, 
constant. 

.2.3.4 constants, 

flexural rigidity = £/i 3 /12(l — v 2 ), 
hydraulic mean diameter, 
derivative w.r.t. time, 
modulus of elasticity, 
in-phase, or storage modulus, 
quadrature, or loss modulus, 
complex modulus = £' 4- j£". 
exciting force amplitude. 

Coulomb (dry) friction force (fi N). 
transmitted force, 
centre of mass, 
modulus of rigidity, 
gain factor. 

mass moment of inertia, 
second moment of area, 
moment of inertia, 
stiffness, 
gain factor, 
length. 

Laplace transform, 
mass, 
moment, 
mobility. 

applied normal force, 
gear ratio, 
force. 

factor of damping, 
flow rate. 
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Q { generalized external force. 

R radius of curvature. 

[S’] system matrix. 

T kinetic energy, 

tension, 
time constant. 

7^ transmissibility = F T /F. 

V potential energy, 

speed. 

X amplitude of motion. 

{A'} column matrix. 

X s static deflection = F/k. 

X/X s dynamic magnification factor. 

Z impedance, 

a coefficient, 

influence coefficient, 
phase angle, 
receptance. 

P coefficient, 

receptance. 
y coefficient, 

receptance. 

S deflection. 

e short time, 

strain. 

£ 0 strain amplitude. 

C damping ratio = c/c e 

>1 loss factor = £"/£'. 

9 angular displacement, 

slope. 

/ matrix eigenvalue, 

lpA(o 2 /Eiy iA 

p coefficient of friction, 

mass ratio - m/M. 
v Poisson’s ratio, 

circular exciting frequency (rad/s). 
p material density, 

a stress. 

a 0 stress amplitude. 

t period of vibration = \/f. 

t d period of dry friction damped vibration. 

t v period of viscous damped vibration. 

4> phase angle, 

function of time, 
angular displacement. 
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\j/ phase angle. 

a) undamped circular frequency (rad/s). 

cd a dry friction damped circular frequency. 

co v viscous damped circular frequency = co,/(l — ( 2 ). 

A logarithmic decrement = In X l /X 11 . 

$ transfer function. 

Q natural circular frequency (rad/s). 
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Introduction 



The vibration which occurs in most machines, vehicles, structures, buildings and 
dynamic systems is undesirable, not only because of the resulting unpleasant motions 
and the dynamic stresses which may lead to fatigue and failure of the structure or 
machine, and the energy losses and reduction in performance which accompany 
vibrations, but also because of the noise produced. Noise is generally considered to 
be unwanted sound, and since sound is produced by some source of motion or 
vibration causing pressure changes which propagate through the air or other 
transmitting medium, vibration control is of fundamental importance to sound 
attenuation. Vibration analysis of machines and structures is therefore often a 
necessary prerequisite for controlling not only vibration but also noise. 

Until early this century, machines and structures usually had very high mass and 
damping, because heavy beams, timbers, castings and stonework were used in their 
construction. Since the vibration excitation sources were often small in magnitude, 
the dynamic response of these highly damped machines was low. However, with the 
development of strong lightweight materials, increased knowledge of material 
properties and structural loading, and improved analysis and design techniques, the 
mass of machines and structures built to fulfil a particular function has decreased. 
Furthermore, the efficiency and speed of machinery have increased so that the 
vibration exciting forces are higher, and dynamic systems often contain high energy 
sources which can create intense noise and vibration problems. This process of 
increasing excitation with reducing machine mass and damping has continued at an 
increasing rate to the present day when few, if any, machines can be designed without 
carrying out the necessary vibration analysis, if their dynamic performance is to be 
acceptable. The demands made on machinery, structures, and dynamic systems are 
also increasing, so that the dynamic performance requirements are always rising. 
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There have been very many cases of systems failing or not meeting performance 
targets because of resonance, fatigue, excessive vibration of one component or another, 
or high noise levels. Because of the very serious effects which unwanted vibrations 
can have on dynamic systems, it is essential that vibration analysis be carried out as 
an inherent part of their design, when necessary modifications can most easily be 
made to eliminate vibration, or at least to reduce it as much as possible. However, 
it must also be recognized that it may sometimes be necessary to reduce the vibration 
of an existing machine, either because of inadequate initial design, or by a change in 
function of the machine, or by a change in environmental conditions or performance 
requirements, or by a revision of acceptable noise levels. Therefore techniques for the 
analysis of vibration in dynamic systems should be applicable to existing systems as 
well as those in the design stage; it is the solution to the vibration or noise problem 
which may be different, depending on whether or not the system already exists. 

There are two factors which control the amplitude and frequency of vibration of 
a dynamic system: these are the excitation applied and the dynamic characteristics 
of the system. Changing either the excitation or the dynamic characteristics will 
change the vibration response stimulated in the system. The excitation arises from 
external sources, and these forces or motions may be periodic, harmonic or random 
in nature, or arise from shock or impulsive loadings. 

To summarize, present-day machines and structures often contain high-energy 
sources which create intense vibration excitation problems, and modern construction 
methods result in systems with low mass and low inherent damping. Therefore careful 
design and analysis is necessary to avoid resonance or an undesirable dynamic 
performance. 

The demands made on automatic control systems are also increasing. Systems are 
becoming larger and more complex, whilst improved performance criteria, such as 
reduced response time and error, are demanded. Whatever the duty of the system, from 
the control of factory heating levels to satellite tracking, or from engine fuel control 
to controlling sheet thickness in a steel rolling mill, there is continual effort to improve 
performance whilst making the system cheaper, more efficient, and more compact. 
These developments have been greatly aided in recent years by the wide availability 
of microprocessors. Accurate and relevant analysis of control system dynamics is 
necessary in order to determine the response of new system designs, as well as to predict 
the effects of proposed modifications on the response of an existing system, or to 
determine the modifications necessary to enable a system to give the required response. 

There are two reasons why it is desirable to study vibration analysis and the 
dynamics of control systems together as dynamic analysis. Firstly, because control 
systems can then be considered in relation to mechanical engineering using mechanical 
analogies, rather than as a specialized and isolated aspect of electrical engineering, 
and secondly, because the basic equations governing the behaviour of vibration and 
control systems are the same: different emphasis is placed on the different forms of 
the solution available, but they are all dynamic systems. Each analysis system benefits 
from the techniques developed in the other. 

Dynamic analysis can be carried out most conveniently by adopting the following 
three-stage approach: 
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Stage I. Devise a mathematical or physical model of the system to be analysed. 
Stage II. From the model, write the equations of motion. 

Stage III. Evaluate the system response to relevant specific excitation. 

These stages will now be discussed in greater detail. 



Stage I. The mathematical model 

Although it may be possible to analyse the complete dynamic system being considered, 
this often leads to a very complicated analysis, and the production of much unwanted 
information. A simplified mathematical model of the system is therefore usually sought 
which will, when analysed, produce the desired information as economically as possible 
and with acceptable accuracy. The derivation of a simple mathematical model to 
represent the dynamics of a real system is not easy, if the model is to give useful and 
realistic information. 

However, to model any real system a number of simplifying assumptions can often 
be made. For example, a distributed mass may be considered as a lumped mass, or 
the effect of damping in the system may be ignored particularly if only resonance 




Fig. 1.1. Rover 800 front suspension. (By courtesy of Rover Group.) 
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frequencies are needed or the dynamic response required at frequencies well away 
from a resonance, or a non-linear spring may be considered linear over a limited 
range of extension, or certain elements and forces may be ignored completely if their 
effect is likely to be small. Furthermore, the directions of motion of the mass elements 
are usually restrained to those of immediate interest to the analyst. 

Thus the model is usually a compromise between a simple representation which 
is easy to analyse but may not be very accurate, and a complicated but more realistic 
model which is difficult to analyse but gives more useful results. Consider for example, 
the analysis of the vibration of the front wheel of a motor car. Fig. 1.1 shows a typical 
suspension system. As the car travels over a rough road surface, the wheel moves up 
and down, following the contours of the road. This movement is transmitted to the 
upper and lower arms, which pivot about their inner mountings, causing the coil 




Fig. 1.2(a). Simplest model - motion in a 
vertical direction only can be analysed. 



Fig. 1.2(b). Motion in a vertical direction 
only can be analysed. 




Body mass 

Damper 
Spring 



r Unsprung mass 
(front wheel) 
Tyre damping 
Tyre stiffness 



Fig. 1.2(c). Motion in a vertical direction, roll, and pitch can be analysed. 
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spring to compress and extend. The action of the spring isolates the body from the 
movement of the wheel, with the shock absorber or damper absorbing vibration and 
sudden shocks. The tie rod controls longitudinal movement of the suspension unit. 

Fig. 1.2(a) is a very simple model of this same system, which considers translational 
motion in a vertical direction only: this model is not going to give much useful 
information, although it is easy to analyse. The more complicated model shown in 
Fig. 1.2(b) is capable of pro ducing some meaningful results at the cost of increased 
labour in the analysis, but the analysis is still confined to motion in a vertical direction 
only. A more refined model, shown in Fig. 1.2(c), shows the whole car considered, 
translational and rotational motion of the car body being allowed. 

If the modelling of the car body by a rigid mass is too crude to be acceptable, a 
finite element analysis may prove useful. This technique would allow the body to be 
represented by a number of mass elements. 

The vibration of a machine tool such as a lathe can be analysed by modelling the 
machine structure by the two degree of freedom system shown in Fig. 1.3. In the 
simplest analysis the bed can be considered to be a rigid body with mass and inertia, 
and the headstock and tailstock are each modelled by lumped masses. The bed is 
supported by springs at each end as shown. Such a model would be useful for 
determining the lowest or fundamental natural frequency of vibration. A refinement 
to this model, which may be essential in some designs of machine where the bed 
cannot be considered rigid, is to consider the bed to be a flexible beam with lumped 
masses attached as before. 




Machine 

tool 




Two degree of 
freedom model 



Fig. 1.3. Machine tool vibration analysis model. 





6 Introduction 



[Ch. 1 




Radio telescope in 
vertical position 



Five degree of freedom 
torsional vibration model 



Fig. 1.4. Radio telescope vibration analysis model. 



To analyse the torsional vibration of a radio telescope when in the vertical position 
a five degree of freedom model, as shown in Fig. 1.4, can be used. The mass and 
inertia of the various components may usually be estimated fairly accurately, but the 
calculation of the stiffness parameters at the design stage may be difficult; fortunately 
the natural frequencies are proportional to the square root of the stiffness. If the 
structure, or a similar one, is already built, the stiffness parameters can be measured. 
A further simplification of the model would be to put the turret inertia equal to zero, 
so that a three degree of freedom model is obtained. Such a model would be easy to 
analyse and would predict the lowest natural frequency of torsional vibration with 
fair accuracy, providing the correct inertia and stiffness parameters were used. It 
could not be used for predicting any other modes of vibration because of the coarseness 
of the model. However, in many structures only the lowest natural frequency is 
required, since if the structure can survive the amplitudes and stresses at this frequency 
it will be able to survive other natural frequencies too. 

None of these models include the effect of damping in the structure. Damping in 
most structures is very low so that the difference between the undamped and the 
damped natural frequencies is negligible. It is usually only necessary to include the 
effects of damping in the model if the response to a specific excitation is sought, 
particularly at frequencies in the region of a resonance. 

A block diagram model is usually used in the analysis of control systems. For 
example, a system used for controlling the rotation and position of a turntable about 







Feedback loop 



Fig. 1.6. Turntable position control system: block diagram model. 

It can be seen that the feedback loop enables the input and output positions to 
be compared, and the error signal, if any, is used to activate the motor and hence 
rotate the turntable until the error signal is zero; that is, the actual position and the 
desired position are the same. 

The model parameters 

Because of the approximate nature of most models, whereby small effects are neglected 
and the environment is made independent of the system motions, it is usually 
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reasonable to assume constant parameters and linear relationships. This means that 
the coefficients in the equations of motion are constant and the equations themselves 
are linear: these are real aids to simplifying the analysis. Distributed masses can often 
be replaced by lumped mass elements to give ordinary rather that partial differential 
equations of motion. Usually the numerical value of the parameters can, substantially, 
be obtained directly from the system being analysed. However, model system 
parameters are sometimes difficult to assess, and then an intuitive estimate is required, 
engineering judgement being of the essence. 

It is not easy to create a relevant mathematical model of the system to be analysed, 
but such a model does have to be produced before Stage II of the analysis can be 
started. Most of the material in subsequent chapters is presented to make the reader 
competent to carry out the analyses described in Stages II and III. A full understanding 
of these methods will be found to be of great help in formulating the mathematical 
model referred to above in Stage I. 

Stage II. The equations of motion 

Several methods are available for obtaining the equations of motion from the 
mathematical model, the choice of method often depending on the particular model 
and personal preference. For example, analysis of the free-body diagrams drawn for 
each body of the model usually produces the equations of motion quickly: but it can 
be advantageous in some cases to use an energy method such as the Lagrange equation. 

From the equations of motion the characteristic or frequency equation is obtained, 
yielding data on the natural frequencies, modes of vibration, general response, and 
stability. 

Stage III. Response to specific excitation 

Although Stage II of the analysis gives much useful information on natural frequencies, 
response, and stability, it does not give the actual system response to specific 
excitations. It is necessary to know the actual response in order to determine such 
quantities as dynamic stress, noise, output position, or steady-state error for a range 
of system inputs, either force or motion, including harmonic, step and ramp. This is 
achieved by solving the equations of motion with the excitation function present. 





Stage I. 


Model 


Remember: 


Stage II. 


Equations 




Stage III. 


Excitation 



A few examples have been given above to show how real systems can be modelled, 
and the principles of their analysis. To be competent to analyse system models it is 
first necessary to study the analysis of damped and undamped, free and forced 
vibration of single degree of freedom systems such as those discussed in Chapter 2. 
This not only allows the analysis of a wide range of problems to be carried out, but 
it is also essential background to the analysis of systems with more than one degree 
of freedom, which is considered in Chapter 3. Systems with distributed mass, such 
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as beams, are analysed in Chapter 4. Some aspects of automatic control system 
analysis which require special consideration, particularly their stability and system 
frequency response, are discussed in Chapter 5. Each of these chapters includes worked 
examples to aid understanding of the theory and techniques described, whilst Chapter 
6 contains a number of problems for the reader to try. Chapter 7 contains answers 
and worked solutions to most of the problems in Chapter 6. A comprehensive 
bibliography and an index are included. 




2 



The vibrations of systems having one degree of 
freedom 



All real systems consist of an infinite number of elastically connected mass elements 
and therefore have an infinite number of degrees of freedom; and hence an infinite 
number of coordinates are needed to describe their motion. This leads to elaborate 
equations of motion and lengthy analyses. However, the motion of a system is often 
such that only a few coordinates are necessary to describe its motion. This is because 
the displacements of the other coordinates are restrained or not excited, so that they 
are so small that they can be neglected. Now, the analysis of a system with a few 
degrees of freedom is generally easier to carry out than the analysis of a system with 
many degrees of freedom, and therefore only a simple mathematical model of a system 
is desirable from an analysis viewpoint. Although the amount of information that a 
simple model can yield is limited, if it is sufficient then the simple model is adequate 
for the analysis. Often a compromise has to be reached, between a comprehensive 
and elaborate multi-degree of freedom model of a system, which is difficult and costly 
to analyse but yields much detailed and accurate information, and a simple few 
degrees of freedom model that is easy and cheap to analyse but yields less information. 
However, adequate information about the vibration of a system can often be gained 
by analysing a simple model, at least in the first instance. 

The vibration of some dynamic systems can be analysed by considering them as 
a one degree or single degree of freedom system; that is a system where only one 
coordinate is necessary to describe the motion. Other motions may occur, but they 
are assumed to be negligible compared to the coordinate considered. 

A system with one degree of freedom is the simplest case to analyse because only 
one coordinate is necessary to completely describe the motion of the system. Some 
real systems can be modelled in this way, either because the excitation of the system 
is such that the vibration can be described by one coordinate although the system 
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could vibrate in other directions if so excited, or the system really is simple, as for 
example a clock pendulum. It should also be noted that a one degree of freedom 
model of a complicated system can often be constructed where the analysis of a 
particular mode of vibration is to be carried out. To be able to analyse one degree 
of freedom systems is therefore an essential ability in vibration analysis. Furthermore, 
many of the techniques developed in single degree of freedom analysis are applicable 
to more complicated systems. 

2.1. FREE UNDAMPED VIBRATION 
2.1.1 Translation vibration 

In the system shown in Fig. 2.1 a body of mass m is free to move along a fixed 
horizontal surface. A spring of constant stiffness k which is fixed at one end is attached 
at the other end to the body. Displacing the body to the right (say) from the equilibrium 
position causes a spring force to the left (a restoring force). Upon release this force 
gives the body an acceleration to the left. When the body reaches its equilibrium 
position the spring force is zero, but the body has a velocity which carries it further 
to the left although it is retarded by the spring force which now acts to the right. 
When the body is arrested by the spring the spring force is to the right so that the 
body moves to the right, past its equilibrium position, and hence reaches its initial 
displaced position. In practice this position will not quite be reached because damping 
in the system will have dissipated some of the vibrational energy. However, if the 
damping is small its effect can be neglected. 




Fig. 2.1. Single degree of freedom model - translation vibration. 

If the body is displaced a distance x 0 to the right and released, the free-body 
diagrams (FBD's) for a general displacement x are as shown in Figs. 2.2(a) and (b). 




Fig. 2.2. (a) Applied force; (b) effective force. 
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The effective force is always in the direction of positive x. If the body is being 
retarded x will be calculated to be negative. The mass of the body is assumed constant: 
this is usually so, but not always, as for example in the case of a rocket burning fuel. 
The spring stiffness k is assumed constant: this is usually so within limits; see section 
2.1.3. It is assumed that the mass of the spring is negligible compared to the mass of 
the body; cases where this is not so are considered in section 4.3.1. 

From the free-body diagrams the equation of motion for the system is 

m'x = — kx or x + ( k/m)x = 0. (2.1) 

This will be recognized as the equation for simple harmonic motion. The solution is 

x = A cos cot + B sin cot, (2.2) 

where A and B are constants which can be found by considering the initial conditions, 
and co is the circular frequency of the motion. Substituting (2.2) into (2.1) we get 

— co 2 (A cos cot + B sin cot) + (/c/m) (A cos cot + B sin cot) = 0. 

Since ( A cos cot + B sin cot)^0 (otherwise no motion), 

co = yj(k/m) rad/s, 

and 

x = A cos yj{klm)t + B sin J(k/m)t. 

Now 

x = x 0 at t = 0, 

thus 

x 0 = A cos 0 + B sin 0, and therefore x 0 = A, 

and 



x = 0 at t = 0, 

thus 

0 = — Aj(k/m) sin 0 + Bj(k/m) cos 0, and therefore B — 0; 

that is, 



x = x 0 cos yj(k/m)t. 



(2.3) 



The system parameters control co and the type of motion but not the amplitude 
x 0 , which is found from the initial conditions. The mass of the body is important, its 
weight is not, so that for a given system, co is independent of the local gravitational field. 
The frequency of vibration,/, is given by 



/ = 



co 

2 / 



or 




(2.4) 
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The motion is as shown in Fig. 2.3. 




The period of the oscillation, t, is the time taken for one complete cycle so that 



T 



1 

7 



= 2n^J(mjk) seconds. 



(2.5) 



The analysis of the vibration of a body supported to vibrate only in the vertical 
or y direction can be carried out in a similar way to that above. Fig. 2.4 shows the 
system. 





+ve 



Fig. 2.4. Vertical motion. 



The spring extension 8 when the body is fastened to the spring is given by kd = mg. 
When the body is given an additional displacement >’ 0 and released the FBDs for a 
general displacement y, are as in Fig. 2.5. 



(a) 



kiy + 5) 

cb 



mg 



□ 



(b) 



Fig. 2.5. (a) Applied forces; (b) effective force. 
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The equation of motion is 
my = mg — k(y + (5) 

= mg — ky — kS = — ky, since mg = kS. 

that is, 

y + ( k/m)y — 0. (2.6) 

This is similar to equation (2.1), so that the general solution can be written as 

>' =, A cos y/(k/m)t + B sin y/(k/m)t. (2.7) 

Note that yj(k/m) = yjig/5), because mg = kd. That is, if <5 is known, then the frequency 
of vibration can be found. 

For the initial conditions y — y 0 at t = 0 and y = 0 at f = 0, 

y = >’o cos J(k/m)t. (2.8) 

Comparing (2.8) with (2.3) shows that for a given system the frequency of vibration 
is the same whether the body vibrates in a horizontal or vertical direction. 

Sometimes more than one spring acts in a vibrating system. The spring, which is 
considered to be an elastic element of constant stiffness, can take many forms in 
practice; for example, it may be a wire coil, rubber block, beam or air bag. Combined 
spring units can be replaced in the analysis by a single spring of equivalent stiffness 
as follows. 

{ 1 ) Springs connected in series 

The three-spring system of Fig. 2.6.(a) can be replaced by the equivalent spring of 
Fig. 2.6(b). 

F 

(b) 

Fig. 2.6. Spring systems. 

If the deflection at the free end, <5, experienced by applying the force F is to be the 
same in both cases, 

5 = F/k e = F/k l + F/k 2 + F/k 3 , 

i /*, = h/k. 





that is. 
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In general, the reciprocal of the equivalent stiffness of springs connected in series 
is obtained by summing the reciprocal of the stiffness of each spring. 

(2) Springs connected in parallel 

The three-spring system of Fig. 2.7(a) can be replaced by the equivalent spring of Fig. 
2.7(b). 




(a) (b) 



Fig. 2.7. Spring systems. 

Since the defection 5 must be the same in both cases, the sum of the forces exerted 
by the springs in parallel must equal the force exerted by the equivalent spring. Thus 

F = kyd + k 2 d + k 3 ( 5 = k e d, 

that is, 

= I K 

i= 1 

In general, the equivalent stiffness of springs connected in parallel is obtained by 
summing the stiffness of each spring. 

2.1.2 Torsional vibration 

Fig. 2.8 shows the model used to study torsional vibration. 

A body with mass moment of inertia / about the axis of rotation is fastened to a 
bar of torsional stiffness k T . If the body is rotated through an angle 9 a and released, 
torsional vibration of the body results. The mass moment of inertia of the shaft about 
the axis of rotation is usually negligible compared with /. 

For a general displacement 9 the FBDs are as given in Figs. 2.9(a) and (b). Hence 
the equation of motion is 

19= -kjd. 



e + = 0 . 



or 
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e 




i 



i 

Fig. 2.8. Single degree of freedom model - torsional vibration. 





Fig. 2.9. (a) Applied torque; (b) effective torque. 



This is of a similar form to equation (2.1). That is, the motion is simple harmonic 
with frequency \/2n y/(k r /I) Hz. 

The torsional stiffness of the shaft, k T , is equal to the applied torque divided by 
the angle of twist. 

Hence 



GJ 

k T = — , for a circular section shaft, 



where G = modulus of rigidity for shaft material, 

J — second moment of area about the axis of rotation, and 
l = length of shaft. 

Hence 



/= S = 2^ (GJ// ° HZ ’ 

and 



9 = 9 0 cos yf(GJ/Il)t, 
when 9 = 9 0 and 6 = 0 at t — 0. 

If the shaft does not have a constant diameter, it can be replaced analytically by 
an equivalent shaft of different length but with the same stiffness and a constant 
diameter. 
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For example, a circular section shaft comprising a length l { of diameter d Y and a 
length 1 2 of diameter d 2 can be replaced by a length l Y of diameter d x and a length 
l of diameter d t where, for the same stiffness, 

{GJ/l) length / 2 diameter d 2 (f**f/Olength //diameter d, 

that is, for the same shaft material, d 2 /l 2 — d^/l- 

Therefore the equivalent length l e of the shaft of constant diameter d x is given by 
/, = /,+ {djd 2 )%. 

It should be noted that the analysis techniques for translational and torsional 
vibration are very similar, as are the equations of motion. 

The torsional vibration of a geared system 

Consider the system shown in Fig. 2.10. The mass moments of inertia of the shafts 



Fixed 

end 




and gears about their axes of rotation are considered negligible. The shafts are 
supported in bearings which are not shown, and the gear ratio is IV: 1. 

From the FBDs, T 2 , the torque in shaft 2 is T 2 = k 2 (9 — (f>) — — M and T u the 
torque in shaft 1, is T, = k l N<f>; since = T 2 ,T 2 = and </> = k 2 0/(k 2 + k^ 2 ). 

Thus the equation of motion becomes 



10 + 



k 2 k 1 N z 
k 2 + k x N : 



0 = 0 , 



and 



/= W 7 Hz ' 

that is, k eq , the equivalent stiffness referred to shaft 2, is {k l k z N 2 )/(k 1 N 2 + k 2 ), or 
l/k eq = l/k 2 + l/ik^N 2 ). 
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2.1.3 Non-linear spring elements 

Any spring elements have a force-deflection relationship which is linear only over a 
limited range of deflection. Fig. 2.11 shows a typical characteristic. 

Force Spring 




Fig. 2.11. Non-linear spring characteristic. 

The non-linearities in this characteristic may be caused by physical effects such as 
the contacting of coils in a compressed coil spring, or by excessively straining the 
spring material so that yielding occurs. In some systems the spring elements do not 
act at the same time, as shown in Fig. 2.12(a), or the spring is designed to be non-linear 
as shown in Figs. 2.12(b) and (c). 





}r\ 

1 

-VWVAVvH 

-JVWWWV- 

-AWAW— 1 




i*^ 

HW/VWW 

-JWWWWV— 

1 — VWWWV 








Fig. 2.12. Non-linear spring systems. 
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Analysis of the motion of the system shown in Fig. 2.12(a) requires analysing the 
motion until the half-clearance a is taken up, and then using the displacement and 
velocity at this point as initial conditions for the ensuing motion when the extra 
springs are operating. Similar analysis is necessary when the body leaves the influence 
of the extra springs. See Example 4. 

2.1.4 Energy methods for analysis 

For undamped free vibration the total energy in the vibrating system is constant 
throughout the cycle. Therefore the maximum potential energy F mai is equal to the 
maximum kinetic energy T max although these maxima occur at different times during 
the cycle of vibration. Furthermore, since the total energy is constant, 

T + V — constant, 

and thus 

A (T + V) = 0. 

dr 

Applying this method to the case, already considered, of a body of mass m fastened 
to a spring of stiffness k, when the body is displaced a distance x from its equilibrium 
position, 

strain energy (SE) in spring = \ kx 2 . 
kinetic energy (KE) of body = \ mx 2 . 

Hence 

V =1 kx 2 , 

and 

T = \mx 2 . 

Thus 

y- (|mx 2 -I- \kx 2 ) = 0, 

dr 



that is 



mxx + kxx — 0, 



or 



x + 




|x = 0, as before, equation (2.1). 



This is a very useful method for certain types of problem in which it is difficult to 
apply Newton’s laws of motion. 

Alternatively, assuming SHM, if x — x 0 cos cot. 
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The maximum SE, V max = ^ kxl 
and the maximum KE, T max = \m(x 0 a>) 2 ' 

Thus, since T max = V max , 

\ kx o = { mxg to 2 , 
or co — N /( k/m ) rad/s. 

Example 1 

A link AB in a mechanism is a rigid bar of uniform section 0.3 m long. It has a mass 
of 10 kg, and a concentrated mass of 7 kg is attached at B. The link is hinged at A 
and is supported in a horizontal position by a spring attached at the mid point of 
the bar. The stiffness of the spring is 2 kN/m. Find the frequency of small free 
oscillations of the system. The system is as shown below. 




For rotation about A the equation of motion is 
I J) = —ka 2 6, 

that is, 

8 + ( ka 2 /I A )d = 0. 

This is SHM with frequency 

^(ka 2 /I A ) Hz. 

In this case 

a = 0.15m, l — 0.3m, k = 2000N/m, 
and 



/ A = 7(0.3) 2 + i x 10(0.3) 2 = 0.93 kg m 2 . 
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Hence 



J_ If 200Q(0.1S) 2 \ _ 

J 2 TtVV 0-93 J ' 



Hz. 



Example 2 

A small turbo-generator has a turbine disc of mass 20 kg and radius of gyration 
0.15 m driving an armature of mass 30 kg and radius of gyration 0.1 m through a steel 
shaft 0.05 m diameter and 0.4 m long. The modulus of rigidity for the shaft steel is 
86 x 10 9 N/m 2 . Determine the natural frequency of torsional oscillation of the system, 
and the position of the node. The shaft is supported by bearings which are not shown. 




The rotors must twist in opposite directions to each other; that is, along the shaft 
there is a section of zero twist: this is called a node. The frequency of oscillation of 
each rotor is the same, and since there is no twist at the node, 




that is, 






Now 



I A = 20(0.15) 2 — 0.45 kg m 2 , and I B = 30(0.1) 2 — 0.3kg m 2 . 

Thus 



U = (0.3/0.45)/ b . 

Since l A + 1 B = 0.4 m, l A — 0.16m, 

that is, the node is 0.16 m from the turbine disc. 
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Hence 



f= 



1 // 86 x 10 9 7r(0.05) 4 

2WV °- 45 x 0.16 X 32 



136 Hz ( = 136 x 60 = 8160 rev/min). 



If the generator is run near to 8160 rev/min this resonance will be excited causing 
high dynamic stresses and probable fatigue failure of the shaft at the node. 



Example 3 

A uniform cylinder of mass m is rotated through a small angle 9 0 from the equilibrium 
position and released. Determine the equation of motion and hence obtain the 
frequency of free vibration. The cylinder rolls without slipping. 




If the axis of the cylinder moves a distance x and turns through an angle 9 so that 
x = r9, 

KE = \ mx 2 + \l8 2 , where / = | mr 2 . 

Hence 

KE = | mr 2 8 2 . 

SE = 2 x \ x k[(r + a)9Y = k(r + a) 2 9 2 - 
Now, energy is conserved, so (f mr 2 & 2 + k(r + a) 2 9 2 ) is constant, that is, 

(| mr 2 9 2 + k(r + a) 2 9 2 ) = 0 
dt 



or 



|mr 2 2M + k(r + a) 2 29& - 0. 



Thus the equation of the motion is 
k(r 4- a) 2 9 



8 + 



o. 



(3/4)mr 2 

Hence the frequency of free vibration is 



1 

2n> 



4 k(r + a) 2 
3mr 2 



Hz. 
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Example 4 

In the system shown in Fig. 2.12(a) find the period of free vibration of the body if it 
is displaced a distance x 0 from the equilibrium position and released, when x 0 is 
greater than the half clearance a. Each spring has a stiffness k, and damping is 
negligible. 

Initially the body moves under the action of four springs so that the equation of 
motion is 

m'x = — 2 kx — 2k(x — a), 

that is, 

m'x + 4 kx — 2ka. 

The general solution comprises a complementary function and a particular integral 
such that 



x — A cos cot + B sin cot + -. 

2 



With initial conditions x = x 0 at t = 0, and x = 0 at t = 0, the solution is 
x = (x 0 — ^ cos cot + where co = I (— ^ rad/s. 



m 



When x = a, t = 1 1 and 

a = ^x 0 — ^ cos COty 

so that 

1 1 = 



— | cos 
4k 



2(x 0 - a/2) 



This is the time taken for the body to move from x = x 0 to x = a. 
If a particular value of x 0 is chosen, x 0 = 2a say, then 



3a 



x = — cos cot + ~ and t x = /( — ) cos 



m 



1 

3 f 



and when x = a, 

3 a 

x = — —to sin coty = — ay/(8k/m). 

For motion from x = a to x = 0 the body moves under the action of two springs 
only, with the initial conditions 

, ,, 8k 

x n = a and x = — a / — 



The equation of motion is m'x + 2 kx = 0 for this interval, the general solution to 
which is x — C cos Q t + D sin Qf, where Q = s J(2kjm) rad/s. 
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Since x 0 — a at t = 0, C = a 
and because 




and 



x — — C(l sin (It + D(i cos (It 
then D — —2a. That is, 

x = a cos Qr — 2 a sin (It. 

When x = 0, t = t 2 and cos (lt 2 = 2 sin (it 2 or tan (it 2 = that is 

f -=V(s) ,an_i G} 

Thus the time for \ cycle is (t, +r 2 ), and the time for one cycle, that is the period 
of free vibration when x 0 = 2a \§ tyt^ + t 2 ) or 

{VSMWSK# 

Example 5 

A uniform wheel of radius R can roll without slipping on an inclined plane. Concentric 
with the wheel, and fixed to it, is a drum of radius r around which is wrapped one 
end of a string. The other end of the string is fastened to an anchored spring, of 
stiffness k, as shown. Both spring and string are parallel to the plane. The total mass 
of the wheel/drum assembly is m and its moment of inertia about the axis through 
the centre of the wheel O is /. If the wheel is displaced a small distance from its 
equilibrium position and released, derive the equation describing the ensuing motion 
and hence calculate the frequency of the oscillations. Damping is negligible. 
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If the wheel is given an anti-clockwise rotation 6 from the equilibrium position, 
the spring extension is (R + r) 9 so that the restoring spring force is k(R + r)6. 

The FBDs are 





The rotation is instantaneously about the contact point A so that taking moments 
about A gives the equation of motion as 

I a 8 = -k(R + r) 2 6. 

(The moment due to the weight cancels with the moment due to the initial spring 
tension.) 

Now I A = I + mR 2 , so 



0 + 



k(R + rf 



0 - 0 , 



I + mR 2 

and the frequency of oscillation is 



1 

2n\ 



k(R + r)‘ 
1 + mR 2 



Hz. 



An alternative method for obtaining the frequency of oscillation is to consider the 
energy in the system. 



Now 



SE, V= ±k(R + r) 2 Q 2 , 



and 



KE, T= \l A t) 2 . 

(Weight and initial spring tension effects cancel.) 
So T + V= |/ A 0 2 + \k{R + r) 2 8 2 . 
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Hence 

I A 0 + k(R + r) 2 G = 0, 

which is the equation of motion. 

Or, we can put V max = T max , and if 9 = 6 0 sin cot is assumed, 

lk(R + r) 2 e 2 0 = i/ A co 2 9 2 0 , 

so that 




where 



/ A = / + mr 2 and / = {w/ln) Hz. 



Example 6 

A uniform building of height 2 h and mass m has a rectangular base a x b which rests 
on an elasic soil. The stiffness of the soil, k, is expressed as the force per unit area 
required to produce unit deflection. 

Find the lowest frequency of free low-amplitude swaying oscillation of the building. 
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The lowest frequency of oscillation about the axis 0-0 through the base of the 
building is when the oscillation occurs about the shortest side, length a. 

I 0 is the mass moment of inertia of the building about axis 0-0. 

The FBDs are 




and the equation of motion for small 9 is given by 
I 0 9 - mghO — M, 



where M is the restoring moment from the elastic soil. 

For the soil, k = force/(area x deflection), so considering an element of the base 
as shown, the force on element = kb dx x xO, and the moment of this force about 
axis 0-0 = kbdx x x6 x x. Thus the total restoring moment M, assuming the soil 
acts similarly in tension and compression is 



M = 2 



'aj 2 



kbx 2 9dx 



= 2kb9 = k ^2 e . 
3 12 



Thus the equation of motion becomes 
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Motion is therefore simple harmonic, with frequency 

An alternative solution can be obtained by considering the energy in the system. In 
this case, 

T = i/ o 0 2 , 

and 



K = i2 [' 



a/2 ., mghO 2 

kbdx x xQ x x6 — , 



where the loss in potential energy of building weight is given by mgh (1 — cos 6) ~ 
mgh0 2 /2, since cos 6 ~ 1 — 0 2 /2 for small values of 6. Thus 



V = 



ka 3 b 



- • 
24 2 ) 



Assuming simple harmonic motion, and putting T mstx = F max , gives 
, fka 3 b/\2 — mgh\ 

" = ( — r. — ) 

as before. 

Note that for stable oscillation, a) > 0, so that 





> 0 . 



That is ka 3 b > \2mgh. 

This expression gives the minimum value of k, the soil stiffness, for stable oscillation 
of a particular building to occur. If k is less that 12 mgh/a 3 b the building will fall 
over when disturbed. 



2.2 FREE DAMPED VIBRATION 

All real systems dissipate energy when they vibrate. The energy dissipated is often 
very small, so that an undamped analysis is sometimes realistic; but when the damping 
is significant its effect must be included in the analysis, particularly when the amplitude 
of vibration is required. Energy is dissipated by frictional effects, for example that 
occurring at the connection between elements, internal friction in deformed members, 
and windage. It is often difficult to model damping exactly because many mechanisms 
may be operating in a system. However, each type of damping can be analysed, and 
since in many dynamic systems one form of damping predominates, a reasonably 
accurate analysis is usually possible. 

The most common types of damping are viscous, dry friction and hysteretic. 
Hysteretic damping arises in structural elements due to hysteresis losses in the material. 
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The type and amount of damping in a system has a large effect on the dynamic 
response levels. 

2.2.1 Vibration with viscous damping 

Viscous damping is a common form of damping which is found in many engineering 
systems such as instruments and shock absorbers. The viscous damping force is 
proportional to the first power of the velocity across the damper, and it always 
opposes the motion, so that the damping force is a linear continuous function of the 
velocity. Because the analysis of viscous damping leads to the simplest mathematical 
treatment, analysts sometimes approximate more complex types of damping to the 
viscous type. 

Consider the single degree of freedom model with viscous damping shown in Fig. 
2.13. 




Fig. 2.13. Single degree of freedom model with viscous damping. 



The only unfamiliar element in the system is the viscous damper with coefficient 
c. This coefficient is such that the damping force required to move the body with a 
velocity x is cx. 

For motion of the body in the direction shown, the free body diagrams are as in 



(a) (b) 

Fig. 2.14. (a) Applied force; (b) effective force. 

Fig. 2.14(a) and (b). The equation of motion is therefore 

mx + cx + kx = 0. (2.9) 

This equation of motion pertains to the whole of the cycle: the reader should verify 
that this is so. (Note: displacements to the left of the equilibrium position are negative, 
and velocities and accelerations from right to left are also negative). 
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Equation (2.9) is a second-order differential equation which can be solved by 
assuming a solution of the form x = Xe sl . Substituting this solution into equation 
(2.9) gives 

(ms 2 + cs + k)Xe s ' = 0. 

Since Ae st ^ 0 (otherwise no motion), ms 2 + cs + k = 0, 



so, 



c J(c 2 — 4 mk) 
Su2= ~2m ± ~~2m ' 



Hence 



x = + X 2 e v , 

where X 1 and X 2 are arbitrary constants found from the initial conditions. The 
system response evidently depends on whether c is positive or negative, and on 
whether c 2 is greater than, equal to, or less than 4mk. 

The dynamic behaviour of the system depends on the numerical value of the radical, 
so we define critical damping as that value of c(c c ) which makes the radical zero: that is, 

c c = 2 yf(km). 

Hence 



cJ2m = yf(k/m) = co, the undamped natural frequency, and 
c c = 2 yj(km) = 2mo). 

The actual damping in a system can be specified in terms of c c by introducing the 
damping ratio £. 

Thus £ = c/c c , 

and s 12 = [ - £ ± yf(£ 2 - 1 )]a>. (2.10) 

The response evidently depends on whether c is positive or negative, and on whether 
£ is greater than, equal to, or less than unity. Usually c is positive, so we need consider 
only the other possibilities. 

Case 1. f < 1 ; that is, damping less than critical 
From equation (2.10) 

•Si . 2 = - ±j v / ( 1 ~ C 2 )“> where j = !)> 

so 

x = e -c “‘[A’ 1 e jv<1 + A r 2 e“ j ' <1_i ' )w '] 

and 

x = A'e - ^' sin (^/(l — £ 2 ) a>t + (p). 
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The motion of the body is therefore an exponentially decaying harmonic oscillation 
with circular frequency co v = co^/fl — £ 2 ), as shown in Fig. 2.15. 

The frequency of the viscously damped oscillation co v , is given by co v = a\/(l — £ 2 ), 
that is, the frequency of oscillation is reduced by the damping action. However, in 
many systems this reduction is likely to be small, because very small values of £ are 
common; for example in most engineering structures £ is rarely greater than 0.02. 
Even if £ = 0.2, co v = 0.98 co. This is not true in those cases where £ is large, for 
example in motor vehicles where £ is typically 0.7 for new shock absorbers. 



Case 2. £ = 1; that is, critical damping 

Both values of s are — co. However, two constants are required in the solution of 
equation (2.9), thus x = (A + Br)e~“' may be assumed. 

Critical damping represents the limit of periodic motion, hence the displaced body 
is restored to equilibrium in the shortest possible time, and without oscillation or 
overshoot. Many devices, particularly electrical instruments, are critically damped to 
take advantage of this property. 

Case 3. £ > 1; that is, damping greater than critical 
There are two real values of s, so x = A^e 51 ' 4- X 2 e s -'. 

Since both values of s are negative the motion is the sum of two exponential decays, 
as shown in Fig. 2.16. 

Logarithmic decrement A 

A convenient way of determining the damping in a system is to measure the rate of 
decay of oscillation. It is usually not satisfactory to measure co v and co, because unless 
£ > 0.2, co i co v . 

The logarithmic decrement, A, is the natural logarithm of the ratio of any two 
successive amplitudes in the same direction, and so from Fig. 2.17. 
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Thus 



Xq ~ 

A = ln Xe = ^ T v 



Since 



2n 



2n 



co v a)J{\ -C 2 )’ 



A = 



2tcC 



= In 



Vd- n v*n 



*i 



For small values of £(> 0.25), A ~ 2nC- 

It should be noted that this analysis assumes that the point of maximum 
displacement in a cycle and the point where the envelope of the decay curve Ae~ ? "' 
touches the decay curve itself, are coincident. This is usually very nearly so, and the 
error in making this assumption is usually negligible, except in those cases where the 
damping is high. 

For low damping it is preferable to measure the amplitude of oscillations many 
cycles apart so that an easily measurable difference exists. 



In this case A = In 




N 



x In 




since 







etc. 



Example 7 

Consider the transverse vibration of a bridge structure. For the fundamental frequency 
it can be considered as a single degree of freedom system. The bridge is deflected at 
mid-span (by winching the bridge down) and suddenly released. After the initial 
disturbance the vibration was found to decay exponentially from an amplitude of 
10 mm to 5.8 mm in three cycles with a frequency of 1.62 Hz. The test was repeated 
with a vehicle of mass 40 000 kg at mid-span, and the frequency of free vibration 
was measured to be 1.54 Hz. 

Find the effective mass, the effective stiffness, and the damping ratio of the structure. 



Let m be the effective mass and k the effective stiffness. Then 



and 



fi = 1-62 = 



f 2 = 1.54 = 



1 



m 



Hz, 



Hz, 



2n\] \m + 40 x 10 3 / 
if it is assumed that £ is small enough for/ v ~f 



Thus 



1.62\ 2 m + 40x 10 3 
L54 ) ~ m ' 
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so 



m — 375 x 10 3 kg. 

Since 

k = (Infrfm, 
k = 38 850 kN/m. 

Now 



A U 

= 0.182. 



X 




Hence 



A = 



V(1 - C 2 ) 



0.182, 



and so £ = 0.029. (This compares with a value of about 0.002 for cast iron material. 
The additional damping originates mainly in the joints of the structure.) This value 
of £ confirms the assumption that f y ~f 



Example 8 

A light rigid rod of length L is pinned at one end O and has a body of mass m 
attached at the other end. A spring and viscous damper connected in parallel are 
fastened to the rod at a distance a from the support. The system is set up in a 
horizontal plane: a plan view is shown. 

Assuming that the damper is adjusted to provide critical damping, obtain the 
motion of the rod as a function of time if it is rotated through a small angle 8 0 and 
then released. Given that 6 0 = 2° and the undamped natural frequency of the system 
is 3 rad/s, calculate the displacement 1 s after release. 

Explain the term logarithmic decrement as applied to such a system and calculate 
its value assuming that the damping is reduced to 80% of its critical value. 
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Let the rod turn through an angle 9 from the equilibrium position. 

Note that the system oscillates in the horizontal plane so that the FBDs are: 



Spring force = kaQ 
Damper force = ca8 




Applied forces Effective forces 



Taking moments about the pivot O gives 
I 0 8 — — ca 2 9 — ka 2 9, 

where I 0 = ml 2 , so the equation of motion is 
ml}Q + ca 2 0 + ka 2 6 — 0. 

Now the system is adjusted for critical damping, so that £ = 1. 
The solution to the equation of motion is therefore of the form 

6 = (A + Br)e~“'. 

Now, 6 = 9 0 when t = 0, and d 9/dt = 0 when t — 0. Hence 

0 o = A, 

and 

0 = Be ~ + (A + Bt){- ce)e~"', 

so that 

B = 9 0 to. 

Hence 

9 = 8 0 {\ + 

If co = 3 rad/s, t = 1 s and 9 0 = 2°, 

9 = 2(1 + 3)e^ 3 = 0.4°. 

The logarithmic decrement 
A . X x 

A = In — — = z-, 

-A 1 1 v ( 1 — s ) 
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so that if C = 0.8, 



A = 



5.027 

lX6~ 



8.38 



Root locus study of damping 

It is often convenient to consider how the roots of equation (2.10) vary as ( increases 
from zero. The roots of this equation are given by 

s i. 2 /°> = - C+jsA 1 - £ 2 ) for 1 > C > 0 

and 

Sm/to = ~ C±yJ(C 2 - 1) for c > 1. 

These roots can be conveniently displayed in a plot of imaginary (s/co) against real 
( s/co ); since for every value of f there are two values of (s/co), the roots when plotted 
form two loci as shown in Fig. 2.18. The position of a root in the (s/co) plane indicates 
the frequency of oscillation (Im (s/co) axis) of the system, if any, and the rate of growth 
of decay of oscillation (Re (s/co) axis). 

Because the Re (s/co) is negative on the left of the Im (s/co) axis, all roots which lie 
to the left of the Im (s/co) axis represent a decaying oscillation and therefore a stable 
system. Roots to the right of the Im (s/co) axis represent a growing oscillation and 
an unstable system. 

This root study of damping is a useful design technique, because the effects of 
changing the damping ratio on the response of a system can easily be seen. It also 
has important applications in the study of control system dynamics. Note that the 
damping ratio C is given by cos 9 in Fig. 2.18. 




Fig. 2.18. Root locus plot. 
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2.2.2 Vibration with Coulomb (dry friction) damping 

Steady friction forces occur in many systems when relative motion takes place between 
adjacent members. These forces are independent of amplitude and frequency; they 
always oppose the motion and their magnitude may, to a first approximation, be 
considered constant. Dry friction can, of course, just be one of the damping 
mechanisms present; however, in some systems it is the main source of damping. In 
these cases the damping can be modelled as in Fig, 2.19 




Fig. 2.19. System with Coulomb damping. 



The constant friction force F d always opposes the motion, so that if the body is 
displaced a distance x 0 to the right and released from rest we have, for motion from 
right to left only, 

mx = F d — kx 

or mx + kx — F d . (2.11) 

The solution to the complementary function is x = A sin cot + B cos cot, and the 
complete solution is 

p 

x — A sin cot + B cos cot + ~r (2-12) 

k 

where to - y/(k/m) rad/s. 

Note. The particular integral may be found by using the D-operator. Thus equation 
(2.11) is 

(D 2 + co 2 )x = FJm 



so 



x = (l/co 2 )[l + (D> 2 )]- 1 F d /m 
= [1 - (D 2 /<u 2 ) + •••]F d /moj 2 = FJk. 

The initial conditions were x = x 0 at t = 0, and x = 0 at t = 0. Substitution into 
equation (2.12) gives 
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p 

A = 0 and B — x 0 — ~. 

k 

Hence 




At the end of the half cycle right to left, <x>t — it and 



x 



(t - n/co) 




(2.13) 



That is, there is a reduction in amplitude of 2F d /k per half cycle. 

From symmetry, for motion from left to right when the friction force acts in the 
opposite direction to the above, the initial displacement is (x 0 — 2FJk) and the final 
displacement is therefore (x 0 — 4 FJk), that is the reduction in amplitude is 4 FJk per 
cycle. This oscillation continues until the amplitude of the motion is so small that 
the maximum spring force is unable to overcome the friction force F d . This can 
happen whenever the amplitude is ^ + (FJk). The manner of oscillation decay is 
shown in Fig. 2.20; the motion is sinusoidal for each half cycle, with successive half 
cycles centred on points distant + (FJk) and — (FJk) from the origin. The oscillation 
ceases with | x | ^ FJk. The zone x = + FJk is called the dead zone. 

To determine the frequency of oscillation we rewrite the equation of motion (2.1 1) as 

mx + k(x ~ (FJk)) = 0. 

Now if x' = x — (FJk), x' = x so that mx' + kx' = 0, from which the frequency of 
oscillation is (l/27t) J(k/m) Hz. That is, the frequency of oscillation is not affected by 
Coulomb friction. 



Example 9 

Part of a structure can be modelled as a torsional system comprising a bar of stiffness 
10 kN m/rad and a beam of moment of inertia about the axis of rotation of 
50 kg m 2 . The bottom guide imposes a friction torque of 10 N m. 

Given that the beam is displaced through 0.05 rad from its equilibrium position 
and released, find the frequency of the oscillation, the number of cycles executed 
before the beam motion ceases, and the position of the beam when this happens. 
Now 



co = 




/ 10 x 10 3 
50 



14.14 rad/s. 



14.14 

2 n 



2.25 Hz. 



Hence 





Fig. 2.20. Vibration decay of system with Coulomb damping 
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Loss in amplitude/cycle 



4£a 

k 



4 x 10 
10 * 



rad 



= 0.004 rad. 



Number of cycles for motion to cease 



0.05 

0.004 



= 12i 



The beam is in the initial (equilibrium) position when motion ceases. The motion 
is shown opposite. 



2.2.3 Vibration with combined viscous and Coulomb damping 

The free vibration of dynamic systems with viscous damping is characterized by an 
exponential decay of the oscillation, whereas systems with Coulomb damping possess 
a linear decay of oscillation. Many real systems have both forms of damping, so that 
their vibration decay is a combination of exponential and linear functions. 

The two damping actions are sometimes amplitude dependent, so that initially the 
decay is exponential, say, and only towards the end of the oscillation does the Coulomb 
effect show. In the analyses of these cases the Coulomb effect can easily be separated 
from the total damping to leave the viscous damping alone. The exponential decay 
with viscous damping can be checked by plotting the amplitudes on logarithmic-linear 
axes when the decay should be seen to be linear. 
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If the Coulomb and viscous effects cannot be separated in this way, a mixture of 
linear and exponential decay functions have to be found by trial and error in order 
to conform with the experimental data. 



2.2.4 Vibration with hysteretic damping 

Experiments on the damping that occurs in solid materials and structures which have 
been subjected to cyclic stressing have shown the damping force to be independent 
of frequency. This internal, or material, damping is referred to as hysteretic damping. 
Since the viscous damping force cx is dependent on the frequency of oscillation, it 
is not a suitable way of modelling the internal damping of solids and structures. The 
analysis of systems and structures with this form of damping therefore requires the 
damping force cx to be divided by the frequency of oscillation a>. Thus the equation 
of motion becomes mx + (e/a>)x + kx = 0. 

However, it has been observed from experiments carried out on many materials 
and structures that under harmonic forcing the stress leads the strain by a constant 
angle, x. 

Thus for an harmonic strain, e — s 0 sin vt, where v is the forcing frequency, the 
induced stress is a = o 0 sin (vr + x). Hence 

a — a 0 cos x sin vt + a 0 sin x cos vt 
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= cr 0 cos a sin vf + a 0 sin a sin ^vf + -J. 

The first component of stress is in phase with the strain e, whilst the second 
component is in quadrature with e and n/ 2 ahead. Putting j = ^( — 1), 

a — c r 0 cos a sin vf + ja 0 sin a sin vf. 

Hence a complex modulus £* can be formulated, where 

(X 0" q , (T o . 

£* = - = — cos a + j — sin a 
e e 0 £ 0 

= F+;£\ 

where £' is the in-phase or storage modulus, and E" is the quadrature or loss modulus. 

The loss factor rj, which is a measure of the hysteretic damping in a structure, is 
equal to E"/E', that is, tan a. 

It is not usually possible to separate the stiffness of a structure from its hysteretic 
damping, so that in a mathematical model these quantities have to be considered 
together. The complex stiffness k* is given by k* = k(l + j rj), where k is the static 
stiffness and rj the hysteretic damping loss factor. 

The equation of free motion for a single degree of freedom system with hysteretic 
damping is therefore nix + k*x — 0. Fig. 2.21 shows a single degree of freedom model 
with hysteretic damping of coefficient c H . 




Fig. 2.21. Single degree of freedom model with hysteretic damping. 



The equation of motion is 
mx + (c H /a»)x + kx = 0. 

Now if x = A r e j "', 

f c h \ 

x - j (ox and ( — I x = jc H x. 

Thus the equation of motion becomes 
mx + (k + jc H )x = 0. 
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Since 

J C H 



we can write 

mx + k*x = 0. 

That is, the combined effect of the elastic and hysteretic resistance to motion can be 
represented as a complex stiffness, k*. 

A range of values of r\ for some common engineering materials is given below. For 
more detailed information on material damping mechanisms and loss factors, see 
Damping of Materials and Members in Structural Mechanics by B. J. Lazan (Pergamon, 
1968). 



T = KU +JW 



k + jc H = k 1 + 



Material 


Loss factor 


Aluminium-pure 


0.00002-0.002 


Aluminium alloy-dural 


0.0004-0.001 


Steel 


0.001-0.008 


Lead 


0.008-0.014 


Cast iron 


0.003-0.03 


Manganese copper alloy 


0.05-0.1 


Rubber-natural 


0. 1-0.3 


Rubber-hard 


1.0 


Glass 


0.0006-0.002 


Concrete 


0.01-0.06 



2.2.5 Energy dissipated by damping 

The energy dissipated per cycle by the viscous damping force in a single degree of 
freedom vibrating system is approximately 
fx 

4 cx d.x. 

Jo 



if x = X sin cor is assumed for the complete cycle. The energy dissipated is therefore 



4 



f*/2 

cX 2 io 2 cos cot dr = nccoX 2 . 



Jo 

The energy dissipated per cycle by Coulomb damping is 4 F d X approximately. Thus 
an equivalent viscous damping coefficient for Coulomb damping c d can be deduced, 
where 



7 zc d coX 2 = 4 F d X, 




44 The vibrations of systems having one degree of freedom 



[Ch. 2 



that is, 

c 

d TtCoX 

The energy dissipated per cycle by a force F acting on a system with hysteretic 
damping is Jf dx, where F = k*x = k{ 1 + j>/)x, and x is the displacement. 

For harmonic motion x = X sin cot. 



so 



F = kX sin cot + j qkX sin cot 
= kX sin cot + rjkX cos cot. 



Now 



sin cot = 



x 

X' 



therefore 



cos cot 



J(X 2 - x 2 ) 
X 



Thus 

F = kx ± t]kyj{X 2 — x 2 ). 

This is the equation of an ellipse as shown in Fig. 2.22. The energy dissipated is 
given by the area enclosed by the ellipse. 




Hence 




(fcx ± r}ksJ(X 2 — x 2 ))dx 



= nX 2 rjk. 
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An equivalent viscous damping coefficient c H is given by 
nc H coX 2 = nt]kX 2 , 

that is 




Note also that c — c H /a>. 

Example 10 

A single degree of freedom system has viscous damping, with ( = 0.02. Find the energy 
dissipated per cycle as a function of the energy in the system at the start of that cycle. 
Also find the amplitude of the 12th cycle if the amplitude of the 3rd cycle is 1.8 mm. 

c 1, so In {XJX 2 ) = 2nC = 0.126. 

Thus 

XJX 2 = e 0 126 = 1.134. 

Energy at start of cycle = jkX 2 (stored as strain energy in spring) 

Energy at end of cycle = jkX 2 2 

Energy dissipated during cycle 
Energy at start of cycle 
I by 2 _ it y 2 

= - — tttA — - = 1 - ( x 2 /x i ) 2 - 1 - 0.7773 = 0.223. 

that is 22.3% of initial energy is dissipated in one cycle. 

Now 

XJX 2 = 1.134, XJX 3 = 1.134... (X n _ l )/X n = 1.134. 

Therefore 

X 3 /X 12 = (1.134) 9 = 3.107 

that is 

= 3157 = °' 579 



Example 11 

A gun is designed so that when fired the barrel recoils against a spring. At the end 
of the recoil a viscous damper is engaged which allows the barrel to return to its 
equilibrium position in the minimum time without overshoot. Determine the spring 
constant and damping coefficient for a gun whose barrel has a mass of 500 kg if the 
initial recoil velocity is 30 m/s and the recoil distance is 1.6 m. 
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Find also the time required for the barrel to return to a position 0.1 m from its 
initial position. 

Since the initial kinetic energy of the barrel must equal the maximum strain energy 
in the spring, 

|fcl.6 2 = |500(30) 2 , 

that is k = 175 kN/m. 

Since critical damping is required, the damping coefficient is 
2^(500 x 175 x 10 3 ) = 18 700 Ns/m. 

For critical damping 
x = (A + Bt)e ~ 

If x = x 0 and x = 0 at t = 0, 

-x = x 0 (l + cof)e“°". 

Now oj = ^/(175 x 10 3 /500) = 18.7 rad/s and x = 0.1 m, so we require the solution 
to 0.1 = 1.6(1 + 18.7f)e“ 18 - 7 '. 

This equation can be solved in many ways, e.g. by trial and error or graphically, 
to give t = 0.24 s. 

2.3 FORCED VIBRATION 

Many real systems are subjected to periodic excitation. This may be due to unbalanced 
rotating parts, reciprocating components, or a shaking foundation. Sometimes large 
motions of the suspended body are desired as in vibratory feeders and compactors, 
but usually we require very low vibration amplitudes over a large range of exciting 
forces and frequencies. Some periodic forces are harmonic, but even if they are not, 
they can be represented as a series of harmonic functions using Fourier analysis 
techniques. Because of this the response of elastically supported bodies to harmonic 
exciting forces and motions must be studied. 

2.3.1 Response of a viscous damped system to a simple harmonic exciting force with 
constant amplitude 

In the system shown in Fig. 2.23, the body of mass m is connected by a spring and 
viscous damper to a fixed support, whilst an harmonic force. of circular frequency v 
and amplitude F acts upon it, in the line of motion. 

The equation of motion is 

m'x + cx + kx = F sin vf. (2-14) 

The solution to m'x + cx + kx = 0, which has already been studied, is the 
complementary function; it represents the initial vibration which quickly dies away. 
The sustained motion is given by the particular solution. A solution x = X (sin vr — </>) 
can be assumed, because this represents simple harmonic motion at the frequency of 
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Fig. 2.23. Single degree of freedom model of a forced system with viscous damping. 



the exciting force with a displacement vector which lags the force vector by </>, that 
is the motion occurs after the application of the force. 

Assuming x = X sin(vf — (ft), 

x = Xv cos (vr — 0) = Xv sin (vf — 4> + \n), 

and 

x = — Xv 2 sin (vf — (j>) — Xv 2 sin (vf — <f> + n). 

The equation of motion (2.14) thus becomes 

mXv 2 sin (vf — (j) + n) + cXv sin (vf — cp + n/2) + kX sin (vf — (f>) 

= F sin vf. 

A vector diagram of these forces can now be drawn, Fig 2.24. 




From the diagram, 

F 2 = (kX - mXv 2 ) 2 + (cX\) 2 , 

or 

X = F/J((k - mv 2 ) 2 + (c v) 2 ), (2.15) 
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and 

tan <fi = cXv/(kX — mXv 2 ). 

Thus the steady-state solution to equation (2.14) is 
F 

X = — 7771 m 7 — ( vt ~ $)> 

yj((k - mv 2 ) 2 + (cv) 2 ) 

where 



(f> = tan" 



cv 



k — mv 2 / 



The complete solution includes the transient motion given by the complementary 
function: 

x = /4e -? “' sin (o\/(l — £ 2 )t -I- a). 

Fig. 2. 25 shows the combined motion. 




Equation (2.15) can be written in a more convenient form if we put 



co = 




rad/s 



and 



X . = 



F 

k' 




Then 



X 

X 



(2.16) 
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and (j) — tan 



' 2c (v/to) 

1 - (v/cu) 2 



X/X s is known as the dynamic magnification factor, because X s is the static 
deflection of the system under a steady force F, and X is the dynamic amplitude. 

By considering different values of the frequency ratio v/cn, we can plot X/X s and 
<fi as functions of frequency for various values of £. Figs. 2.26 and 2.27 show the results. 




Fig. 2.27. Phase-frequency response for system of Fig. 2.23. 
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Fig. 2.28. Forced vibration vector diagrams: (a) v/w « 1 Exciting force approximately equal to spring 
force; (b) v/cj = 1 Exciting force equal to damping force, and inertia force equal to spring force; (c) v/w » 1 
Exciting force nearly equal to inertia force. 



The effect of the frequency ratio on the vector diagram is shown in Fig. 2.28. 

The importance of mechanical vibration arises mainly from the large values of 
X/X s experienced in practice when v/co has a value near unity: this means that a 
small harmonic force can produce a large amplitude of vibration. The phenomenon 
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known as resonance occurs when the forcing frequency is equal to the natural 
frequency, that is when v/co = 1. The maximum value of X/X s actually occurs at 
values of v/co less than unity: the value can be found by differentiating equation (2.16) 
with respect to v/co. Hence 

(v/<o\ XIXs)max = V(! - 2 £ 2 ) - 1. for C small > 



and 



(X/XJ max = i/(2CVd - C 2 )). 

For small values of (, (A'/X s ) max ~ 1/2C which is the value pertaining to v/co = 1. 
l/2£ is a measure of the damping in a system and is referred to as the Q factor. 
Note. An alternative solution to the equation of motion can be obtained by putting 
F sin vt = Im(Fe jw ). 

Then mx + cx + kx = Fe ,vt , and a solution x = A r e jv< can be assumed. 

F 

Thus ( k — mv 2 )X + j cvX = F, or X = — T : — . 

(k — mv ) + jcv 

F 

Hence X — -j— — — — - 5 -, as before, equation (2.15). 

>/(( fc “ mv ) + ( cv ) ) 

Both reciprocating and rotating unbalance in a system produce an exciting force of 
the inertia type and result in the amplitude of the exciting force being proportional 
to the square of the frequency of excitation. 

For an unbalanced body of mass m r , at an effective radius r, rotating at an angular 
speed v, the exciting force is therefore m/v 2 . If this force is applied to a single degree 
of freedom system such as that in Fig. 2.23, the component of the force in the direction 
of motion is m/v 2 sin vt, and the amplitude of vibration is 

_ ( m r /m)rjv/co ) 2 

V((l ~ {v/co) 2 ) 2 + (2£vco) 2 ) 1 j 

(see equation (2.16)). 

The value of v/co for maximum X found by differentiating equation (2.17) is given by 
( v Mrma* = WO - 2 C 2 ) 

that is, the peak of the response curve occurs when v > co. This is shown in Fig. 2.29 
Also, 



X mai = (m»r/ 2fV(l - C 2 ). 
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It can be seen that away from the resonance condition (v/eo = 1) the system response 
is not greatly affected by damping unless this happens to be large. Since in most 
mechanical systems the damping is small (( < 0.1) it is often permissible to neglect 
the damping when evaluating the frequency for maximum amplitude and also the 
amplitude-frequency response away from the resonance condition. 




Example 12 

A turbo-charger rotor consists of a steel shaft 10 mm in diameter which is simply 
supported in two self-aligning bearings placed 100 mm apart. A turbine rotor of mass 
0.2 kg is located at the mid-point of the shaft. There is negligible damping. The rotor 
is carefully balanced, the maximum out-of-balance being 10“ 5 kg m. In previous 
applications the turbo-charger has been run successfully at 35 000 rev/min and 55 000 
rev/min. Calculate the amplitude of the deflection of the mid-point of the shaft, and 
its phase with respect to the out-of-balance excitation when it is running at these speeds. 

In a new application it was proposed to run the turbo-charger at 45 000 rev/min. 
Comment on the implications of this proposal. 



The stiffness, k of the shaft at its mid-point is 48 £/// 3 , where E = 2 x 10 11 N/m 2 , 
I — 7 t( 0.01) 4 4/64 m 4 and l = 0.1 m. Hence 



k = 



48 x 2 x 10 11 
OT 3 



7r(0.01) 4 

64 



= 4.71 x 10 6 N/m. 



The effective mass, m eff , of the system can be taken to be the mass of the rotor plus 
one half of the mass of the shaft. Thus 
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m cf[ = 0.2 + j x 7850 x 0.1 



tt(O.OI) 2 

4 



= 0.231 kg. 



The natural frequency of the system is therefore 



4.71 x 10 6 
0.231 



= 4515 rad/s. 



which corresponds to 43 100 rev/min. 
The amplitude of shaft deflection is 



F 




where F = mrv 2 and v — 2n x 35 000/60 = 3663 rad/s. Thus 

F= 10" 5 x 3663 2 = 134.2 N, and - = — = 0.811. 

co 4515 



Hence the amplitude at 35 000 rev/min 



134.2 

4.71 x 10 6 (1 - 0.81 1 2 ) 



= 0.083 mm. 



and is in phase with the excitation because the excitation is below the resonance 
frequency, and the amplitude is therefore positive. 

At 55 000 rev/min, the amplitude is 



331 

4.71 x 10 6 (1 - 1.27 2 ) 



— 0.112 mm, 



which is out of phase with the excitation. 

Running the turbo-charger at 45 000 rev/min is potentially dangerous because of 
the proximity to resonance which occurs at 43 100 rev/min. 

It can be seen from Fig. 2.26, 2.28 and 2.29 that the system response at low 
frequencies (« co) is stiffness-dependant, and that in the region of resonance the 
response is damping-dependant, whereas at high frequencies ( » <x>) the response is 
governed by the system mass. It is most important to realize this when attempting 
to reduce the vibration of a system. For example the application of increased damping 
will have little effect if the excitation and response frequencies are in a region well 
away from resonance such as in that controlled by the system mass. 
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Before attempting to reduce the vibration levels in a machine or structure by 
increasing its damping, every effort should be made to reduce the vibration excitation 
at its source. It has to be accepted that many machines and processes generate a 
disturbing force of one sort or another, but the frequency of the disturbing force 
should not be at or near a natural frequency of the structure otherwise resonance 
will occur, with the resulting high amplitudes of vibration and dynamic stresses, and 
noise and fatigue problems. Resonance may also prevent the system fulfilling the 
desired function. 

Some reduction in excitation can often be achieved by changing the machinery 
generating the vibration, but this can usually only be done at the design stage. Resiting 
equipment may also effect some improvement. However, structural vibration caused 
by external excitation sources such as ground vibration, cross winds or turbulence 
from adjacent buildings can only be controlled by damping. 

In some machines vibrations are deliberately excited as part of the process, for 
example, in vibratory conveyors and compactors, and in ultrasonic welding. Naturally, 
nearby machines have to be protected from these vibrations. 

Rotating machinery such as fans, turbines, motors and propellers can generate 
disturbing forces at several different frequencies such as the rotating speed and blade 
passing frequency. Reciprocating machinery such as compressors and engines can 
rarely be perfectly balanced, and an exciting force is produced at the rotating speed 
and at harmonics. Strong vibration excitation in structures can also be caused by 
pressure fluctuation in gases and liquids flowing in pipes; and intermittent loads such 
as those imposed by lifts in buildings. 

There are two basic types of structural vibration; steady-state vibration caused by 
continually running machines such as engines, air-conditioning plants and generators 
either within the structure or situated in a neighbouring structure, and transient 
vibration caused by a short-duration disturbance such as a lorry or train passing 
over an expansion joint in a road or over a bridge. 

Some relief from steady-state vibration excitation can often be gained by moving 
the source of the excitation, since the mass of the vibration generator has some effect 
on the natural frequencies of the supporting structure. For example, in a building it 
may be an advantage to move mechanical equipment to a lower floor, and in a ship 
re-siting propulsion or service machinery may prove effective. The effect of local 
stiffening of the structure may prove to be disappointing, however, because by 
increasing the stiffness the mass is also increased, so that the change in the 
may prove to be very small. 

Occasionally a change in the vibration generating equipment can reduce vibration 
levels. For example a change in gear ratios is a mechanical drive system, or a change 
from a four-bladed to a three-bladed propeller in a ship propulsion system will alter 
the excitation frequency provided the speed of rotation is not changed. However, in 
many cases the running speeds of motors and engines are closely controlled as in 
electric generator sets, so there is no opportunity for changing the excitation Irequency. 

If vibration excitation cannot be reduced to acceptable levels, so that the system 
response is still too large, some measure of vibration isolation may be necessary; see 
section 2.3.2. 1 below. 
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2.3.2 Response of a viscous damped system supported on a foundation subjected to 
harmonic vibration 

The system considered is shown in Fig. 2.30. The foundation is subjected to harmonic 
vibration A sin vt and it is required to determine the response, x, of the body. 




= A sin vt 



Fig. 2.30. Single degree of freedom model of a vibrated system with viscous damping. 



The equation of motion is 

mx = c{y — x) + k(y — x). (2.18) 

If the displacement of the body relative to the foundation, u, is required, we may 
write u — x — y, and equation (2.18) becomes 

mu + cii + ku = — my — mv 2 A sin vt. 

This equation is similar to (2.14) so that the solution may be written directly as 

A(v/(o) 2 . 2{(v/oj) 

“ ^/((l — (v/co) 2 ) 2 + (2Cv/a>) 2 ) Sm ( tan 1 — (v/a)) 2 

If the absolute motion of the body is required we rewrite equation (2.18) as 
mx + cx + kx = cy + ky 

— cAv cos vt + kA sin vt 
= AyJ{k 2 + (cv) 2 ) sin (vt + a) 

where 



a = tan 




Hence, from the previous result, 



^y ; (fe 2 + (cv) 2 ) 

sj'((k - mv 2 ) 2 + (cv) 2 ) 



sin (vt — cf> + a). 
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The motion transmissibility is defined as the ratio of the amplitude of the absolute 
body vibration to the amplitude of the foundation vibration. Thus, 

X 

motion transmissibility = — 

A 




2.3. 2.1 Vibration isolation 

The dynamic forces produced by machinery are often very large. However, the force 
transmitted to the foundation or supporting structure can be reduced by using flexible 
mountings with the correct properties; alternatively a machine can be isolated from 
foundation vibration by using the correct flexible mountings. Fig. 2.31 shows a model 
of such a system. 

The force transmitted to the foundation is the sum of the spring force and the 




Fig. 2.31. Single degree of freedom system with foundation. 

damper force. Thus the transmitted force = kx + cx and F T , the amplitude of the 
transmitted force is given by 

Ft = y/UkX) 2 + ( cv *) 2 ]- 

The force transmission ratio or transmissibility, T R , is given by 

„ F T Xjlk 2 + (cv) 2 ] 

/r “ F ~ F 



since 
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Therefore the force and motion transmissibilities are the same. 

The effect of v/co on 7^ is shown in Fig. 2.32. It can be seen that for good isolation 
v/co > yjl, hence a low value of co is required which implies a low stiffness; that is, 
a flexible mounting: this may not always be acceptable in practice where a certain 
minimum stiffness is usually necessary to satisfy operating criteria. 




It is particularly important to be able to isolate vibration sources because structure 
borne vibration can otherwise be easily transmitted to parts which radiate well, and 
serious noise problems can occur. Theoretically, low stiffness isolators are desirable 
to give a low natural frequency. However, this often results in isolators which are 
too soft and stability problems may arise. The system can be attached rigidly to a 
large block which effectively increases its mass so that stiffer isolators can be used. 
The centre of mass of the combined system is also lowered, giving improved stability. 
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For the best response a mounting system may be designed with snubbers which 
control the large amplitudes while providing little or no damping when the amplitudes 
are small. 

There are four types of spring material commonly used for resilient mountings and 
vibration isolation: air, metal, rubber and cork. Air springs can be used for very 
low-frequency suspensions: resonance frequencies as low as 1 Hz can be achieved 
whereas metal springs can only be used for resonance frequencies greater than about 
1.3 Hz. Metal springs can transmit high frequencies, however, so rubber or felt pads 
are often used to prohibit metal-to-metal contact between the spring and the structure. 
Different forms of spring element can be used such as coil, torsion, cantilever and 
beam. Rubber can be used in shear or compression but rarely in tension. It is important 
to determine the dynamic stiffness of a rubber isolator because this is generally much 
greater than the static stiffness. Also rubber possesses some inherent damping although 
this may be sensitive to amplitude, frequency and temperature. Natural frequencies 
from 5 Hz upwards can be achieved. Cork is one of the oldest materials used for 
vibration isolation. It is usually used in compression and natural frequencies of 25 
Hz upwards are typical. 

It is shown in section 2.3.2 that good vibration isolation, that is low force and 
motion transmissibility, can be achieved by supporting the vibration generator on a 
flexible low-frequency mounting. Thus although disturbing forces are generated, only 
a small proportion of them are transmitted to the supporting structure. However, 
this theory assumes that a mode of vibration is excited by a harmonic force passing 
through the centre of mass of the installation; although this is often a reasonable 
approximation it rarely actually occurs in practice because, due to a lack of symmetry 
of the supported machine, several different mountings may be needed to achieve a 
level installation, and the mass centre is seldom in the same plane as the tops of the 
mountings. 

Thus the mounting which provides good isolation against a vertical exciting force 
may allow excessive horizontal motion, because of a frequency component close to 
the natural frequency of the horizontal mode of vibration. Also a secondary exciting 
force acting eccentrically from the centre of mass can excite large rotation amplitudes 
when the frequency is near to that of a rocking mode of an installation. 

To limit the motion of a machine installation that generates harmonic forces and 
moments the mass and inertia of the installation supported by the mountings may 
have to be increased; that is an inertia block may have to be added to the installation. 
If non-metallic mountings are used the dynamic stiffness at the frequencies of interest 
will have to be found, probably by carrying out further dynamic tests in which the 
mounting is correctly loaded; they may also possess curious damping characteristics 
which may be included in the analysis by using the concept of complex stiffness, as 
discussed in section 2.2.4. 

Air bags or bellows are sometimes used for very low-frequency mountings where 
some swaying of the supported system is acceptable. This is an important consideration 
because if the motion of the inertia block and the machinery is large, pipework and 
other services may be overstressed, which can lead to fatigue failure of these 
components. Approximate analysis shows that the natural frequency of a body 
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supported on bellows filled with air under pressure is inversely proportional to the 
square root of the volume of the bellows, so that a change in natural frequency can 
be effected simply by a change in bellows volume. This can easily be achieved by 
opening or closing valves connecting the bellows to additional receivers, or by adding 
a liquid to the bellows. Natural frequencies of 0.5 Hz, or even less, are obtainable . 
An additional advantage of air suspension is that the system can be made self-levelling, 
when fitted with suitable valves and an air supply. Air pressures of about 5 to 10 
times atmospheric pressure are usual. 

Greater attenuation of the exciting force at high frequencies can be achieved by 
using a two-stage mounting. In this arrangement the machine is set on flexible 
mountings on an inertia block, which is itself supported by flexible mountings. This 
may not be too expensive to install since in many cases an existing sub-frame or 
structure can be used as the inertia block. If a floating floor in a building is used as 
the inertia block, some allowance must be made for the additional stiffness arising 
from the air space below it. This can be found by measuring the dynamic stiffness of 
the floor by means of resonance tests. 

Naturally, techniques used for isolating structures from exciting forces arising in 
machinery and plant can also be used for isolating delicate equipment from vibrations 
in the structure. For example, sensitive electrical equipment in ships can be isolated 
from hull vibration, and operating tables and metrology equipment can be isolated 
from building vibration. 

The above isolation systems are all passive; an active isolation system is one in 
which the exciting force or moment is applied by an externally powered force or 
couple. The opposing force or moment can be produced by means such as hydraulic 
rams, out-of-balance rotating bodies or electro-magnetism. Naturally it is essential 
to have accurate phase and amplitude control, to ensure that the opposing force is 
always equal, and opposite, to the exciting force. Although active isolation systems 
can be expensive to install, excellent results are obtainable so that the supporting 
structure is kept almost completely still. However it must be noted that force actuators 
such as hydraulic rams must react on another part of the system 

If, after careful selection and design of machinery and equipment, careful installation 
and commissioning, and carrying out isolation as necessary the vibration levels in 
the system are still too large, then some increase in the damping is necessary. This 
is also the case when excitation occurs from sources beyond the designers’ control 
such as cross winds, earthquakes and currents. 



Example 13 

A spring-mounted body moves with a velocity v along an undulating surface, as shown. 

The body has a mass m and is connected to the wheel by a spring of stiffness k, 
and a viscous damper whose damping coefficient is c. The undulating surface has a 
wavelength Land an amplitude h. 

Derive an expression for the ratio of amplitudes of the absolute vertical displacement 
of the body to the surface undulations. 
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Hence the equation of motion is 

mx — — k(x — y) — c(x — y), 
or 

mx + cx + kx — cy + ky. 

Now 

y = h cos vt, so 

mx + cx + kx — J[k 2 + (cv) 2 ]h sin (vt + 4 >)• 
Hence, if 

x = X 0 sin (vt 4- a), 




Example 14 

The vibration on the floor in a building is SHM at a frequency in the range 15-60 Hz. 
It is desired to install sensitive equipment in the building which must be insulated 
from floor vibration. The equipment is fastened to a small platform which is supported 
by three similar springs resting on the floor, each carrying an equal load. Only vertical 
motion occurs. The combined mass of the equipment and platform is 40 kg, and the 
equivalent viscous damping ratio of the suspension is 0.2. 

Find the maximum value for the spring stiffness, if the amplitude of transmitted 
vibration is to be less than 10% of the floor vibration over the given frequency range. 
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that is 




17.84 




2 

-99 = 0. 



Hence 



co 

When 

v = 15.271 rad/s, co = 19.97 rad/s. 

Since 



co = 




and m = 40 kg, 



total k = 15 935 N/m, 

that is, stiffness of each spring = 15 935/3 N/m = 5.3 kN/m. 

The amplitude of the transmitted vibration will be less than 10% at frequencies 
above 15 Hz. 



Example 15 

A machine of mass m generates a disturbing force F sin vt; to reduce the force 
transmitted to the supporting structure, the machine is mounted on a spring of stiffness 
k with a damper in parallel. Compare the effectiveness of this isolation system for 
viscous and hysteretic damping. 




Viscous damping. From section 2.3.2. 1, 




Hysteretic damping. From section 2.2.5, 

CV V 

Putting rj = — = 2(— , 
k co 
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The effectiveness of these isolators can be compared using these expressions for T R . 
The results are given in the table below. 

It can be seen that the isolation effects are similar for the viscous and hysteretically 
damped isolators, except at high frequency ratios when the hysteretic damping gives 
much better attentuation of 7^. At these frequencies it is better to decouple the viscous 
damped isolator by attaching small springs or rubber bushes at each end. 





Viscously damped 
isolator 


Hysteretically damped 
isolator 


Value of T r when v = 0 + 


1 


1 


Frequency ratio v/co for 






resonance 


1 


1 


Value of T R at resonance 


Vn + (20 2 ] i 


Vd + '7 2 )_ 1 




2£ 2£ 


n v 


Value of T r when 






v/co = J2 


1 


l 


Frequency ratio v/co for 






isolation. 


>V2 


>y/2 


High frequency, v/co 1, 


2C 


1 


attenuation of 7^ 


v/co 


(v/co) 2 



Example 16 

A single degree of freedom vibrational system of very small viscous damping (£ < 0.1) 
is excited by an harmonic force of frequency v and amplitude F. Show that the Q 
factor of the system is equal to the reciprocal of twice the damping ratio (. The Q 
factor is equal to (X/X s ) max . 

It is sometimes difficult to measure £ in this way because the static deflection X s 
of the body under a force F is very small. Another way is to obtain the two frequencies 
Pi and p 2 (one either side of the resonance frequency co) at the half- power points. 

Show that Q — l/2£ = co/(p 2 — Pi) = co/Aco. 

(The half-power points are those points on the response curve with an amplitude 
1/^/2 times the amplitude at resonance). 

From equation (2.16) above 

x = {F/kyJd 1 - (v/co) 2 ) 2 + (2£v/co) 2 ) 

If v = 0, X s = F/k, and at resonance v — a>, 
so 



*ma, = (F/k)/ 2£ = X S /2C, 
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that is 



Q = = m- 

If X s cannot be determined, the Q factor can be found by using the half-power point 
method. This method requires very accurate measurement of the vibration amplitude 
for excitation frequencies in the region of resonance. Once X max and co have been 
located, the so called half-power points are found when the amplitude is X p = X ma Jj2 
and the corresponding frequencies either side of 10 , and p 2 determined. Since the 
energy dissipated per cycle is proportional to X 2 , the energy dissipated is reduced 
by 50% when the amplitude is reduced by a factor of l/y/2. 




Now 



X = 



F/k 



[‘-G)T + fcJ} 



Thus 



F/k ( v 

*max = C small, so X max occurs at - ~ I }, 

2C \ <° 
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and 



y _* m »_ f/fc 

p V 2 V 2 -2C 



F/Jk 



1 - 



+ \ltr 

co / J L CO 



Hence 



1 - 



CO 



+ [«i J = «> 



and 



J = (1 - 2C 2 ) ± 2CV(1 - C 2 ). 

That is 

2 _ 2 

= 4C/(1 - c 2 ) - 4£, if C is small. 

CD 

Now 

P2 - Pi = / Pa ~ Pi V Pa + PA = J Pz ~Pi \ 

CO 2 \ CO )\ 00 J \ c of 

because (p x + p 2 )/oo = 2. That is, a symmetrical response curve is assumed for small 
c. Thus 

Pa - Pi = 2 ; = — = 1 
to ' co Q' 

where Aco is the frequency bandwidth at the half-power points. 

Therefore, for light damping, the damping ratio (, and hence the Q factor, associated 
with any mode of vibration can be found from the amplitude-frequency measurements 
at resonance and the half-power points. Care is needed to ensure that the exciting 
device does not load the system and alter the frequency response and the damping. 
It should be noted that some difficulty is often encountered in measuring X max 
accurately. 

In real systems and structures a very high Q at a low frequency, or a very low Q 
at a high frequency, seldom occur, but it can be appreciated from the above that very 
real measuring difficulties can be encountered when trying to measure bandwidths 
of only a few Hz accurately, even if the amplitude of vibration can be determined. 
The table below shows the relationship between Q and A/ for different values of 
frequency. 
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Q factor 


500 


50 


5 


Resonance frequency (Hz) 
10 


0.02 


Frequency bandwidth (Hz) 
0.2 


2 


100 


0.2 


2 


20 


1000 


2 


20 


200 



An improvement in accuracy in determining Q can often be obtained by measuring 
both amplitude and phase of the response for a range of exciting frequencies. Consider 
a single degree of freedom system under forced excitation Fe JVI . The equation of 
motion is 

m'x + cx + kx = Fe jv '. 

A solution x = Xe^' can be assumed, so that 
— mv 2 X + j cvX + kX = F. 

Hence 



X 1 

F (k — mv 2 ) + jcv 

k — mv 2 cv 

(k — mv 2 ) 2 + (cv) 2 \k — mv 2 ) 2 4- (cv) 2 

That is, X/F consists of two vectors, R e(X/F) in phase with the force, and Im(A7F) 
in quadrature with the force. The locus of the end point of vector X/F as v varies is 
shown opposite for a given value of c. This is obtained by calculating real and 
imaginary components of X/F for a range of frequencies. 

Experimentally this curve can be obtained by plotting the measured amplitude and 
phase of (X/F) for each exciting frequency. 

Since 



k — mv 2 

tan (p = , 

cv 



when (f) = 45° and 135°, 
k — mp\ 



1 = 



cp i 



and 



Hence 



k — mp\ 
cp 2 



mp\ + cp i — k — 0 and mp\ — cp 2 — k = 0. 
Subtracting one equation from the other gives p 2 ~ Pi = c/m, or 
Pi -Pi 1 



to 



to 








Sec. 2.31 



Forced vibration 67 




That is, X/F at resonance lies along the imaginary axis, and the half-power points 
occur when 4> = 45° and 135°. If experimental results are plotted on these axes a 
smooth curve can be drawn through them so that the half-power points can be 
accurately located. 

The method is also effective when the damping is hysteretic, because in this case 

X _ 1 

F (k — mv 2 ) + j r]k' 

so that 





68 The vibrations of systems having one degree of freedom 



[Ch. 2 



That is, the locus of (X/F) as v increases from zero is part of a circle, centre 
{0, — l/2tik) and radius 1/2 rjk, as shown. 




OJ 



In this case therefore it is particularly easy to draw an accurate locus from a few 
experimental results, and p j and p 2 are located on the horizontal diameter of the circle. 

This technique is known variously as a frequency locus plot, Kennedy-Pancu 
diagram, or Nyquist diagram. 

It must be realized that the assessment of damping can only be approximate. It is 
difficult to obtain accurate, reliable, experimental data, particularly in the region of 
resonance; the analysis will depend on whether viscous or hysteretic damping is 
assumed, and some non-linearity may occur in a real system. These effects may cause 
the frequency-locus plot to rotate and translate in the Re(X/F), Im(2f/F) plane. In 
these cases the resonance frequency can be found from that part of the plot where 
the greatest rate of change of phase with frequency occurs. 

Example 17 

A machine of mass 550 kg is flexibly supported on rubber mountings which provide 
a force proportional to displacement of 210 kN/m, together with a viscous damping 
force. The machine gives an exciting force of the form Rv 2 cos vt, where R is a 
constant. At very high speeds of rotation, the measured amplitude of vibration is 
0.25 mm, and the maximum amplitude recorded as the speed is slowly increased from 
zero is 2 mm. Find the value of R and the damping ratio. 
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Now, 

X = Rv 2 /(yJ(k — mv 2 ) 2 + c 2 v 2 ) 

If v is large, 

X -*■ Rv 2 /(7(m 2 v 4 )) = R/m. 

Hence 

R = ml = (550 x 0.25)/1000 = 0.1375 kg m. 

For maximum X, dX/dv = 0, hence v 2 = 2k 2 /(2mk — c 2 ), and 
R/2m 

mai cva - c 2 )' 

So 

CV(1 - C 2 ) = R/(2mX m J = 0.1375/(2 x 550 x 2 x 10" 3 ) 

- 0.062, 

that is, 

c = 0.062. 

23.3 Response of a Coulomb damped system to a simple harmonic exciting force 
with constant amplitude 




Fig. 2.33. Single degree of freedom model of a forced system with Coulomb damping. 

In the system shown in Fig. 2.33 the damper relies upon dry friction. 

The equation of motion is non-linear because the constant friction force F d always 
opposes the motion: 

mx + kx ± F d = F sin vt. 

If F d is large compared to F, discontinuous motion will occur, but in most systems 
F d is usually small so that an approximate continuous solution is valid. The 
approximate solution is obtained by linearizing the equation of motion; this can be 
done by expressing F d in terms of an equivalent viscous damping coefficient, c d . From 
section 2.2.5, 
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c d 



4 



The solution to the linearized equation of motion gives the amplitude X of the 
motion as 



F 

X ~ yjl(k ~ ™ 2 ) 2 + (c d v) 2 ]' 

Hence 



F 

X -JiT-mv 2 ) 2 +(4 FJkX) 2 J 
That is, 

X V(l-(4 FJnF) 2 ) 

X s 1 - (v/o;) 2 

This expression is satisfactory for small damping forces, but breaks down if 
AFJnF < 1, that is F d > (jt/4)F. 

At resonance the amplitude is not limited by Coulomb friction. 

23.4 Response of a hysteretically damped system to a simple harmonic exciting force 
with constant amplitude 

In the single degree of freedom model shown in Fig. 2.34 the damping is hysteretic. 




Fig. 2.34. Single degree of freedom model of a forced system with hysteric damping. 



The equation of motion is 
mx + k*x = F sin vt. 

Since 

k* = k( 1 + j rj), 

F sin vf 
(k — mv 2 ) + j tjk' 



x = 
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and 



X 1 

Y s ~ jdi-iv/coyy + ^r 

This result can also be obtained from the analysis of a viscous damped system by 
substituting c = rjk/v. 

It should be noted that if c = t]k/v, at resonance c = r]yj(km), that is rj = 2£ = 1 /Q. 

2.3.5 Response of a system to a suddenly applied force 



rWWV- 





+ve 






m 




9 r 



Fig. 2.35. Single degree of freedom model with constant exciting force. 



Consider a single degree of freedom undamped system, such as the system shown in 
Fig. 2.35, which has been subjected to a suddenly applied force F. The equation of 
motion is m'x + kx = F. The solution to this equation comprises a complementary 
function A sin cot + B cos cot, where c o = ,/(k/m) rad/s together with a particular 
solution. The particular solution may be found by using the D-operator. Thus the 
equation of motion can be written 



and 




F 

k’ 



x = 






F 

k' 



That is, the complete solution to the equation of motion is 



x — A sin cot + B cos cot + — . 

k 



If the initial conditions are such that x = x = 0 at t = 0, then B = — F/k and A = 0. 
Hence 
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The motion is shown in Fig. 2.36. It will be seen that the maximum dynamic 
displacement is twice the static displacement occurring under the same load. This is 
an important consideration in systems subjected to suddenly applied loads. 

If the system possesses viscous damping of coefficient c, the solution to the equation 
of motion is x = Xe _£ "' sin (co v t + </>) + F/k. 

With the same initial conditions as above, 

* = sin (“s/o - « ' + ,a,r ' 

This reduces to the undamped case if £ = 0. The response of the damped system is 
shown in Fig. 2.37. 




Fig. 2.36. Displacement-time response for system shown in Fig. 2.35. 




Fig. 2.37. Displacement-time response for single degree of freedom system with viscous damping. 



2.3.6 Shock excitation 

Some systems are subjected to shock or impulse loads arising from suddenly applied, 
non-periodic, short duration exciting forces. 

The impulsive force shown in Fig. 2.38 consists of a force of magnitude F max /e 
which has a time duration of e. 

The impulse is equal to 
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Fig. 2.38. Impulse. 



When F max is equal to unity, the force in the limiting case e-»0 is called either the 
unit impulse or the delta function, and is identified by the symbol S(t — £,), where 

r x 

d(t - f)d£ = 1. 

Jo 

Since F df = m dv, the impulse F max acting on a body will result in a sudden change 
in its velocity without an appreciable change in its displacement. Thus the motion 
of a single degree of freedom system excited by an impulse F max corresponds to free 
vibration with initial conditions x = 0 and x = v 0 = F ma Jm at t = 0. 

Once the response g(t ) say, to a unit impulse excitation is known, it is possible to 
establish the equation for the response of a system to an arbitrary exciting force F{t). 
For this the arbitrary pulse in considered to comprise a series of impulses as shown 
in Fig. 2.39. 




Fig. 2.39. Force-time pulse. 
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If one of the impulses is examined which starts at time £, its magnitude is F(£) A£, 
and its contribution to the system response at time f is found by replacing the time 
with the elapsed time (t — £) as shown in Fig. 2.40. 

If the system can be assumed to be linear, the principle of superposition can be 
applied, so that 



x(f) = 



F(Og(t - f)d£. 



This is known as the Duhamel integral. 




Fig. 2.40. Displacement-time response to impulse. 



2.3.7 Harmonic analysis 

A function which is periodic but not harmonic can be represented by the sum of a 
number of terms, each term representing some multiple of the fundamental frequency. 
In a linear system each of these harmonic terms acts as if it alone were exciting the 
system, and the system response is the sum of the excitation of all the harmonics. 

For example, if the periodic forcing function of a single degree of freedom undamped 
system is 

F 2 sin (vt + <Xj) + F 2 sin (2vt + a 2 ) + F 3 sin (3vr + a 3 ) + • • • + F n sin (nvt + a„), 
the steady-state response to F 1 sin (vr + aj is 







sin (vt + aj. 



and the response to F 2 sin (2vt + a 2 ) is 
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x 2 




sin (2vf 4- a 2 ). 



and so on, so that 



x 



n 



= z 




sin ( nvt + a„). 



Clearly that harmonic which is closest to the system natural frequency will most 
influence the response. 

A periodic function can be written as the sum of a number of harmonic terms by 
writing a Fourier series for the function. A Fourier series can be written 



F(t) = -y + Yj ( fl n COS nVt + S i n nV 0> 
^ n = 1 



where 



2 f r 

“"-Jo 



2 

= ~ 
x 



F(t)dt, 

F{t ) cos vtdt, and 
F(t) sin vfdf. 



For example, consider the first four terms of the Fourier series representation of the 
square wave shown in Fig. 2.41 to be required; r = 2?r so v = 1 rad/s. 




Fig. 2.41. Square wave. 
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Uq 

F(t) = — + a 1 cos vt + a 2 cos 2vr + 



+ by sin vt + b 2 sin 2vt + • • • 



-SI 

■hi 



-I Fm= Tn 



1 dr + - — 
In 



1 dr = 0. 



F(t) cos vtdt 



2 _ 

2n 



Similarly 

a 7 = a , 



cos vtdt H 

2n 



0. 



— cos vrdr = 0. 



by = - F(t) sin vtdt 



— sin vrdr + 



_ 2 r* 

Jo 



2n 



— sin vrdr 
Jn 

1 1,4 

= — [- cos vr]o -\ — [cos vr]„ - — . 
7iv nv nv 



Since v = 1 rad/s, 



It is found that b 2 — 0, b, = 4/3 n and so on. 



Thus F(r) = 



1 . , 1 . . \ . n 
sin r + - sin 3r + - sin 5r + - sin 7r + • • 






the series representation 



of the square wave shown. 

If this stimulus is applied to a simple undamped system with to = 4 rad/s, say, the 
steady-state response is given by 



4 . 4 . 4 4 

- sin r — sin 3r — sin 5r — sin It 
n 3 n 5 n In 

x= y^w + 1-® 2 + 1 - (it 2 + 1-© 2 ' 

that is, x = 1.36 sin r + 0.97 sin 3t — 0.45 sin 5t — 0.09 sin It — ■ ■ - 
Usually three or four terms of the series dominate the predicted response. 

It is worth sketching the components of F(t) above to show that they produce a 
reasonable square wave, whereas the components of x do not. This is an important 
result. 
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2.3.8 Random vibration 

If the vibration response parameters of a dynamic system are accurately known as 
functions of time, the vibration is said to be deterministic. However, in many systems 
and processes responses cannot be accurately predicted; these are called random 
processes. Examples of a random process are turbulence, fatigue, the meshing of 
imperfect gears, surface irregularities, the motion of a car running along a rough 
road, and building vibration excited by an earthquake. Fig. 2.42 shows a random 
process. 




A collection of sample functions Xjfr), x 2 (f), x 3 (r), ..., x„(r) which make up the 
random process x(f) is called an ensemble, as shown in Fig. 2.43. These functions may 
comprise, for example, records of noise, pressure fluctuations or vibration levels, taken 
under the same conditions but at different times. 

Any quantity which cannot be precisely predicted is non-deterministic and is known 
as a random variable or a probabilistic quantity. That is, if a series of tests is conducted 
to find the value of a particular parameter x and that value is found to vary in an 
unpredictable way that is not a function of any other parameter, than x is a random 
variable. 

2. 3. 8.1 Probability distribution 

If n experimental values of a variable x are x l5 x 2 , x 3 , . . . , x„, the probability that the 
value of x will be less than x' is n'/n, where n is the number of x values which are 
less than or equal to x'. That is, 

Prob (x ^ x ) = n'/n. 

When n approaches x this expression is the probability distribution function of x, 
denoted by P(x), so that 
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x„<r) 




Fig. 2.43. Ensemble of a random process. 



P(x ) = Lt (n'/n) 

n~*o o 

The typical variation of P(x) with x is shown in Fig. 2.44. Since x(t) denotes a physical 
quantity, 

Prob (x< — oo) = 0, and Prob (x< + co) = 1. 

The probability density function is the derivative of P(x) with respect to x and this 
is denoted by p(x). That is, 



P(x) 



dPjx) 

d(x) 



= Lt 

Ajc-» 0 



P(x + Ax) - P(x f 
Ax 



where P(x + Ax) — P(x) is the probability that the value of x{t) will lie between x 
and x + Ax (Fig. 2.44). Now, 



p(x) = 



dP(x) 

dx 
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E\_x(t j)] = E[x(ty + t)] 
for any value of t. 

If ^x) = .x, the expected value or mean value of x, 



£[x] or x, is 



E[x] 



-c 



xp(x)dx. 



In addition, if^x) = x 2 , the mean square value x 2 of x is 



£[x 2 ] 






x 2 p(x) dx. 



The variance of x, a 2 is the mean square value of x about the mean, that is, 



a 2 = E[(x - x) 2 ] - 



(x — x) 2 p(x)dx = (x 2 ) — (x) 2 . 



a is the standard deviation of x, hence 

variance = (standard deviation) 2 = (mean square -(mean) 2 } 

If two or more random variables x x and x 2 , represent a random process at two 
different instants of time, then 



E\J(x i,x 2 )] = 



/(x 1 ,x 2 )p(x 1 ,x 2 )dx 1 dx 2 , 



and if t x and t 2 are the two instants of time, 

E[x{t 1 ),x(t 2 )'] = ,t 2 \ 

which is the auto-correlation function for the random process (Fig 2.45). 

For random processes which are stationary, 

£[x(fj, x(r 2 )] = R(t^ t 2 ) = R(t 2 - tj) = R(r), 

say, since the average depends only on time differences. If the process is also ergodic, 
then 



= l« i r 

r-ao J -i 



R( T) 

It is worth noting that 

R( 0) = £[x(t) 2 ] 



x(f)x(f + r)d t. 



T-*oc J -1 



x 2 (t)dr, 



which is the average power in a sample function. 



2. 3. 8. 2 Random processes 

The most important random process is the Gaussian, or normal random process. This 
is because a wide range of physically observed random waveforms can be represented 
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as Gaussian processes, and the process has mathematical features which make analysis 
relatively straightforward. 

The probability density function of a Gaussian process x(f) is 



p(x) = 



V(^ exp 




where a is the standard deviation of x, and x is the mean value of x. 

The values of o and x may vary with time for a non-stationary process but are 
independant of time if the process is stationary. 

One of the most important features of the Gaussian process is that the response 
of a linear system to this form of excitation is usually another, but still Guassian, 
random process. The only changes are that the magnitude and standard deviation 
of the response may differ from those of the excitation. 

A Gaussian probability density function is shown in Fig. 2.46. It can be seen to 
be symmetric about the mean value x, and the standard deviation a controls the spread. 




— Xo 0 +Xa 



Fig. 2.47. Gaussian probability density function with zero mean. 
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The probability that x(r) lies between — act and + act, where A is a positive number, 
can be found since, if x = 0, 



Prob { — act ^ x(r) ^ + act} 



7<^ exp 




Fig. 2.47 shows the Gaussian probability density function with zero mean. 

This integral has been calculated for a range of values of a; the results are given 
below. The probability that x(f) lies outside the range — act to + act is 1 minus the 
value of the above integral; this probability is also given. 



Value of a 


Prob [ — act ^ x(f) < act] 


Prob [|x(t)| > act] 


0 


0 


1.0000 


0.2 


0.1585 


0.8415 


0.4 


0.3108 


0.6892 


0.6 


0.4515 


0.5485 


0.8 


0.5763 


0.4237 


1.0 


0.6827 


0.3173 


1.2 


0.7699 


0.2301 


1.4 


0.8586 


0.1414 


1.6 


0.8904 


0.1096 


1.8 


0.9281 


0.0719 


2.0 


0.9545 


0.0455 


2.2 


0.9722 


0.0278 


2.4 


0.9835 


0.0164 


2.6 


0.9907 


0.0093 


2.8 


0.9949 


0.0051 


3.0 


0.9973 


0.0027 


3.2 


0.9986 


0.00137 


3.4 


0.9993 


0.00067 


3.6 


0.9997 


0.00032 


3.8 


0.9998 


0.00014 


4.0 


0.9999 


0.00006 
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2.3.8.3 Spectral density 

The spectral density S(a> ) of a stationary random process is the Fourier transform of 
the autocorrelation function R(t). It is given by 

K(t)e- j<Bt dr. 

The inverse, which also holds true, is 



m = 



* 00 

S(m)e“ j “ r d(w 

J — oc 



Ift - 0 



R(0) = 



I 



00 

S(o;)doj = E[x 2 ] 

- GO 



That is, the mean square value of a stationary random process x is the area under 
the S(co) against frequency curve. A typical spectral density function is shown in Fig. 



2.48. 



S(«) 




Fig. 2.48. Typical spectral density function. 



A random process whose spectral density is constant over a very wide frequency 
range is called white noise. If the spectral density of a process has a significant value 
over a narrower range of frequencies, but one which is nevertheless still wide compared 
with the centre frequency of the band, it is termed a wide-band process (Fig. 2.49). If 
the frequency range is narrow compared with the centre frequency it is termed a 
narrow-band process (Fig. 2.50). Narrow-band processes frequently occur in engin- 
eering practice because real systems often respond strongly to specific exciting 
frequencies and thereby effectively act as a filter. 
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2.3.9 The measurement of vibration 

The most commonly used device for vibration measurement is the piezoelectric 
accelerometer, which gives an electric signal proportional to the vibration acceleration. 
This signal can readily be amplified, analysed, displayed, recorded, and so on. The 
principles of this device can be studied by referring to Fig. 2.51 which shows a body 
of mass m supported by an elastic system of stiffness k and effective viscous damping 
of coefficient c. 

This dynamic system is usually enclosed in a case which is fastened to the surface 
whose vibration is to be measured. The body has a pointer fixed to it, which moves 
over a scale fastened to the case: that is it measures u, the motion of the suspended 
body relative to that of the vibrating surface so that u = x — y. 

Now from section 2.3.2, the amplitude of u is 




so that if a) is low and v $> a>. 




Fig. 2.51. Vibration measuring device. 
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U ~ - A 

~ (v/o>) 2 A 

that is, the device measures the input vibration amplitude; when operating in this 
mode it is called a vibrometer , and if u> is high so that <x>$>v, then 

AMa>) ‘ = -U 2 . 



U 



1 



co 



that is, the device measures the input vibration acceleration amplitude; when operating 
in this mode it is called an accelerometer. 

By adjusting the system parameters correctly it is possible to make 




have a value close to unity for exciting frequencies v up to about 0.3to. Commercial 
accelerometers usually have piezoelectric elements instead of a spring and damper, 
so that the electric signal produced is proportional to the relative motion u, above. 




3 



The vibrations of systems having more than one 
degree of freedom 



Many real systems can be represented by a single degree of freedom model. However, 
most actual systems have several bodies and several restraints and therefore several 
degrees of freedom. The number of degrees of freedom that a system possesses is 
equal to the number of independent coordinates necessary to describe the motion of 
the system. Since no body is completely rigid, and no spring is without mass, every 
real system has more than one degree of freedom, and sometimes it is not sufficiently 
realistic to approximate a system by a single degree of freedom model. Thus, it is 
necessary to study the vibration of systems with more than one degree of freedom. 

Each flexibly connected body in a multi-degree of freedom system can move 
independently of the other bodies, and only under certain conditions will all bodies 
undergo an harmonic motion at the same frequency. Since all bodies move with the 
same frequency, they all attain their amplitudes at the same time even if they do not 
all move in the same direction. When such motion occurs the frequency is called a 
natural frequency of the system, and the motion is a principal mode of vibration: the 
number of natural frequencies and principal modes that a system possesses is equal 
to the number of degrees of freedom of that system. The deployment of the system 
at its lowest or first natural frequency is called its first mode, at the next highest or 
second natural frequency it is called the second mode and so on. 

A two degree of freedom system will be considered initially. This is because the 
addition of more degrees of freedom increases the labour of the solution procedure 
but does not introduce any new analytical principles. 

Initially, we will obtain the equations of motion for a two degree of freedom model, 
and from these find the natural frequencies and corresponding mode shapes. 

Some examples of two degree of freedom models of vibrating systems are shown 
in Figs 3.1(a)-(h)- 
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Fig. 3.1(b) Linear undamped system, vertical motion. Coordinates >’ x and y 2 . 




1 1 h 

Fig. 3.1(c) Torsional undamped system. Coordinates 6, and 0 2 . 
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Fig. 3.1(d) Coupled pendula. Coordinates 0, and <f> 2 . 





Figt. 3.1(e) System with combined translation and rotation. Coordinates x and 6 




Fig. 3.1(f) Shear frame. Coordinates x, and x 2 . 
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3.1 THE VIBRATION OF SYSTEMS WITH TWO DEGREES OF FREEDOM 

3.1.1 Free vibration of an undamped system 

Of the examples of two degree of freedom models shown in Figs. 3.1(a) — (h), 
consider the system shown in Fig. 3.1(a). If Xj > x 2 the FBDs are as shown in 
Figure 3.2. 



Vi 











1 




1 





"Vi 



k(x. — x„ 

t < 



k 2 x 2 



"V 2 



(a) (b) 

Fig. 3.2. (a) Applied forces; (b) effective forces. 



The equations of motion are therefore, 

m 1 x 1 = — k 1 x l — k(x 1 — x 2 ) for body 1, (3.1) 

and 

m 2 x 2 — k(x t — x 2 ) — k 2 x 2 for body 2. (3.2) 

The same equations are obtained if x x < x 2 is assumed because the direction of 
the central spring force is then reversed. 

Equations (3.1) and (3.2) can be solved for the natural frequencies and corresponding 
mode shapes by assuming a solution of the form 

Xi = j4jsin(a>f + ip) and x 2 = A 2 sin(cut + ip). 

This assumes that Xj and x 2 oscillate with the same frequency to and are either in 
phase or n out of phase. This is a sufficient condition to make to a natural frequency. 
Substituting these solutions into the equations of motion gives 

— Wj/lj^sin^f + ip) = — kjAjSi^cuf + i p) — k(A { — ^ 2 )sin(a)t + i p) 

and 

— m 2 /4 2 a) 2 sin(cut + t p) = k(A t — A 2 )sin(a>t + i p) — k 2 A 2 s\n(a>t + i p). 

Since these solutions are true for all values of t, 

A x (k + — mjeo 2 ) + A 2 ( — k) = 0 (3.3) 

and 



A t ( — k) + A 2 (k 2 + k — m 2 to 2 ) = 0. 



(3.4) 
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A x and A 2 can be eliminated by writing 



k + k 1 — WiO) 2 —k 

— k k + k 2 — m 2 o) 2 



= 0 



(3.5) 



This is the characteristic or frequency equation. Alternatively, we may write 
AJA 2 = k/(k + k x — mjco 2 ) from (3.3) 'l 

and 



Thus 



and 



AJA 2 = (k 2 + k — m 2 a> 2 )/k from (3.4) 



k/(k + k 1 — m t a> 2 ) = (k 2 4- k — m 2 a> 2 )/k 



(3.6) 



(fc + /cj — mjOj 2 Xfc 2 + k — m 2 co 2 ) — k 2 = 0. 



(3.7) 



This result is the frequency equation and could also be obtained by multiplying out 
the above determinant, equation (3.5). 

The solutions to equation (3.7) give the natural frequencies of free vibration for 
the system in Fig. 3.1(a). The corresponding mode shapes are found by substituting 
these frequencies, in turn, into either of equations (3.6). 

Consider the case when k 1 = k 2 = k, and m t = m 2 = m. The frequency equation is 
(2k — mco 2 ) 2 — k 2 = 0; that is, 

m 2 co 4 — 4 mkco 2 + 3k 2 = 0, or (mco 2 — k)(ma) 2 — 3k) - 0. 

Therefore, either mco 2 — k = 0, or mco 2 — 3k = 0. 

Thus 

co ! = yj(k/m) rad/s and co 2 = ^/pk/m) rad/s. 

If 

co = yj(k/m) rad/s, (A l /A 2 ) a >- j ik /m) = + 1 , 

and if 

(o = yJ(3k/m) rad/s, (A 1 /A 2 ) 0J = j Wm) = - 1. (from 3.6) 

This gives the mode shapes corresponding to the frequencies co t and co 2 . Thus, the 
first mode of free vibration occurs at a frequency (l/2n)yJ(k/m)Hz and (A j /Aj = 1, 
that is, the bodies move in phase with each other and with the same amplitude as if 
connected by a rigid link, Fig. 3.3. The second mode of free vibration occurs at a 
frequency (l/27t) v /(3k/m)Hz and {A x j A 2 ) 11 = — 1, that is, the bodies move exactly out 
of phase with each other, but with the same amplitude, see Fig. 3.3. 
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Fig. 3.3. Natural frequencies and mode shapes for two degree of freedom translation vibration system. 
Bodies of equal mass and springs of equal stiffness. 



3.1.2 Free motion 

The two modes of vibration can be written 



1 

A 



sin(cOjt + y>i), 



and 





'"-1 


\M 


l* 2 J 


1 1 


[Ail 



sin(<y 2 f + ij/ 2 ). 



where the ratio A 1 /A 2 is specified for each mode. 

Since each solution satisfies the equation of motion, the general solution is 



f*il 


L-J 


\A,\ 


l* 2 J 


M 


t/f 2 j 



sin^f + i/tj) + 






sin(co 2 r + \f/ 2 ). 



where A v A 2 , 1 A 1 , 4>i are found from the initial conditions. 
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For example, for the system considered above, if one body is displaced a distance 
X and released, 

x^O) = X and x 2 (0) = Xj(0) = x 2 (0) = 0 

where x t (0) means the value of Xl when t = 0, and similarly for x 2 (0), X[(0) and x 2 (0). 
Remembering that in this system cu, = yj(k/m), to 2 = yj(3k/m), and 




we can write 

Xj = sin Q(k/m)t + i p x ) + s\n(yj{3k/m)t + \j/ 2 ), 



and 



x 2 = s\n(yj(k/m)t + i fry) — sin (yj{3k/m)t + ip 2 ). 

Substituting the initial conditions x^O) — X and x 2 (0) = 0 gives 
X = sin t/fj 4- sin ip 2 . 



and 



0 = sin i/f, — sin i j/ 2f 



that is, 



sin (//j = sin i p 2 = X/2. 

The remaining conditions give cos = cos i jj 2 = 0. 
Hence 



Xj = (A'/2)cos y/(k/m)t + {X /2) cos yj(3k/m)t. 



and 



x 2 = <X/2)cos ^{k/nijt — (X/2)cos s /(3k/m)t. 

That is, both natural frequencies are excited and the motion of each body has two 
harmonic components. 
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3.13 Coordinate coupling 

In some systems the motion is such that the coordinates are coupled in the equations 
of motion. Consider the system shown in Fig. 3.1(e); only motion in the plane of the 
figure is considered, horizontal motion being neglected because the lateral stiffness 
of the springs is assumed to be negligible. The coordinates of rotation, 6, and 
translation, x, are coupled as shown in Fig. 3.4. G is the centre of mass of the rigid 
beam of mass m and moment of inertia I about G. 




Fig. 3.4. Two degree of freedom model, rotation plus translation. 

The FBDs are shown in Fig. 3.5; since the weight of the beam is supported by the 
springs, both the initial spring forces and the beam weight may be omitted. 




Fig. 3.5. (a) Applied forces, (b) Effective force and moment. 



For small amplitudes of oscillation (so that sin 9 ~ 9) the equations of motion are 
mx = — fc x (x — Lj0) — k 2 {x + L 2 6) 

and 



l9 = fcj(x — L 1 9)L 1 — k 2 (x + L 2 9)L 2 , 



that is. 



mx + (k t + k 2 )x — ( k 1 L l — k 2 L 2 )9 = 0 
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and 

I& - (k 1 L 1 - k 2 L 2 )x + {^L] + k 2 L\)e = 0. 

It will be noticed that these equations can be uncoupled by making k i L l = k 2 L 2 , 
if this is arranged, translation (x motion) and rotation (0 motion) can take place 
independently. Otherwise translation and rotation occur simultaneously. 

Assuming x t = A 2 sin(cut + i//) and 0 = A 2 sin (cut + i j/), substituting into the 
equations of motion gives 

— mcoMj 4- (k, + k 2 )Aj — (k 1 L 1 — k 2 L 2 )A 2 = 0, 
and 

— Ia> 2 A 2 - (kjL, - k 2 L 2 )A , + (kjLj + k 2 L 2 2 )A 2 = 0. 

That is, 

A l (k i + k 2 — ma> 2 ) + A 2 ( — (k 1 L 1 — k 2 L 2 )) = 0, 
and 

A l (-(k 1 L 1 - k 2 L 2 )) + A 2 {k 2 L\ + k 2 L 2 2 - Ico 2 ) = 0. 

Hence the frequency equation is 

ki + k 2 — mw 2 —(k^L^ — k 2 L 2 ) 

— (kjL, - k 2 L 2 ) kjLi + k 2 L 2 2 — Ia> 2 = ' 

For each natural frequency, there is a corresponding mode shape, given by A x /A 2 . 



Example 18 

A system is modelled as a straight link AB lm long of mass 10 kg, supported 
horizontally by springs A and B of stiffness k l and k 2 respectively. The moment of 
inertia of AB about its centre of mass G is 1 kg m 2 , and G is located at a distance a 
from A and b from B. 

Find the relationship between k lt k 2 , a, and b so that one mode of free vibration 
shall be translational motion only and the other mode rotation only. 

If a = 0.3 m and k x = 13kN/m, find k 2 to give these two modes of vibration and 
calculate the two natural frequencies. 

The frequency equation is (section 3.1.3) 

k t + k 2 — mco 2 — (kjfl — k 2 b) 

— {k^ — k 2 b) k^a 2 + k 2 b 2 — Ia> 2 

For the modes to be uncoupled k y a = k 2 b and then 
cOj = V(( fe i + k: 2 )/m), (bouncing mode) 

and 

a> 2 = yjiik^a 2 + k 2 b 2 )/I) (rocking mode). 

Since a = 0.3 m, b = 0.7 m and k x = 13 kN/m, 
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k 2 = (13 x 0.3)/0.7 = 5.57 kN/m and then 
w, = ^(18570/10) = 43.1 rad/s; so f x = 6.86 Hz. 
and io 2 = 7((13 x 0.3 2 + 5.57 x 0.7 2 )10 3 ) = 62.45 rad/s; so f 2 = 9.94 Hz. 

Example 19 

When transported, a space vehicle is supported in a horizontal position by two 
springs, as shown. The vehicle can be considered to be a rigid body of mass m and 
radius of gyration h about an axis normal to the plane of the figure through the mass 
centre G. The rear support has a stiffness k 1 and is at a distance a from G while the 
front support has a stiffness k 2 and is at a distance b from G. The only motions 
possible for the vehicle are vertical translation and rotation in the vertical plane. 

Write the equations of small amplitude motion of the vehicle and obtain the 
frequency equation in terms of the given parameters. 

Given that k y a — k 2 b, determine the natural frequencies of the free vibrations of 
the vehicle and sketch the corresponding modes of vibration. Also state or sketch 
the modes of vibration if k r a # k 2 b. 




The FBDs are as below: 
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The equations of motion are 

k^y + ad ) + k 2 (y — bd) = —my, 

and 

k t (y + ad)a — k 2 (y — bd)b = —mh 2 d. 
Assuming 

y = Y sin vt and d = © sin vf, 
these give 

Y(k 1 + k 2 — mv 2 ) + ®(k t a — k 2 b) = 0, 

and 



Tf/cjO — k 2 b ) + ©(/qa 2 + k 2 b 2 — mh 2 \ 2 ) = 0. 

The frequency equation is, therefore, 

(k l + k 2 — mv 2 )(/cja 2 + k 2 b z — mh 2 v 2 ) — — k 2 b) 2 — 0. 

If fcjfl = k 2 b, motion is uncoupled so 

% + k 2 



m 



^ rad/s and v 2 = J ra d/s. 



Vj is the frequency of a bouncing or translation mode (no rotation): 



v 2 is the frequency of a rotation mode (no bounce) 
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If k x a ^ k 2 b, the modes are coupled: 




Example 20 

In a study of earthquakes, a building is idealized as a rigid body of mass M supported 
on two springs, one giving translational stiffness k and the other rotational stiffness 
k T , as shown. 

Given that I G is the mass moment of inertia of the building about its mass centre 
G, write down the equations of motion using coordinates x for the translation from 
the equilibrium position, and 9 for the rotation of the building. 

Hence determine the frequency equation of the motion. 
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The FBDs are as follows. 



m 

k T e 

m kX 

Assume small 0 (earthquakes), hence 
m(x 4- h&) = — kx, 

and 

l a 9 + m(x + h&)h = — k r 6 + mghd. 

The equations of motion are therefore 
mhS + mx + kx = 0, 

and 

mh'x + ( mh 2 + I G )S — ( mgh — k r )9 = 0. 

If 

0 = A 1 sino)f and x = A 2 sin cot, 

— mhco 2 A 1 — mco 2 A 2 + kA 2 = 0, 

and 

+ mha> 2 A 2 + (mh 2 4- I G )a) 2 A l + (mgh — k J )A i = 0. 

The frequency equation is 

— mho 2 k — moj 2 _ ^ 

(mh 2 + I G )a> 2 4- (mgh — k r ) mho} 2 

That is, (mha> 2 ) 2 4 - (k — maj 2 )[(mh 2 + l G )w 2 4- (mgh — k T )] = 0, or, 

mI G o ) 4 — a> 2 \_mkh 2 4- I G k — m 2 gh 4- mkj] — mghk 4- kk T = 0. 
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3.1.4 Forced vibration 

Harmonic excitation of vibration in a system may be generated in a number of ways, 
for example by unbalanced rotating or reciprocating machinery, or it may arise from 
periodic excitation containing a troublesome harmonic component. 

A two degree of freedom model of a dynamic system excited by an harmonic force 
F sin vt is shown in Fig. 3.6. Damping is assumed to be negligible. The force has a 
constant amplitude F and a frequency v/2 n Hz. 



+ve +ve 





where 
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A — (k 2 + k — m 2 v 2 )(k l + k — m Y v z ) — k 2 , 

and 

A = 0 is the frequency equation. 

Hence the response of the system to the exciting force is determined. 

Example 21 

A two-wheel trailer is drawn over an undulating surface in such a way that the vertical 
motion of the tyre may be regarded as sinusoidal, the pitch of the undulations being 
5 m. The combined stiffness of the tyres is 170kN/m and that of the main springs is 
60kN/m; the axle and attached parts have a mass of 400 kg, and the mass of the 
body is 500 kg. Find (a) the critical speeds of the trailer in km/h, and (b) the amplitude 
of the trailer body vibration if the trailer is drawn at 50 km/h, and the amplitude of 
the undulations is 0.1 m. 




The equations of motion are 
mjx, = - x 2 ), 

and 

m 2 x 2 = k,(x„ - x 2 ) - k 2 (x 2 - x 3 ). 

Assuming x t = A t sin vt, x 2 = A 2 sin vr, and x 3 = A 3 sin vr, 
A 1 (k i — mjv 2 ) + A 2 ( — kj - 0 

and 

+ A 2 ( k i ~ k 2 m 2 v 2 ) = k 2 A 3 . 
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The frequency equation is 

(k 1 +k 2 - m z v 2 )(k 1 - m x v 2 ) - k\ = 0. 

The critical speeds are those which correspond to the natural frequencies and hence 
excite resonances. The frequency equation simplifies to 

m,m 2 v 4 — ( m 1 k 1 + m l k 2 + m 2 ky)v 2 + k t k 2 = 0. 

Hence substituting the given data, 

500 x 400 x v 4 - (500 x 60 + 500 x 170 + 400 x 60) 10 3 v 2 + 60 x 170 x 10 6 = 0. 

That is 0.2v 4 — 139v 2 + 10200 = 0, which can be solved by the formula. Thus 
v = 16.3 rad/s or 20.78 rad/s, and / = 2.59 Hz or 3.3 Hz. 

Now if the trailer is drawn at v km/h, or v/3.6 m/s, the frequency is i’/(3.6 x 5)Hz. 
Therefore the critical speeds are 

v 1 = 18 x 2.59 = 46.6 km/h, 

and 



v 2 = 18 x3.3 = 59.4 km/h. 

Towing the trailer at either of these speeds will excite a resonance in the system. 
From the equations of motion, 

Mi \ A 

(*i + k 2 - m 2 v 2 )(ky - m^ 2 ) - k\) 3 ’ 

10200 

0.2v 4 - 139v 2 + 10200 

At 50 km/h, v = 17.49 rad/s. 

Thus Ay = —0.749 A 3 . Since A 3 = 0.1 m, the amplitude of the trailer vibration is 
0.075 m. This motion is n out of phase with the road undulations. 

3.1.5 The undamped dynamic vibration absorber 

If a single degree of freedom system or mode of a multi-degree of freedom system is 
excited into resonance, large amplitudes of vibration result with accompanying high 
dynamic stresses and noise and fatigue problems. In most mechanical systems this is 
not acceptable. 

If neither the excitation frequency nor the natural frequency can conveniently be 
altered, this resonance condition can often be successfully controlled by adding a 
further single degree of freedom system. Consider the model of the system shown in 
Fig. 3.7, where K and M are the effective stiffness and mass of the primary system 
when vibrating in the troublesome mode. 

The absorber is represented by the system with parameters k and m. From section 
3.1.4 it can be seen that the equations of motion are 

MX = — KX — k(X — x) + F sin vt, for the primary system 
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and 

mx - k(X — x), 

for the vibration absorber. Substituting 
X = X 0 sin vt and x = x 0 sin vf 

gives 

X 0 (K + k - Mv 2 ) + x 0 {-k) = F, 

and 

X 0 ( — k) + x 0 (k — mv 2 ) — 0. 

Thus 

v F(k - mv 2 ) 

A°- — , 

and 

Fk 

where 

A = (k — mv 2 \K + k — Mv 2 ) — k 2 , 
and A = 0 is the frequency equation. 

It can be seen that not only does the system now possess two natural frequencies, 
Qj and Q 2 instead of one, but by arranging for k — mv 2 = 0, X 0 can be made zero. 
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Thus if yj(k/m ) = yJ{K/M), the response of the primary system at its original 
resonance frequency can be made zero. This is the usual tuning arrangement for an 
undamped absorber because the resonance problem in the primary system is only 
severe when v ~ y/(K/M) rad/s. This is shown in Fig. 3.8. 




w =/*/*#) = 



Fig. 3.8. Amplitude-frequency response for system with and without tuned absorber. 



When X 0 = 0, x 0 = —F/k, so that the force in the absorber spring, kx 0 is — F; 
thus the absorber applies a force to the primary system which is equal and opposite 
to the exciting force. Hence the body in the primary system has a net zero exciting 
force acting on it and therefore zero vibration amplitude. 

If an absorber is correctly tuned to 2 = K/M = k/m, and if the mass ratio /i = m/M, 
the frequency equation A = 0 is 




— (2 + n) 




+ 1 = 0 . 



This is a quadratic equation in (v/cy) 2 . Hence 
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and the natural frequencies and Q 2 are found to be 




For a small fi, Qj and H 2 are very close to each other, and near to w; increasing 
/i gives better separation between and as shown in Fig. 3.9. 




Fig. 3.9. Effect of absorber mass ratio on natural frequencies. 



This effect is of great importance in those systems where the excitation frequency 
may vary; if n is small, resonances at or Q 2 may be excited. It should be noted 
that since 




and 




Then 



nln 



2 

2 



CO 



4 





= 1 . 



2 
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That is, 

Qj Q 2 = a)2 - 

Also 





= 2 + n- 



These relationships are very useful when designing absorbers. If the proximity of 
Qj and fi 2 to a> is likely to be a hazard, damping can be added in parallel with the 
absorber spring, to limit the response at these frequencies. Unfortunately, if damping 
is added, the response at frequency a> will no longer be zero. A design criterion that 
has to be carefully considered is the possible fatigue and failure of the absorber spring: 
this could have severe consequences. In view of this, some damped absorber systems 
dispense with the absorber spring and sacrifice some of the absorber effectiveness. 
This is particularly important in torsional systems, where the device is known as a 
Lanchester damper. 



Example 22 

The figure represents a pump of mass which rests on springs having a stiffness 
fcj, so that only vertical motion can occur. Given that the damping is negligible and 
the mass m 2 is ignored, derive an expression for the frequency of the harmonic 
disturbing force at which the pump will execute vertical oscillations of very 
large — theoretically infinite — amplitude. 

Given that an undamped dynamic absorber of mass m 2 is then connected to the 
pump by a spring of stiffness k 2 , as shown, prove that the amplitude of the oscillations 
of the pump is reduced to zero when 



m 2 



where v is the natural frequency of the free vibrations of the pump in the absence of 
the dynamic vibration absorber. 

The pump has a mass of 1 30 kg and rotates at a constant speed of 2400 revolutions 
per minute but due to a rotating unbalance very large amplitudes of pump vibration 
on the spring supports result. An undamped vibration absorber is to be fitted so that 
the nearest natural frequency of the system is at least 20 per cent removed from the 
running speed of 2400 revolutions per minute. Find the smallest absorber mass 
necessary and the corresponding spring stiffness. 




The equation of motion is 



mjX 

so that if 
*1 = 



j 4- = Fsin vr. 



X x sin vt, 



*1 ~ 



F 

k t — mjV 2 
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When 




That is, resonance occurs when v = u> = y l{k l lm l ). With a vibration absorber added, 
the system is 




The FBDs are therefore, if x 2 > is assumed, 



M 4 

I m 2 

k 2 (x 2 - x ,) 

Fsinvf k 2^ x 2~ x \) 








m 1 






t 

k 1 X 1 




m 2 


X 2 



The equations of motion are thus 
m 2 x 2 = -k 2 (x 2 - Xj), 
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or 

m 2 x 2 + k 2 x 2 — k 2 x l - 0, 

and 

m 1 x 1 = k 2 (x 2 — Xj) — kyX j + F sin vt, 
or 

m l x l + (/cj + k 2 )x j — k 2 x 2 = F sin vt. 

Assuming x, = X 2 sin vt and x 2 = X 2 sin vt, these equations give 
*i[*i + k 2 - mjv 2 ] + X 2 [-/c 2 ] = F, 
and 

*i[-k 2 ] + X 2 [k 2 - w 2 v2 ] = 0. 

That is 

x = f(fc 2 ~ m 2 y2 ) 

1 [(^1 + fc 2 ) - "V’ 2 ][fc 2 - m 2 v2 ] - k 2 

Thus if 




Now the frequency equation is 

[(^i + k 2 ) - m l v 2 '][k 2 - m 2 v 2 ] -k\ = 0. 
If we put 




this becomes 

v 4 - v 2 (/iQ 2 + 2Q 2 ) + Q 4 = 0, 
or 

©‘-©e + * + '- 0 

so that 

The limiting condition for the smallest absorber mass is (Vj/Q) = 0.8 because then 
(v 2 /Q) = 1.25, which is acceptable. Thus 
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+ 2\ 2 n 2 + 4/t 

~Y~) ~~ 4 ~ 

and 

n = 0 . 2 . 

Hence 

m 2 = 0.2 x 130 = 26 kg, 

and 

k 2 = (807r) 2 m — 1642 kN/m. 




Example 23 

A system has a violent resonance at 79 Hz. As a trial remedy a vibration absorber 
is attached which results in a resonance frequency of 65 Hz. How many such absorbers 
are required if no resonance is to occur between 60 and 120 Hz? 



Since 



— Y + (— Y =2+ & 



and 



= 0)2 

in the case of one absorber, with to = 79 Hz and = 65 Hz, 
79 2 

Q 2 = — = 96 Hz. 

65 



Also 




so n — 0.154. 

In the case of n absorbers, if 

79 2 

flj = 60 Hz, Q 2 = — - = 104 Hz (too low). 
60 



So require Q 2 = 120 Hz and then Qj = (79 2 /120) = 52 Hz. Hence 




Thus 

0.74 

n = 0.74 - n.n and n = — = 4.82. 



Thus five absorbers are required. 
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Example 24 

A machine tool of mass 3000 kg has a large resonance vibration in the vertical direction 
at 120 Hz. To control this resonance, an undamped vibration absorber of mass 600 kg 
is fitted tuned to 120 Hz. Find the frequency range in which the amplitude of the 
machine vibration is less with the absorber fitted than without. 

If (X 0 ) with absorber = (A 0 ) without absorber, 



F(k — mv 2 ) 

(K + k — Mv 2 )(k — mv 2 ) — k 2 




(phase requires — ve sign). 



Multiplying out and putting jj. = m/M gives 

2 £)‘- <4 + ^© i + 2 = 0 ' 

Since 



600 

3000 



= 0 . 2 , 




4 + n 
4 



± ivV + M = 1.05 ± 0.32. 



Thus 



- = 1.17 or 0.855, 
co 

and 

/ = 102 Hz or 140 Hz, where v = 2nf. 
Thus the required frequency range is 102—140 Hz. 



3.1.6 System with viscous damping 

If a system possesses damping of a viscous nature, the damping can be modelled 
similarly to that in the system shown in Fig. 3.10. 



I— 



PMWH 



m . 



n 



HVyA/VWWWH 



n 




Fig. 3.10. Two degree of freedom viscous damped model with forced excitation. 
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For this system the equations of motion are 

+ k 1 x 1 + k 2 (x l - x 2 ) + qxj + c 2 (x 1 - x 2 ) =/ lf 

and 

m 2 x 2 + k 2 (x 2 - xj + k 3 x 2 + c 2 (x 2 - x,) + c 3 x 2 =f 2 . 

Solutions of the form x } = >t,e s! and x 2 = A 2 e st can be assumed, where the Laplace 
operator s is equal to a + jb, j = v /(— 1), and a and b are real— that is each solution 
contains an harmonic component of frequency b, and a vibration decay component 
of damping factor a. By substituting these solutions into the equations of motion a 
frequency equation of the form 

s 4 + as 3 + /fa 2 + ys + 5 = 0 

can be deduced, where a, /5, y, and S are real coefficients. From this equation four 
roots and thus four values of s can be obtained. In general the roots form two complex 
conjugate pairs such as a { ± j b x , and a 2 ± }b 2 . These represent solutions of the form 
x = Ref-Xe 0 ' . e jfcI ) = Xe ttt cos bt. That is, the motion of the bodies is harmonic, and 
decays exponentially with time. The parameters of the system determine the magnitude 
of the frequency and the decay rate. 

It is often convenient to plot these roots on a complex plane as shown in Fig. 3.11. 
This is known as the s-plane. 



imls! 




Fig. 3.11. s-plane. 
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For light damping the damped frequency for each mode is approximately equal 
to the undamped frequency, that is, b x — co x and b 2 ^o) 2 . 

The right-hand side of the s-plane (Re(s) -I- ve) represents a root with a positive 
exponent, that is, a term which grows with time, so unstable motion may exist. 

The left-hand side contains roots with a negative exponent so stable motion exists. 
See also Fig. 2.20. 

All passive systems have negative real parts and are therefore stable but some 
systems such as rolling wheels and rockets can be become unstable, and thus it is 
important that the stability of a system is considered. This can be conveniently done 
by plotting the roots of the frequency equation on the s-plane, as above. This technique 
has particular application in the study of control system stability. 



3.2 THE VIBRATION OF SYSTEMS WITH MORE THAN TWO DEGREES 
OF FREEDOM 

The vibration analysis of a dynamic system with three or more degrees of freedom 
can be carried out in the same way as the analysis given above for two degrees of 
freedom. However, the method becomes tedious for many degrees of freedom, and 
numerical methods may have to be used to solve the frequency equation. A computer 
can, of course, be used to solve the frequency equation and determine the corre- 
sponding mode shapes. Although computational and computer techniques are 
extensively used in the analysis of multi-degree of freedom systems, it is essential for 
the analytical and numerical bases of any program used to be understood, to ensure 
its relevance to the problem considered, and that the program does not introduce 
unacceptable approximations and calculation errors. For this reason it is necessary 
to derive the basic theory and equations for multi-degree of freedom systems. 
Computational techniques are essential, and widely used, for the analysis of the 
sophisticated structural models often devised and considered necessary, and computer 
packages are available for routine analyses. However, considerable economies in 
writing the analysis and performing the computations can be achieved, by adopting 
a matrix method for the analysis. Alternatively an energy solution can be obtained 
by using the Lagrange equation, or some simplification in the analysis achieved by 
using the receptance technique. The matrix method will be considered first. 



3.2.1 The matrix method 

The matrix method for analysis is a convenient way of handling several equations 
of motion. Furthermore, specific information about a system, such as its lowest natural 
frequency, can be obtained without carrying out a complete and detailed analysis. 
The matrix method of analysis is particularly important because it forms the basis 
of many computer solutions to vibration problems. The method can best be 
demonstrated by means of an example. For a full description of the matrix method 
see Mechanical Vibrations : Introduction to Matrix Methods, by J. M. Prentis & F. A. 
Leckie (Longmans, 1963). 
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Example 25 

A structure is modelled by the three degree of freedom system shown. Determine the 
highest natural frequency of free vibration and the associated mode shape. 
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The equations of motion are 

2 mx 1 + 2&Xj 4- fc(xj — x 2 ) = 0, 

2 mx 2 + k(x 2 — x,) + k(x 2 — x 3 ) = 0, and 



mx 3 + k{x 3 — x 2 ) = 0. 

If Xj, x 2 and x 3 take the form X sin art and ). = mio 2 /k, these equations can be written 

l x i ~ \ X 2 =^i, 

- \X l + X 2 - jX 3 = /.X 2 , and 

- x 2 + x 3 = ;.a 3 . 

that is, 




or 

[S]{x} = ;.{x} 

where [S] is the system matrix, {X} is a column matrix, and the factor A is a scalar 
quantity. 

This matrix equation can be solved by an iteration procedure. This procedure is 
started by assuming a set of deflections for the column matrix {X} and multiplying 
by [S]; this results in a new column matrix. This matrix is normalized by making 
one of the amplitudes unity and dividing each term in the column by the particular 
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amplitude which was put equal to unity. The procedure is repeated until the amplitudes 
stabilize to a definite pattern. Convergence is always to the highest value of X and 
its associated column matrix. Since /. = ma> 2 k, this means that the highest natural 
frequency is found. Thus to start the iteration a reasonable assumed mode would be 



Now 





1.5 


-0.5 


0 


f o 


-0.5 
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-0.5 


- 1 


_ 0 


-1 


1 _ 


l 2 J 



Using this new column matrix gives 




= 3 



0.67') 
— 0.83 2 

i.ooJ 



1.5 -0.5 0 

-0.5 1 -0.5 

_ 0 -1 1 _ 



f 0.67] f 1.415] 

< -0.83 > = \ — 1.665 > 

l l.OOj L 1.83 3 




and eventually, by repeating the process the following is obtained: 
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-0.5 
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f 1 
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-0.5 
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l lJ 


1 l 1 3 



Hence X — 2 and <x> 2 = 2k/m. X is an eigenvalue of [S], and the associated value of 
{X} is the corresponding eigenvector of [S]. The eigenvector gives the mode shape. 

Thus the highest natural frequency is 1/2 7t yj(2k/m) Hz, and the associated mode 
shape is 1: — 1:1. Thus if X 1 — 1, X 2 = — 1 and X 3 = 1. 

If the lowest natural frequency is required, it can be found from the lowest eigenvalue. 
This can be obtained directly by inverting [S] and premultiplying [S]{X} — /{X} by 

;. -1 [S] -1 . 

Thus [S] '{X} = /. '{X}. Iteration of this equation yields the largest value of X 1 
and hence the lowest natural frequency. A reasonable assumed mode for the first 
iteration would be 




Alternatively, the lowest eigenvalue can be found from the flexibility matrix. The 
flexibility matrix is written in terms of the influence coefficients. The influence 
coefficient a pq of a system is the deflection (or rotation) at the point p due to a unit 
force (or moment) applied at a point q. Thus, since the force each body applies is the 
product of its mass and acceleration. 

X 1 = a,, 2mX 1 a> 2 + a 12 2mX 2 io 2 + ot 13 mX 3 a> 2 , 

X 2 — a 21 2mX l co 2 + a 22 2 mX 2 o 2 4- a 23 mX 3 a> 2 , 
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and 



X 3 = a 31 2mA' 1 c/j 2 + a 32 2mX 2 oj 2 + oc 33 mX 3 a> 2 . 
or 

2a n 2a 12 
2a 21 2a 22 

_2a 31 2a 32 

The influence coefficients are calculated by applying a unit force or moment to 
each body in turn. Since the same unit force acts between the support and its point 
of application, the displacement of the point of application of the force is the sum of 
the extensions of the springs extended. The displacements of all points beyond the 
point of application of the force are the same. 

Thus 




1 

*11 = *12 =*13 = *21 = *31 =2fc> 

and 

1 1 1 _ 5 

Ci33 ~2k + k + k~2k' 

Iteration causes the eigenvalue k/mco 2 to converge to its highest value, and hence 
the lowest natural frequency is found. The other natural frequencies of the system 
can be found by applying the orthogonality relation between the principal modes of 
vibration. 



3.2. 1.1 Orthogonality of the principal modes of vibration 

Consider a linear elastic system that has n degrees of freedom, n natural frequencies, 
and n principal modes. 

The orthogonality relation between the principal modes of vibration for an n degree 
of freedom system is 

X 'M i (r)/l i (s) = 0, 

i = 1 

where A^r) are the amplitudes corresponding to the rth mode, and A ; (s) are the 
amplitudes corresponding to the sth mode. 

This relationship is used to sweep unwanted modes from the system matrix, as 
illustrated in the following example. 

Example 26 

Consider the three degree of freedom model of a system shown. 




Sec. 3.2] 



The vibration of systems with more than two degrees of freedom 119 











■ I 


+ve 


■ 




mm 




91 










. I 






4m 




2m 




m 



p . - a 



a " a 



The equations of motion in terms of the influence coefficients are 
= 4oc ll mX l to 2 + 2 oc 12 mX 2 a} 2 + <x i3 mX i co 2 , 

X 2 — 4ot 2l mX 1 u ) 2 + 2a 22 mX 2 (o 2 + x 23 mX 3 u) 2 , 

and 

X 3 = 4cc 3l mX 3 a> 2 + 2ot 32 mX 2 u) 2 + oc 33 mX 3 oj 2 . 

That is, 



Now, 
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and 



^ = Tk- 



Hence, 



l ^ 1 ( co 2 m 

= ^k 



4 2 1 

4 8 4 
L4 8 7. 



I, 



To start the iteration a reasonable estimate for the first mode is 

b 

this is inversely proportional to the mass ratio of the bodies. 
Eventually iteration for the first mode gives 
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or u> 1 = 0A6yJ(k/m) rad/s. 

To obtain the second principal mode, use the orthogonality relation to remove the 
first mode from the system matrix: 



m l A l A 2 + m 2 B x B 2 + m 3 C l C 2 — 0. 

Thus 



4m(1.0)/l 2 4- 2m(3.18)B 2 + m(4.0)C 2 — 0, 
or A 2 = — 1.59 B 2 — C 2 , since the first mode is 




Hence, rounding 1.59 up to 1.6, 




When this sweeping matrix is combined with the original matrix equation, iteration 
makes convergence to the second mode take place because the first mode is swept 
out. Thus, 




Now estimate the second mode as 




Hence oj 2 = ^J(k/m) rad/s, and the second mode was evidently estimated correctly 
as 1 :0: — 1. 

To obtain the third mode, write the orthogonality relation as 
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m 1 A 2 A 3 + m 2 B 2 B 3 + m 3 C 2 C 3 = 0, 

and 

m l A l A 3 + m 2 B x B 3 + m 3 C x C 3 = 0. 

Substitute 

A x = 1.0, B x = 3.18, C, = 4.0, 
and 

A 2 = 1.0, B 2 = 0, c 2 =- 1.0, 

as found above. Hence 

(a 3 ) To o 0.25 I r^ 3 'j 

<B 3 >= 0 0 -0.78 <B 3 >. 

Lc 3 J Lo o i J tc 3 J 

When this sweeping matrix is combined with the equation for the second mode 
the second mode is removed, so that it yields the third mode on iteration: 




An estimate for the third mode shape now has to be made and the iteration 
procedure carried out once more. In this way the third mode eigenvector is found to be 




and a> 3 = 1.32 yj(k/m) rad/s. 

The convergence for higher modes becomes more critical if impurities and 
rounding-off errors are introduced by using the sweeping matrices. It is well to check 
the highest mode by the inversion of the original matrix equation, which should be 
equal to the equation formulated in terms of the stiffness influence coefficients. 

3.2.2 The Lagrange equation 

Consideration of the energy in a dynamic system together with the use of the Lagrange 
equation is a very powerful method of analysis for certain physically complex systems. 
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It is an energy method which allows the equations of motion to be written in terms 
of any set of generalized coordinates. Generalized coordinates are a set of independent 
parameters which completely specify the system location and which are independent 
of any constraints. The fundamental form of Lagrange’s equation can be written in 
terms of the generalized coordinates q t as follows: 

<j/d(7Tv c(T) | d(V) | d(DE) 
d AcqJ dq { dq, dq t ^ 

where T is the total kinetic energy of the system, Fis the total potential energy of 
the system, DE is the energy dissipation function when the damping linear (it is half 
the rate at which energy is dissipated so that for viscous damping DE = \cx 2 X Qi is 
a generalized external force (or non-linear damping force) acting on the system, and 
q t is a generalized coordinate that describes the position of the system. 

The subscript i denotes n equations for an n degree of freedom system, so that the 
Lagrange equation yields as many equations of motion as there are degrees of freedom. 
For a free conservative system Q t and DE are both zero, so that 

d fd(T)\ d(T) | d(V) _ Q 
dr V 8q, ) £q t cq i 

The full derivation of the Lagrange equation can be found in Vibration Theory 
and Applications by W. T. Thomson (Allen & Unwin, 1989). 

Example 27 

A solid cylinder has a mass M and radius R. Pinned to the axis of the cylinder is an 
arm of length l which carries a bob of mass m. Obtain the natural frequency of free 
vibration of the bob. The cylinder is free to roll on the fixed horizontal surface shown. 

The generalized coordinates are x 1 and x 2 . They completely specify the position 
of the system and are independent of any constraints. 

T = \Mx \ + &MR 2 )& 2 + \mx\ 

= \Mx\ + &Mxi) + \mx\. 

V = mg 1(1 - cos <f>) = (mgl/2)<f) 2 = (mg/2l)(x 2 - x t )\ 

for small values of <£. Apply the Lagrange equation with q f = x l : 

(d/drxer/dxj) = Mx l 4- \Mx j 

dV/dx l ~ (mg/2l)( — 2x 2 + 2x 1 ). 

Hence f Mx 1 + ( mg/l)(x l — x 2 ) — 0 is an equation of motion. 

Apply the Lagrange equation with q t = x 2 : 

(d/dt)(ST/dx 2 ) = mx 2 

dV/dx 2 = ( mg/2l)(2x 2 — 2x x ). 



Hence 




Sec. 3.2] 



The vibration of systems with more than two degrees of freedom 1 23 




mx 2 + (mg/l)(x 2 — xj = 0 is an equation of motion. 

These equations of motion can be solved by assuming that Xj = X t sin cot and 
x 2 — X 2 sintof. Then 

Xi((mg/t) ~ (3M/2)co 2 ) + X 2 (-mg/l) = 0, 

and X^ — mg/l) + X 2 ((mg/l) — mo) 2 ) — 0. 

The frequency equation is therefore 

(3M/2)o> 4 - (g/l)co 2 (m + (3M/2)) = 0. 

Thus either to = 0, or to = ^/((l + 2m/3M)g/l) rad/s, and XJX 2 = —2m/3M. 

Example 28 

Find the equations of motion for free vibration of the triple pendulum shown overleaf. 

T = + \m 2 L 2 & 2 2 + \m 3 l}b\ 

V = m^gL{ 1 — cos + m 2 gL( 1 — cos d 2 ) + m 3 gL( 1 — cos 9 3 ) 

+ sin 6 2 — a sin 9 J 2 + jk 2 (a sin 9 3 — a sin d 2 ) 2 

Hence, by applying the Lagrange equation with qt — 9 1 ,9 2 ,9 3 in turn, the equations 
of motion are obtained. Assume 6 is small so that sin 9^9, and cos 9 ~ 1 — 9 2 / 2. 
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Example 29 

To isolate a structure from the vibration generated by a machine, the machine is 
mounted on a large block. The block is supported on springs as shown. Find the 
equations which describe the motion of the block in the plane of the figure. 
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The coordinates used to describe the motion are q x ,q 2 , and q 3 . These are generalized 
coordinates because they completely specify the position of the system and arq 
independent of any constraints. If the mass of the block and machine is M, and the 
total mass moment of inertia about G is / G , then 

T = \Mq\ + \Mq\ + \l G q\, and 

V = strain energy stored in the springs 

= ^k 1 {q l + bq 3 ) 2 + { k l (q 1 - d q 3 ) 2 

2^2^q 2 — a q$) "T 2k 2 (q 2 ~ a q^) ■ 

Now apply the Lagrange equation with q t = q x . 



and 



cq i 

cT 

— = Mq lf 

cqi 



d 

so — x 
dr 



cT 
— - 

q i 



cV 

— = kMi + bq 3 ) + k 1 (q l - dq 3 ). 
cq x 

Thus the first equation of motion is 



Mq x + 2k x q x + k x (b — d)q 3 = 0. 

Similarly by putting q t = q 2 and q t = q 3 , the other equations of motion are obtained as 
Mq 2 + lk x q 2 — lak 2 q 3 = 0, 

and 



(b — d)q x — 2ak 2 q 2 + (( b 2 + d 2 )k x + 2a 2 k 2 )q 3 = 0. 

The system therefore has three coordinate-coupled equations of motion. The natural 
frequencies can be found by substituting q x = A t sin cor, and solving the resulting 
frequency equation. It is usually desirable to have all natural frequencies low so that 
the transmissibility is small throughout the range of frequencies excited. 



3.2.3 Receptances 

Some simplification in the analysis of multi-degree of freedom undamped dynamic 
systems can often be gained by using receptances, particularly if only the natural 
frequencies are required. If an harmonic force F sin vr acts at some point in a system 
so that the system responds at frequency v, and the point of application of the force 
has a displacement x = X sin vt, then if the equations of motion are linear, 
x = xF sin vr, where x, which is a function of the system parameters and v but not 
a function of F, is known as the direct receptance at x. If the displacement is determined 
at some point other than that at which the force is applied, x is known as the transfer 
or cross receptance. 
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It can be seen that the frequency at which a receptance becomes infinite is a natural 
frequency of the system. Receptances can be written for rotational and translational 
coordinates in a system, that is the slope and deflection at a point. 

Thus, if a body of mass m is subjected to a force F sin vf and the response of the 
body is x = X sin vf, 

F sin vf = mx = m( — Xv 2 sin vf) = — mv 2 x. 

Thus 

1 

x = ~ F sin vf, 

mv 



and 



1 

a = 2- 

mv 

This is the direct receptance of a rigid body. 

For a spring, a = l/k. This is the direct receptance of a spring. 

In an undamped single degree of freedom model of a system, the equation of motion 
is 



mx + kx = F sin vf, 

If x — X sin vf, a = 1 /(k — mv 2 ). This is the direct receptance of a single degree of 
freedom system. 

In more complicated systems, it is necessary to be able to distinguish between 
direct and cross receptances and to specify the points at which the receptances are 
calculated. This is done by using subscripts. The first subscript indicates the coordinate 
at which the response is measured, and the second indicates that at which the force 
is applied. Thus a pq , which is a cross receptance, is the response at p divided by the 
harmonic force applied at q, and a pp and a qq are direct receptances at p and q 
respectively. 

Consider the two degree of freedom system shown in Fig. 3.12 




HW\Ah m i hWM/WV“ 



Q D 




Fig. 3.12. Two degree of freedom system with forced excitation. 

The equations of motion are 

m 1 x 1 + (fc, + k 2 )x 1 - k 2 x 2 =/j. 



and 
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m 2 x 2 + (k 2 + k 2 )x 2 — k 2 x 1 — 0 . 

Let /j = Fj sin vf, and assume that x x = X x sin vr and x 2 — X 2 sin vt. 
Substituting into the equations of motion gives 

(k { + k 2 - m^v^X , +(-k 2 )X 2 = F x , 

and 

( — k 2 )X l + (k 2 + k 3 — m 2 v 2 )X 2 — 0 . 

Thus 



X i k 2 + k 3 — m 2 v 2 

aii = 7 ; = a 5 

where 

A = (fcj + k 2 — m 1 v 2 Xfe 2 + k 3 — m 2 v 2 ) — /c 2 . 

3t n is a direct receptance, and A = 0 is the frequency equation. 
Also the cross receptance 




This system has two more receptances, the responses due to f 2 applied to the second 
body. Thus a 12 and a 22 may be found. It is a fundamental property that a 12 = a 21 
{principle of reciprocity) so that symmetrical matrices result. 

A general statement of the system response is 

X l = SjjFj + 3«i 2 T 2 , 

and 



X 2 — Xjik l + a 22^2’ 
That is, 




Some simplification in the analysis of complex systems can be achieved by 
considering the complex system to be a number of simple systems (whose receptances 
are known) linked together by using conditions of compatibility and equilibrium. The 
method is to break the complex system down into subsystems and analyse each 
subsystem separately. Find each subsystem receptance at the point where it is 
connected to the adjacent subsystem, and ‘join’ all subsystems together, using the 
conditions of compatibility and equilibrium. 

For example, to find the direct receptance y u of a dynamic system C at a single 
coordinate the system is considered as two subsystems A and B, as shown in 
Fig. 3.13. 
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Ft sin v? 

Xi sin vt 

System C 
receptance y 



|F a sinvf F b sinvt| 

A — ► — ► B 

|X a sin vf X b sin vfj 

Sub-system A Sub-system B 

receptance a receptance 

Fig. 3.13. Dynamic systems. 



By definition 



Because the systems are connected. 
X a = X b = X l , (compatibility) 



ft 

P it — F ■ 



F j = F a -f F h , (equilibrium). 



Hence 

_L-_L _L 

Vn *n Pu 

that is, the system receptance y can be found from the receptances of the subsystems. 

In a simple spring — body system, subsystems A and B are the spring and body 
respectively. Hence a n = 1 jk and j? n — — I/mv 2 and l/y u = k — mv 2 . as above. 
The frequency equation is a n 4- )S n = 0, because this condition makes y n = oo. 

Consider the undamped dynamic vibration absorber application shown in Fig. 3.14. 
The system is split into subsystems A and B. 




Fig. 3.14. Absorber sub systems. 
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For subsystem A, 



(— *A =*A Sin,,f 

U‘ F S' nvt 



/ A = Mx A + Kx a ; hence a = 1 /(K — Mv 2 ). 

For subsystem B, 

H“ X 8 =X B SinW I ^>'B = y B sinW 




f B = k(x B - y B ) = my B = — mv 2 T B ; hence = —(k - mv 2 )/kmv 2 . 

Thus the frequency equation a + /? = 0 gives 
Mmv 4 — ( mK + Mk + mk)v 2 + Kk = 0. 

It is often convenient to solve the frequency equation a + = 0, or a = —[}, by a 

graphical method. In the case of the absorber, both a and — /? can be plotted as a 
function of v, and the intersections give the natural frequencies and Q 2 (Fig. 3.15). 

It can be seen that the effect of adding different absorbers to the primary system 
can readily be determined without re-analysing the whole system. It is merely necessary 
to sketch the receptance of each absorber on Fig. 3.15 to find and fl 2 for the 
complete system. 

The receptance technique is particularly useful when it is required to investigate 
the effects of adding a dynamic system to an existing system, for example an extra 
floor, or an air-conditioning plant to a building. Once the receptance of the original 
system is known, it is only necessary to analyse the additional system, and then to 
include this in the original analysis. Furthermore, sometimes the receptances of 
dynamic systems are measured, and available only in graphical form. Example 30. 
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Some subsystems, such as those shown in Fig. 3.16, are linked by two coordinates, 
for example deflection and slope at the common point. 




System C Subsytem A 

Receptance y Receptance a 



F a1 sin vfF b1 sin vt 
X a1 sin vtXi,, sin vt 
Fa2 sin vt F b 2 sin vt 
Xaz sin vt Xb2 sin vt 




Subsystem B 
Receptance p 



Fig. 3.16. Applied forces and system responses. 

Now in this case, 

-^al = + a l2^a2> 

-^a2 = a 2t-^ al + a 22^a2’ 

-^bl = Pi 1 ^bl + Pi 2^b2 5 

-^b2 = Pzi^bi + P 22 F bl- 



and 
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The applied forces or moments are F 1 sin vf and F 2 sin vf where 
F\ = F a j + F bl , 

and 

^2 = F F b2 - 

Since the subsystems are linked 

and 

* 2 =X a2 = X b2 . 

Hence if excitation is applied at x x only, F 2 — 0 and 

_ -X'l _ a n(^li^22 ~ ^12 2 ) F ^n(^ii g 22 ~ a 12 ) 

yn ~T x ~ A 

where 

A = (^11 + ^uX a 22 F P 22 ) ~~ ( a i 2 "h P 12 )’ 

and 

_ %2 _ a ! 2 (^ 1 1/^22 ~ ^ 12 ^ 12 ) ~ ^12( g li a 22 ~ a i2^12) 

721 ~J~ 1 ~ A 

If 

F X — 0, 

_ %2 _ ' X 22^P 11^22 ~ Pl2 2 ) ~ ^ 22 (^ 11^22 ~ a l2 ) 

722 ~ "" A 

Since A = 0 is the frequency equation, the natural frequencies of the system C are 
given by 

X ll F *\2 F ft 12 _ Q 
a 2 i + P 21 ^22 "h P 22 

This is an extremely useful method for finding the frequency equation of a system 
because only the receptances of the subsystems are required. The receptances of many 
dynamic systems have been published in The Mechanics of Vibration by R. E. D. 
Bishop & D. C. Johnson (CUP 1960/1979). By repeated application of this method, 
a system can be considered to consist of any number of subsystems. This technique 
is, therefore, ideally suited to a computer solution. 

It should be appreciated that although the receptance technique is useful for writing 
the frequency equation, it does not simplify the solution of this equation. 
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Example 30 

An engine with a flywheel is shown. The engine is to drive a propeller of inertia I 
through a shaft of torsional stiffness k as indicated. The receptance for torsional 



Engine Flywheel 




■* n 

-TT-'J’J 

k 



Propeller 

' “\ 

I • 




oscillation of the engine-flywheel system has been measured at point A over a limited 
frequency range which does not include any internal resonances of the system. The 
figure shows the receptance at A as a function of (frequency) 2 . 

Calculate the lowest non-zero natural frequency for the engine, flywheel and 
propeller system if / — 0.9 kg m 2 and k = 300 kN m/rad. 
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The engine — propeller system can be considered to be two subsystems, A and B: 




The frequency equation is a n -f /? X1 = 0, where a n is given in graphical form and 
Pn is found as follows. For the propeller system, 




Now 

h = k(9 B - (f > B ) = 1 $ : B = -/<V’ 2 sin vt, 
so 

k(&B ® b ) = — I®B v2 

and 

n 1 ^ 2 ^ (k — Iv 2 

Ob = $b - l <V - ^b( fc 



Thus 
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0n 



©B 

7b 




1 

<Vv 2 




If a u and — ^ u -are plotted as functions of (frequency) 2 , the intersection gives the 
value of the frequency which is a solution of 2 n +/J n =0, that is, the natural 
frequency of free vibration is found. The table below can be calculated for — /3 n 
because k = 300 x 10 3 Nm/rad, and / = 0.9 kg m 2 . 



V 2 




/v 2 


k~Iv 2 


klv 2 


-fin 


X 

r*l 

© 


10 6 


0.27 x 10 6 


0.03 x 10 6 


0.081 x 10 12 


-0.37 x 10 6 


X 

■'3- 

d 


10 6 


0.36 x 10 6 


0.06 x 10 6 


0.108 x 10 12 


-0.55 x 10“ 6 


0.5 x 


I0 6 


0.45 x 10 6 


0.15 x 10 6 


0.135 x 10 12 


-1.11 x 10" 6 



The receptance /? n can now be plotted against (frequency) 2 as below: 




The intersection occurs at (frequency) 2 = 0.39 x 10 6 (rad/s) 2 , that is, a frequency of 
624 rad/s or 99.4 Hz. This is the natural frequency of the combined engine — propeller 
system. It can be seen that the effect of using different propellers with different k and 
I values is easily found, without having to re-analyse the whole engine — propeller 
system. 
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3.2.4 Impedance and Mobility 

Impedance and mobility analysis techniques are frequently applied to systems and 
structures with many degrees of freedom. However, the method is best introduced 
by considering simple systems initially. 

The impedance of a body is the ratio of the amplitude of the harmonic exciting 
force applied, to the amplitude of the resulting velocity. The mobility is the reciprocal 
of the impedance. It will be appreciated, therefore, that impedance and mobility 
analysis techniques are similar to those used in the receptance analysis of dynamic 
systems. 

For a body of mass m subjected to an harmonic exciting force represented by Fe jv ' 
the resulting motion is x = Xe jvl , Thus 

Fe jM = mx = -mv 2 Xe jw , 
and the receptance of the body, 

X _ 1 

F mv 2 ' 



Now 



Fe jw = -mv 2 Xe jv ' 

= mjv(jvXe jvI ) = mjv v, 

where v is the velocity of the body, and v = Fe jvI . 

Thus the impedance of a body of mass m is Z m , where 




and the mobility of a body of mass m is M m , where 



M 



m 



F jmv 



Putting s = jv so that x = Xe st gives 



z m = m s, 



M m = — , 
ms 



and 



V = sX. 

For a spring of stiffness k , Fe jvt = kXe ivt and thus Z k — F/V— k/s and M k = s/k, 
whereas for a viscous damper of coefficient c, Z c = c and M c = 1/c. 

If these elements of a dynamic system are combined so that the velocity is common 
to all elements, then the impedances may be added to give the system impedance, 
whereas if the force is common to all elements the mobilities may be added. This is 
demonstrated below by considering a spring— mass single degree of freedom system 
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Fig. 3.17. Single degree of freedom system with elements connected in parallel. 



with viscous damping, as shown in Fig. 3.17. 

The velocity of the body is common to all elements, so that the force applied is 
the sum of the forces required for each element. The system impedance, 

^ F F m + F k + F c 

V~ V 



~ Z m + Z k + Z c 



Hence 



Z = ms H 1- 

s 



c. 



That is 



F = (ms 2 + cs + k)X, 



or 



F = (k — mv 2 + }cv)X. 



Hence 



7[(fc - mv 2 ) 2 + (cv) 2 ] ' 

Thus when system elements are connected in parallel their impedances are added 
to give the system impedance. 

In the system shown in Fig. 3.18, however, the force is common to all elements. 
In this case the force on the body is common to all elements so that the velocity at 
the driving point is the sum of the individual velocities. The system mobility. 
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+ve 



x 

Fj vt 



Fig. 3.18. Single degree of freedom system with elements connected in series. 



m = v = K± K±K' 

F F 

- + M c , 

1 s 1 

= b 7 + “• 

ms k c 

Thus when system elements are connected in series their mobilities are added to give 
the system mobility. 

In the system shown in Fig. 3.19, the system comprises a spring and damper 
connected in series with a body connected in parallel. 





Fig. 3.19. Single degree of freedom system and impedance analysis model. 



Thus the spring and damper mobilities can be added or the reciprocal of their 
impedances can be added. Hence the system driving point impedance Z is given by 
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Z = Z + 



_ 1 _ _ 1 _ 

a + z c 



_i_ i 

[_fc/s + cj 



= ms + 

mcs 2 + mks + kc 



cs + k 



Consider the system shown in Fig. 3.20. The spring k j and the body m 1 are connected 
in parallel with each other and are connected in series with the damper c l . 

Thus the driving point impedance Z is 

Z — Z + Zit Z z < 

m 2 *2 c 2 1 




where 



and 



M 1 = M +M 2 , M 2 =^, 
Z 2 — Z ki + Z mj . 



Thus 
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Z — Z + Z k Z H — — 
m 2 *2 c 2 | 



+ 



z r z t + z„ 



Hence 

m 1 m 2 s 4 4- (wjC 2 4- 4- mjcjs 3 + (m^! + m 2 k 2 + C!C 2 )s 2 + 

z = {cik 2 + c 2 k x + + kyk 2 

sim^s 2 + CjS + /q) 

The frequency equation is given when the impedance is made equal to zero or 
when the mobility is infinite. Thus the natural frequencies of the system can be found 
by putting s = jco in the numerator above and setting it equal to zero. 

To summarize, the mobility and impedance of individual elements in a dynamic 
system are calculated on the basis that the velocity is the relative velocity of the two 
ends of a spring or a damper, but the absolute velocity of the body. Individual 
impedances are added for elements or subsystems connected in parallel, and individual 
mobilities are added for elements or subsystems connected in series. 

Example 31 

Find the driving point impedance of the system shown in Fig. 3.6, and hence obtain 
the frequency equation. 

The system of Fig. 3.6 can be redrawn as shown. 



Fe ]vt 
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The driving point impedance is therefore 

1 



Z - + “j“ 



— + 



1 



k l 



= m.s H 1- 



z m2 + z k 

1 



s 1 1 

+ 



k/s m 2 s + k 2 /s 

_ ( m t s 2 + ki)(m 2 s 2 + k + k 2 ) + (m 2 s 2 + k 2 )k 
s(m 2 s 2 + k + k 2 ) 

The frequency equation is obtained by putting Z = 0 and s — jro, thus: 
(k j — m 1 a> 2 )(k + k 2 — m 2 a> 2 ) + k(k 2 — m 2 a> 2 ) = 0. 




4 



The vibration of systems with distributed mass 
and elasticity 



Continuous systems such as beams, rods, cables and strings can be modelled by 
discrete mass and stiffness parameters and analysed as multi degree of freedom systems, 
but such a model is not sufficiently accurate for most purposes. Furthermore, mass 
and elasticity cannot always, be separated in models of real systems. Thus mass and 
elasticity have to be considered as distributed parameters. 

For the analysis of systems with distributed mass and elasticity it is necessary to 
assume an homogeneous, isotropic material which follows Hooke’s law. 

Generally free vibration is the sum of the principal modes. However, in the unlikely 
event of the elastic curve of the body in which motion is excited coinciding exactly 
with one of the principal modes, only that mode will be excited. In most continuous 
systems the rapid damping out of high-frequency modes often leads to the fundamental 
mode predominating. 



4.1 WAVE MOTION 



4.1.1 Transverse vibration of a string 

Consider a vibrating flexible string of mass p per unit length and stretched under a 
tension T as shown in Fig. 4.1. 

The FBDs of an element of length dx of the string assuming small deflections and 
slopes are given in Fig. 4.2. If the lateral deflection y is assumed to be small, the 
change in tension with deflection is negligible, and sin 6 ~ 6. 

The equation for motion in the y direction is 



86 

dx 



^2 
o y 



6 + — dx ) — TO = p dx 
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Fig. 4.2. (a) Applied forces; (b) effective force. 



SO 



' T 8t 2 



Since 



6 = dy/dx. 

where c ~^™ (41) 

This is the wave equation, the solution of which is given in section 4.1.4. The velocity 
of propagation of waves along the string is c. 

4.1.2 Longitudinal vibration of a thin uniform bar 

Consider the longitudinal vibration of a thin uniform bar of cross-sectional area S, 
density p and modulus E under an axial force P, as shown in Fig. 4.3. 
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dx. -^d* 
dx 

Fig. 4.3. Longitudinal beam vibration 



C 2 U 



From Fig. 4.3, 

8P 

d x — p S d.x v . 
cx dr 

Now strain cu/cx — P/SE, 

cP/cx — SE(c 2 u/dx 2 

Thus 



so 



c 2 u/ct 2 = (E!p)(c 2 u/cx 2 ), 



or 



c 2 u/cx 2 = (l/c 2 )(c 2 u/ct 2 ), where c = yJ(E/p). (4.2) 

This is the wave equation. The velocity of propagation of the displacement or stress 
wave in the bar is c. 

4.1.3 Torsional vibration of a uniform shaft 

In the uniform shaft shown in Fig. 4.4 the modulus of rigidity in torsion is G, the 
material density is p, and the polar second moment of area about the axis of tw’ist is J. 
From Fig. 4.4, 

cT J c 2 0 

— d.x = p J d.x . 
cx ct 



The elastic equation is 

c6 
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or 




This is the wave equation. The velocity of propagation of the shear stress wave in 
the shaft is c. 



4.1.4 Solution of the wave equation 

The wave equation 

ex 2 \c 2 J\et 2 ) 

can be solved by the method of separation of variables and assuming a solution of 
the form 

y(x, t) = F(x)G(t). 



Substituting this solution into the wave equation gives 



d 2 F(x) 

dx 2 



G(r) = 



i e 2 G{t) 



F(x), 



that is, 

1 8 2 F(x) _ 1 1 8 2 G{t) 

F(x) dx 2 c 2 G(r) dt 2 

The LHS is a function of x only, and the RHS is a function of t only, so partial 
derivatives are no longer required. Each side must be a constant, — (co/c) 2 say. (This 
quantity is chosen for convenience of solution.) Then 
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d 2 F(x) 
dx 2 




F(x) = 0, 



and 



d 2 G(r) 



+ co 2 G(t) = 0. 



Hence 



CO 



F(x) — A sin^— ) x + B cos( — 



and 



G{t) = C sin cot 4- D cos c ot. 

The constants A and B depend on the boundary conditions, and C and D on the 
initial conditions. The complete solution to the wave equation is therefore 



u = (yA x + B cos^jxj^C sin wt + D cos 



(4.4) 



Example 32 

Find the natural frequencies and mode shapes of longitudinal vibrations for a 
free — free beam with initial displacement zero. 

Since the beam has free ends, cu/cx = 0 at x = 0 and x = /. Now 



~ = { a( — \ cos( — ] x — B ( — ) sin( — )x ) ( C sin cot + D cos cot 

cx 



Hence 



and 



— ] = A\ — 1 (C sin cot + D cos cot) = 0, so that A =. 0, 

cx)x= o \ c 



CO 



col 



= ( — ]( —5 sin — I C sin cot + D cos cot | = 0. 



Thus sin(co//c) = 0, since B ^ 0, and therefore 
col col 



C yJiE/p) 



— n, 2n, ■ ■ ■ nn ■ ■ ■ 



That is 



nn // E 



I V \p 



rad/s, 
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where a> = c/wavelength. These are the natural frequencies. 
If the initial displacement is zero, D = 0 and 




where B' — B x C. Hence the mode shape is determined. 

Example 33 

A uniform vertical rod of length l and cross-section S is fixed at the upper end and 
is loaded with a body of mass M on the other. Show that the natural frequencies of 
longitudinal vibration are determined by 

a)lyJ{p/E) tan a>lyJ(p/E) = Spl/M 




A F 



M 



At 

x = 0, u = 0, and at x = /, F — SE (du/dx). 

Also 

F = SE (du/dx) = — M(d z u/dt 2 ). 

The general solution is 

u = (A sin(co/c)x + B cos(oj/c)x)(C sin cot + D cos cot). 
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Now, 

« x =o = 0, so B — 0, 

thus 

u = (A sin (cu/c)xXC sin oot + D cos cut), 

( du/dx) x =i = (A(ro/c) cos(cu//c)XC sin cot + D cos art), 

and 

( d 2 u/dt 2 ) x==i = ( — Aco 2 sin(col/c)XC sin art + D cos cot), 
so 

F = SEA (co/c) cos(o>l/c)(C sin cot + D cos art) 

= MAoo 2 sin(a>//c){C sin cot + D cos cot). 

Hence ( col/c ) tan(oJ//c) = SIE/Mc 2 , and 

a)ly/{p/E) tan colyJip/E) = Spl/M, since c 2 - Ejp. 

4.2 TRANSVERSE VIBRATION 

4.2.1 Transverse vibration of a uniform beam 

The transverse or lateral vibration of a thin uniform beam is another vibration 
problem in which both elasticity and mass are distributed. Consider the moments 
and forces acting on the element of the beam shown in Fig. 4.5. The beam has a 
cross-sectional area A, flexural rigidity El, and material of density p. 




For element Q 

4 



M 



( 






dx 



Q*T~ dx 
dx 



Applied force 



♦ a 

pAdx — \ 
it 2 

Effective force 



Fig. 4.5. Transverse beam vibration. 
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Then for the element, neglecting rotary inertia and shear of the element, taking 
moments about O gives 

dx dx dQ , dx dM , 

M + e T + e T + ^d* T -M + — <■*. 



That is, 

Q = dM/dx. 

Summing forces in the y direction gives 

8 2 y 



(4.5) 



8Q 

dx 



dx — pA dx — 5 - . 

dt 2 



Hence 



d 2 M 

dx 2 



= pA 



8 2 y 

dt 2 ' 



Now El is a constant for a prismatical beam, so 



M = - El 



d 2 y 

dx 2 



and 



d 2 M d 4 y 

dx 2 dx 4 



Thus 



d 4 y ( py4\ d 2 y 
dx 4 + \El)dt 2 



0 . 



(4.6) 



This is the general equation for the transverse vibration of a uniform beam. 

When a beam performs a normal mode of vibration the deflection at any point of 
the beam varies harmonically with time, and can be written 



y = X (B l sin cot + B 2 cos a>t), 

where X is a function of x which defines the beam shape of the normal mode of 
vibration. Hence 



d*X 

dx 4 




cu 2 X = a 4 X, 



(4.7) 



where 

A 4 — pAco 2 /EI. This is the beam equation. 

The general solution to the beam equation (4.7) is 

X = C 1 cos Ax + C 2 sin Ax + C 3 cosh Ax + C 4 sinh Ax, 

where the constants C 1 2 3A are determined from the boundary conditions. 

For example, consider the transverse vibration of a thin prismatical beam of length l, 
simply supported at each end. The deflection and bending moment are therefore zero 
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at each end, so that the boundary conditions are X = 0 and d 2 X/dx 2 = 0 at x = 0 
and x — l. 

Substituting these boundary conditions into the general solution above gives at 



X = 0, 


* 

II 

o 


thus 0 = Cj + C 3 , 


x = 0, 


d 2 X n 
dx 2 ~ ° : 


; thus 0 = Cj — C 3 



That is, 

C 1 = C 3 — 0 and X = C 2 sin kx + C 4 sinh kx. 

Now at 

x = /, X = 0 so that 0 — C 2 sin ?J + C 4 sinh kl, 

and at 

d 2 X 

x = /, — — r- = 0, so that 0 = C 2 sin kl — C 4 sinh kl. 

dx 

That is, 

C 2 sin kl = C 4 sinh kl = 0. 

Since 



kl ^ 0, sinh kl j= 0 and therefore C 4 = 0. 

Also C 2 sin kl = 0. Since C 2 / 0, otherwise X - 0 for all x, then sin kl = 0. Hence 
X = C 2 sin kx, and the solutions to sin kl — 0 give the natural frequencies. These are 



so that 






. . rad/s. 



/. = 0, co = 0 is a trivial solution because the beam is at rest. The lowest or first 
natural frequency is therefore aq = (n/l) 2 J(EI/Ap) rad/s, and the corresponding mode 
shape is X = C 2 sin nx/l. This is the first mode. <x> 2 = ( 2n/l) 2 y/{EI/Ap ) rad/s is the 
second natural frequency, and the second mode is X = C 2 sin 2nx/l, and so on. The 
mode shapes are drawn in Fig. 4.6. 

These sinusoidal vibrations can be superimposed so that any initial conditions can 
be represented. Other end conditions give frequency equations with the solution 





rad/s. 



where the values of x are given in Table 4.1. 
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1st mode shape, one half-wave: 



y = C 2 sin n ( - {(l?! sin wj + B 2 cos cojt); 



"i 




rad/s. 






Fig. 4.6. Transverse beam vibration mode shapes and frequencies. 



Sec. 4.2] 



Transverse vibration 151 



Table 4.1 



End conditions 


Frequency equation 


1st mode 


2nd mode 


3rd mode 


4th mode 


5th mode 


Clamped-free 


cos A/ cosh /./ = — 1 


3.52 


22.4 


61.7 


121.0 


199.9 


Pinned-pinned 

Clamped- 


sin Af = 0 


9.87 


39.5 


88.9 


157.9 


246.8 


pinned 
Clamped- 
clamped or 


tan A/ = tan /.l 


15.4 


50.0 


104.0 


178.3 


272.0 


Free-free 


cos A/ cosh A/ = 1 


22.4 


61.7 


121.0 


199.9 


298.6 



The natural frequencies and mode shapes of a wide range of beams and structures 
are given in Formulas for Natural Frequency and Mode Shape by R. D. Blevins (Van 
Nostrand, 1979). 

4.2.2 The whirling of shafts 

An important application of the theory for transverse beam vibration is to the whirling 
of shafts. If the speed of rotation of a shaft is increased, certain speeds will be reached 
at which violent instability occurs. These are the critical speeds of whirling. Since the 
loading on the shaft is due to centrifugal effects the equation of motion is exactly the 
same as for transverse beam vibration. The centrifugal effects occur because it is 
impossible to make the centre of mass of any section coincide exactly with the axis 
of rotation, because of a lack of homogeneity in the material and other practical 
difficulties. 

Example 34 

A uniform steel shaft which is carried in long bearings at each end has an effective 
unsupported length of 3 m. Calculate the first two whirling speeds. 

Take I/A = 0.1 x 10' 3 m 2 , E = 200 GN/m 2 , and p - 8000 kg/m 3 . 

Since the shaft is supported in long bearings, it can be considered to be ‘built in’ 
at each end so that, from Table 4.1, 

£7 
Ap 

where atj = 22.4 and x 2 = 61.7. For the shaft, 

M)-A 

so that the first two whirling speeds are: 

22.4 

oj 1 = 50 = 124.4 rad/s, 



200 x 10 9 x 0.1 x 10' 
8000 



= 50 m 2 /s, 





so 
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to 124.4 

/, = — - = - — = 19.8 cycle/s, and N t = 1188 rev/min, 
2 n 2n 



and 



N, = 1188 = 3272 rev/min. 

‘ 22.4 

Rotating this shaft at speeds at or near to the above will excite severe resonance 
vibration. 



4.2.3 Rotary inertia and shear effects 

When a beam is subjected to lateral vibration so that the depth of the beam is a 
significant proportion of the distance between two adjacent nodes, rotary inertia of 
beam elements and transverse shear deformation arising from the severe contortions 
of the beam during vibration make significant contributions to the lateral deflection. 
Therefore rotary inertia and shear effects must be taken into account in the analysis 
of high-frequency vibration of all beams, and in all analyses of deep beams. 

The moment equation can be modified to take into account rotary inertia by 
adding a term pi d 3 y/(cx dt 2 ), so that the beam equation becomes 



d 4 y 

El —^7 — pi 
dx 4 P 



° z y , A d *y -n 

T — -~r~; + pA -r-j — 0. 

dx dt 2 dt 2 



Shear deformation effects can be included by adding a term 

EIp d 4 ) ! 
kg dx 2 ct 2 ' 



where k is a constant whose value depends on the cross-section of the beam. Generally, 
k is about 0.85. The beam equation then becomes 



El 



O 



dx 4 



EIp 
~kg dx 2 dt 2 



^ab- 



solutions to these equations are available, which generally lead to a frequency a 
few percent more accurate than the solution to the simple beam equation. However, 
in most cases the modelling errors exceed this. In general, the correction due to shear 
is larger than the correction due to rotary inertia. 



4.2.4 The effect of axial loading 

Beams are often subjected to an axial load, and this can have a significant effect on 
the lateral vibration of the beam. If an axial tension T exists, which is assumed to be 
constant for small amplitude beam vibrations, the moment equation can be modified 
by including a term T d 2 y/cx 2 , so that the beam equation becomes 
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Tension in a beam will increase its stiffness and therefore increase its natural 
frequencies; compression will reduce these quantities. 



Example 35 

Find the first three natural frequencies of a steel bar 3 cm in diameter, which is simply 
supported at each end, and has a length of 1.5 m. Take p = 7780 kg/m 3 and 
E = 208 GN/m 2 . 



For the bar, 

lfEI\ /208 x 10 9 7 c( 0.03) 4 /64 /2 , c 0 /2 

\l\ApJ V rc{0.03/2) 2 7780 ' ' 

Thus 

a>! = 38.8 = 170.2 rad/s and /, — 27.1 FIz. 



Hence 

f 2 = 27.1 x 4 = 108.4 Hz, 



and 



/ 3 = 27.1 x 9 = 243.8 Hz. 

If the beam is subjected to an axial tension T, the modified equation of motion 
leads to the following expression for the natural frequencies: 

aitA 2 T f wr\ 4 El 
l ) Ap + \ l ) Ap' 

For the case when T = 1000 the correction to aj 2 is oi 2 c , where 
a ( 1000 \ 

U.03/2)^ 7780 / ~ 795 (rad/S > 2 
4.5 Hz. Hence /j - V(4.5 2 + 27. 1 2 ) = 27.5 Hz. 



iol - 



Oil — 



It 

L5 



That is,/ c 



4.2.5 Transverse vibration of a beam with discrete bodies 

In those cases where it is required to find the lowest frequency of transverse vibration 
of a beam which carries discrete bodies, Dunkerley’s method may be used. This is a 
simple analytical technique which enables a wide range of vibration problems to be 
solved using a hand calculator. Dunkerley’s method uses the following equation: 



1 



co 2 , P, 



1 1 
2 F 



1 



P 2 + Pi 2 + 



1 

P 2 

r 4 



+ .. 



where tOj is the lowest natural frequency of a system and P 1 , P 2 , F 3 ,... are the 
frequencies of each body acting alone. 

This equation may be obtained for a two degree freedom system by writing the 
equations of motion in terms of the influence coefficients as follows: 
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y i = a n mi<y 2 yi + a l2 m 2 co 2 y 2 , 

and 

y 2 = + a 2 2 m 2" 2 >'2- 

The frequency equation is given by 
a ll m l a) 2 — 1 a 12 m 2 co 2 
a 2i m i a)2 a 2 2 m 2 ft>2 — 1 

By expanding this determinant, and solving the resulting quadratic equation, it is 
found that 



"l,2 2 



Ml + «22 ffl 2) ± >/[(«ll«l + «22 m l) 2 - 4 C«11«22 ~ «21 g 12)] 

2(^11*22 — a 2i a i2) 



Hence it can be shown that 



1 1 

2 ^ 2 = a H m i + a 22 W 2- 

(Uj co 2 

Now Pj is the natural frequency of body 1 acting alone, hence 



1 m i aum,' 
Similarly 



2 — ~ m ‘ 
a 22 W 2 

Thus 



1 1 _ 1 1 

ft)! 2 + 0) 2 2 ~ P," 1 + 

A similar relationship can be derived for systems with more than two degrees of 
freedom. 

If ft> 2 the left-hand side is approximately 1/co 2 , hence 
1 1 1 

W^p? + K 2 ' 

Example 36 

A steel shaft (p = 8000 kg/m 3 , E = 210 GN/m 2 ) 0.055 m diameter, running in self- 
aligning bearings 1.25 m apart, carries a rotor of mass 70 kg, 0.4 m from one bearing. 
Estimate the lowest critical speed. 

For the shaft alone 

Pj = (n/l) 2 J(EI/Ap) rad/s = 141 rad/s = 1350 rev/min. 
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This is the lowest critical speed for the shaft without the rotor. For the rotor alone, 
neglecting the mass of the shaft, 

P 2 = yj{ k/m ) rad/s 

k = 3 £///(x 2 (/ - x) 2 ) 

where x = 0.4 m and / = 1.25 m. 

Thus 

k = 3.07 MN/m 

and 

P 2 - 209 rad/s = 2000 rev/m in. 

Now 

1/N t 2 = 1/1350 2 + 1/2000 2 , hence N t = 1135 rev/min 





L 


1 


/ 



4.26 Receptance analysis 

Many dynamic systems can be considered to consist of a number of beams fastened 
together. Thus if the receptances of each beam are known, the frequency equation 
of the system can easily be found by carrying out a subsystem analysis (Section 3.2.3). 
The required receptances can be found by inserting the appropriate boundary 
conditions in the general solution to the beam equation. 

It will be appreciated that this method of analysis is ideal for computer solutions 
because of its repetitive nature. 

For example, consider a beam which is pinned at one end (x = 0) and free at the 
other end (x = /). This type of beam is not commonly used in practice, but it is useful 
for analysis purposes. With an harmonic moment of amplitude M applied to the 
pinned end, 



at x = 0, X = 0 (zero deflection) and 



d 2 X _ M 
dx 2 ~ £7 



(bending moment M) 



and at x = /, 



d 2 X 

dx 2 



0 



and 




(zero bending moment) 



(zero shear force). 



Now, in general, 

X = Cj cos /.x + C 2 sin /.x + C 3 cosh Ax + C 4 sinh ax. 



Thus applying these boundary conditions, 
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0 = c, + c 3 



and 



M 

Tl 



-C x X 2 + C 3 X 2 . 



Also 

0 = — C x X 2 cos XI — C 2 X 2 sin XI + C 3 X 2 cosh A/ + C 4 A 2 sinh A/. 

and 

0 = C t A 3 sin A/ — C 2 A 3 cos XI + C 3 A 3 sinh XI + C 4 A 3 cosh XL 

By solving these four equations C x 2 3 4 can be found and substituted into the 
general solution. It is found that the receptance moments/slope at the pinned end is 

(1 + cos XI cosh AO 
EIX (cos XI sinh XI — sin XI cosh XI) 

and at the free end is 

2 cos XI cosh XI 

EIX (cos XI sinh XI — sin XI cosh XI) 

The frequency equation is given by 

cos XI sinh X.l — sin A/ cosh XI = 0. 

That is, tan XI = tanh XL 

Moment/deflection receptances can also be found. 

By inserting the appropriate boundary conditions into the general solution, the 
receptance due to an harmonic moment applied at the free end, and harmonic forces 
applied to either end, can be deduced. Receptances for beams with all end conditions 
are tabulated in The Mechanics of Vibration by R. E. D. Bishop & D. C. Johnson 
(CUP, 1960/79), thereby greatly increasing the ease of applying this technique. 

Example 37 

A hinged beam structure is modelled by the array shown below: 



Beam Hinge 




The hinges are pivots with torsional stiffness k T and their mass is negligible. All 
hinges and beams are the same. 

It is required to find the natural frequencies of free vibration of the array, so that 
the excitation of these frequencies, and therefore resonance, can be avoided. 
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Since all the beams are identical, the receptance technique is relevant for finding 
the frequency equation. This is because the receptances of each subsystem are the 
same, which leads to some simplification in the analysis. 

There are two approaches: 

(i) to split the array into subsystems comprising torsional springs and beams, or 

(ii) to split the array into subsystems comprising spring-beam assemblies. 

This approach results in a smaller number of subsystems. 

Considering the first approach, and only the first element of the array, the 
subsystems could be either 

(a) 




or 




For (a) the frequency equation is a n + Pu = 0, whereas for (b) the frequency 
equation is 

a n + Pll a 12 + P 12 _ q 
a 21 + Pi 1 a 22 + P 22 

where a n is the moment/slope receptance for A, /J u is the moment/slope receptance 
for B, /? 12 is the moment/deflection receptance for B, f$ 22 is the force/deflection 
receptance for B, and so on. 

For (a), either calculating the beam receptances as above, or obtaining them from 
tables, the frequency equation is 

1 cos /.l cosh // + 1 

k T El /.(cos /J sinh /.l — sin /./ cosh /.l) 

where 




For (b), the frequency equation is 
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cos /.l sinh /./ + sin /.l cosh a/ 
EI /.( cos /./ cosh /./ — 1) 



— sin /.l sinh // 
£/a 2 ( cos l\ cosh a/ — 1) 



I —sin a/ sinh /./ 

' EI /. 2 ( cos a/ cosh a/ — 1) 

which reduces to the equation given by method (a). 

The frequency equation has to be solved after inserting the structural parameters, 
to yield the natural frequencies of the structure. 

For the whole array it is preferable to use approach (ii), because this results in a 
smaller number of subsystems than (i), with a consequent simplification of the 
frequency equation. However, it will be necessary to calculate the receptances of the 
spring pinned — free beam if approach (ii) is adopted. 



— (cos /./ sinh /J — sin a/ cosh a/) 
£/a 3 ( cos a/ cosh a/ — 1) 



The analysis of structures such as frameworks can also be accomplished by the 
receptance technique, by dividing the framework to be analysed into beam substruc- 
tures. For example, if the in-plane natural frequencies of a portal frame are required, 
it can be divided into three substructures coupled by the conditions of compatibility 



ABC 







Deflection, 2 Deflection, 4 

Fig. 4.7. Portal frame substructure analysis. 



and equilibrium, as shown in Fig. 4.7. 

Substructures A and C are cantilever beams undergoing transverse vibration, 
whereas B is a free-free beam undergoing transverse vibration. Beam B is assumed 
rigid in the horizontal direction, and the longitudinal deflection of beams A and C 
is assumed to be negligible. 

Because the horizontal member B has no coupling between its horizontal and 
flexural motion ft 12 = /? 14 = /? 23 = /? 34 = 0, so that the frequency equation becomes 
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^11 + Pll a i2 Pl3 

a 21 a 2 2 + P 22 0 

P 3 I 0 /’ 3 3 + Pl3 

0 P 42 743 

4.3 THE ANALYSIS OF CONTINUOUS SYSTEMS BY RAYLEIGH’S 
ENERGY METHOD 

Energy methods can be used in the analysis of continuous systems if the mode shape 
of the system is known; this is necessary for evaluating the kinetic and potential 
energies. It has been shown by Rayleigh that the lowest natural frequency of such 
systems can be fairly accurately found by assuming any reasonable deflection curve 
for the shape of the first mode. Naturally the curve must be compatible with the end 
conditions of the system, and since any deviation from the true mode shape puts 
additional constraints on the system, the frequency determined by Rayleigh’s method 
is never less than the exact frequency. Generally the static deflection curve of the 
system gives a frequency within a few per cent of the exact value. The frequency of 
vibration is found by considering the conservation of energy in the system; the natural 
frequency is determined by dividing the expression for potential energy in the system 
by the expression for kinetic energy. 

4.3.1 The vibration of systems with heavy springs 

The mass of the spring elements can have a considerable effect on the frequency of 
vibration of those systems in which relatively heavy springs are used, such as in valve 
mechanisms and in some rotational systems. 

Consider the translational system shown in Fig. 4.8, where a rigid body of mass 
M is connected to a fixed frame by a spring of mass m, length /, and stiffness k. The 
body moves in the x direction only. If the dynamic deflected shape of the spring is 
assumed to be the same as the static shape under an applied axial force, the velocity 
of the spring element is y — (y//)x, and the mass of the element is (m/l)dy. 



0 

P 24 

P 34 

7 44 P 44 



= 0 . 




Fig. 4.8. Single degree of freedom system with heavy spring. 
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Thus 



1 



Mx 2 + 



7 1 (m 



2 V / 






d y 



= - M + — x 



and 



V = \kx 2 . 



Assuming simple harmonic motion and putting T mM = F max gives the frequency of 
free vibration as 



/ = 



1 



2n \M + (m/3) 



Hz. 



That is, if the system is to be modelled with a massless spring, one third of the actual 
spring mass must be added to the mass of the body in the frequency calculation, 
d 

Alternatively, — (T + V) = 0 can be used for finding the frequency of oscillation, 
dt 



4.3.2 Transverse vibration of a beam 

For the beam shown in Fig. 4.9, if m is the mass/length and y is the amplitude of the 
assumed deflection curve, then 



i 



dm 



* 

= W 

* 



y 2 dm, 



where a> is the natural circular frequency of the beam. 

The strain energy of the beam is the work done on the beam which is stored as 
elastic energy. If the bending moment is M and the slope of the elastic curve is 9, 

F=^j*Md0. 



Usually the deflection of beams is small so that the following relationships can be 
assumed to hold: 



thus 



dy 

9 = — - - and 
dx 



1 d0 d 2 y 

R dx dx 2 ’ 



R 69 = dx, 



From beam theory, M/I — E/R, where R is the radius of curvature and El is the 
flexural rigidity. Thus 



F = i 



2 J R 



6x 



El 



d 2 y 

dx' 



dx. 
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Beam 



Fig. 4.9. Beam deflection. 



Since 



T — V 

max max’ 



j £, lS 



dx 



J y 2 dm 



This expression gives the lowest natural frequency of transverse vibration of a beam. 
It can be seen that to analyse the transverse vibration of a particular beam by this 
method requires y to be known as a function of x. For this the static deflected shape 
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or a part sinusoid can be assumed, provided the shape is compatible with the beam 
boundary conditions. 

Example 38 

A simply supported beam of length l and mass m 2 carries a body of mass mj at its 
mid-point. Find the lowest natural frequency of transverse vibration. 




The boundary conditions are y = 0 and d 2 y/dx 2 = 0 at x = 0 and x = /. These 
conditions are satisfied by assuming that the shape of the vibrating beam can be 
represented by a half sine wave. A polynomial expression can be derived for the 
deflected shape, but the sinusoid is usually easier to manipulate. 

y — y 0 sin(7tx//) is a convenient expression for the beam shape, which agrees with 
the boundary conditions. Now 



y = j 0 sinl 




Hence 



and 



d 2 y 

dx 2 





Jo 





dx 



y 2 dm — 




sin 2 





+ yl «i 



y o 




Thus 



and 
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2 EIWlfl/2, v2 
(o — __ — — (rad/s) . 



(Wi + mjl) 



If 



m 2 = 0 , 



, El n 4 El 

a) 2 = — — = 48.7 



2 l 3 m 1 ' mj 3 ' 

The exact solution is 48 El/myl 3 , so the Rayleigh method solution is 1.4% high. 
Alternatively the Dunkerley method can be used. Here, 






Thus 



48 El 
m x l 3 

mj 3 



and 



P, = 



El n* 
m 2 l 3 . 



(O 



+ 



m 2 l 3 



Hence 



48 El tcEI 



El T r 



1.015m! + -~ 



(rad/s) 2 



which is very close to the value determined by the Rayleigh method. 



Example 39 

Find the lowest natural frequency of transverse vibration of a cantilever of mass m, 
which has rigid body of mass M attached at its free end. 




The static deflection curve is y = (y 0 /2l 3 )(Vx 2 — x 3 ). Alternatively 
y = y 0 ( 1 — cos nx /2 1) could be assumed. Hence 



'e,(^Ux 

dx' 



El 



y<> 



3EI 



2/ 3 ) (6/-6x) 2 dx- /3 



To. 



and 



y 2 dm = 



y 2 jdx + y 0 2 M, 
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r ‘ y 0 2 m 



4/ 6 / 



— (3/x 2 — .x 3 ) 2 dx + j^ 0 M, 






Thus 




3 El 






(rad/s) 2 . 



Example 40 

A pin-ended strut of length / has a vertical axial load P applied. Determine the 
frequency of free transverse vibration of the strut, and the maximum value of P for 
stability. The strut has a mass m and a second moment of area /, and is made from 
material with modulus of elasticity E. 



P 

o 




The deflected shape can be expressed by 
y = >’o sin n 

since this function satisfies the boundary conditions of zero deflection and bending 
moment at x = 0 and x = /. Now 
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1 

2 




dx — Pz, 



where 




4 l ' 

Thus 

v 2 
ma ‘ V 4 / 3 4 / / 0 



Now, 



Thus 



T =- 

*max 2 



y 2 dm = - 



m 



y y dx. 



0 

If' 2-2 x m J m 2 

>o sin 71 7 T dx = T • 
(0 / f 4 



El n* P n 2 

4 ~ 1* ~ 4 T 

m 

4 



(rad/s) 2 . 



and 





Hz. 
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For stability / > 0, hence EI(n/l) 2 > P. That is, F mai = EI(n/l) 2 . EIn 2 /! 2 is known as 
the Euler buckling load. 

Example 41 

Part of an industrial plant incorporates a horizontal length of uniform pipe, which 
is rigidly embedded at one end and is effectively free at the other. Considering the 
pipe as a cantilever, derive an expression for the frequency of the first mode of 
transverse vibration using Rayleigh’s method. 

Calculate this frequency, given the following data for the pipe: 



Modulus elasticity 200 GN/m 2 

Second moment of area about bending axis 0.02 m 4 

Mass 6 x 10 4 kg 

Length 30m 

Outside diameter 1 m 




For a cantilever, assume 



y - y. (' - cos 

This is compatible with zero deflection and slope when x = 0, and zero shear force 
and bending moment when x — l. Thus 




and 
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y 2 dm = 



nx \ m 



y,[ \ - cos — \ jdx 



= yi . 



Hence, assuming the structure to be conservative, that is the total energy remains 
constant throughout the vibration cycle, 



co 



2 




/ 

2 



yi 

,2 n 




Thus 



co 



3.66 




In this case 



and 



El 200 x 10 9 x 0.02 / 2 
ml 2 ~ 6 x 10 4 x 30 3 / S 



3.66 

2n 




Hz. 



Hence 



co — 5.75 rad/s and f = 0.92 Hz. 

4.3.3 Wind or current excited oscillation 

A structure exposed to a fluid stream is subjected to an harmonically varying force 
in a direction perpendicular to the stream. This is because of eddy or vortex shedding 
on alternate sides of the structure on the leeward side. Tall structures such as masts, 
bridges, and chimneys are susceptible to excitation from steady winds blowing across 
them. The Strouhal number relates the excitation frequency / s , to the velocity of fluid 
flow t’(m/s) and the hydraulic mean diameter D (m) of the structure as follows: 

Strouhal number =^— . 

v 

If the frequency f s is close to the natural frequency of the structure, resonance may 
occur. 

For a structure, 

4 x area of cross-section 
circumference 

so that for a chimney of circular cross-section and diameter d. 
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r , 4W4) </ ! , 

D = -^- =d ’ 

and for a building of rectangular cross-section a x b, 

Aab lab 
2 (a + b) (a + b)' 

Experimental evidence suggests a value of 0.2 — 0.24 for the Strouhal number for 
most flow rates and wind speeds encountered. This value is valid for Reynolds numbers 
in the range 3 x 10 5 — 3.5 x 10 6 . 

Example 42 

For constructing a tanker terminal in a river estuary a number of cylindrical concrete 
piles were sunk into the river bed and left free standing. Each pile was 1 m diameter 
and protruded 20 m out of the river bed. The density of the concrete was 2400 kg/m 3 
and the modulus of elasticity 14 x 10 6 kN/m 2 . Estimate the velocity of the water 
flowing past a pile which will cause it to vibrate transversely to the direction of the 
current, assuming a pile to be a cantilever and taking a value for the Strouhal number 



— = 0.22, 

v 



where f s is the frequency of flexural vibrations of a pile, D is the diameter, and v is 
the velocity of the current. 

Consider the pile to be a cantilever of mass m, diameter D, and length /; then the 
deflection y at a distance x from the root can be taken to be, y — >',(1 — cos nx/2l), 
where y, is the deflection at the free end. 
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Thus 



d 2 y / 7t V / 7tx 

d^ = y <2/J COS U 

n 



El 



o 



d y 
dx 2 



dx = £/ y, - cos — dx 



21 



= \2, 2 



nx 



21 



\ y2 im ~ I 



( nx 



* 1 1 - cos l 27 ) ) l T ) djc 



m 



m /3 4 



-"7 5-;' 



Hence 



£/y, 



co 



2 _ 



2 / 






m 



l \2 



l 



Substituting numerical values gives co = 5.53 rad/s, that is, / = 0.88 Hz. When 
f s = 0.88 Hz resonance occurs; that is, when 

f s D 0.88 ^ , 

C “ 0Tl = M2 = 4 m/S - 



4.4 THE STABILITY OF VIBRATING SYSTEMS 



If a system is to vibrate about an equilibrium position, it must be stable about that 
position. That is, if the system is disturbed when in an equilibrium position, the elastic 
forces must be such that the system vibrates about the equilibrium position. Thus 
the expression for co 2 must be positive if a real value of the frequency of vibration 
about the equilibrium position is to exist, and hence the potential energy of a stable 
system must also be positive. 

The principle of minimum potential energy can be used to test the stability of 
systems which are conservative. Thus a system will be stable at an equilibrium position 
if the potential energy of the system is a minimum at that position. This requires that 



dV 

dq 



= 0 



and 



d 2 T 

dq 2 



^~2>0, 



where q is an independent or generalized coordinate. Hence the necessary conditions 
for vibration to take place are found, and the position about which the vibration 
occurs is determined. 
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4.5 THE FINITE ELEMENT METHOD 

Many structures, such as a ship hull or engine crankcase, are too complicated to be 
analysed by classical techniques, so that an approximate method has to be used. It 
can be seen from the receptance analysis of complicated structures that breaking a 
dynamic structure down into a large number of substructures is a useful analytical 
technique, provided that sufficient computational facilities are available to solve the 
resulting equations. The finite element method of analysis extends this method to the 
consideration of continuous structures as a number of elements, connected to each 
other by conditions of compatibility and equilibrium. Complicated structures can 
thus be modelled as the aggregate of simpler structures. 

The principal advantage of the finite element method is its generality; it can be 
used to calculate the natural frequencies and mode shapes of any linear elastic system. 
However, it is a numerical technique which requires a fairly large computer, and care 
has to be taken over the sensitivity of the computer output to small changes in input. 

For beam type systems the finite element method is similar to the lumped mass 
method, because the system is considered to be a number of rigid mass elements of 
finite size connected by massless springs. The infinite number of degrees of freedom 
associated with a continuous system can thereby be reduced to a finite number of 
degrees of freedom, which can be examined individually. 

The finite element method therefore consists of dividing the dynamic system into 
a series of elements by imaginary lines, and connecting the elements only at the 
intersections of these lines. These intersections are called nodes. It is unfortunate that 
the word node has been widely accepted for these intersections; this meaning should 
not be confused with the zero vibration regions referred to in vibration analysis. The 
stresses and strains in each element are then defined in terms of the displacements 
and forces at the nodes, and the mass of the elements is lumped at the nodes. A series 
of equations are thus produced for the displacement of the nodes and hence the 
system. By solving these equations the stresses, strains, natural frequencies, and mode 
shapes of the system can be determined. The accuracy of the finite element method 
is greatest in the lower modes, and increases as the number of elements in the model 
increases. The finite element method of analysis is considered in Techniques of Finite 
Elements by B. Irons & S. Ahmad (Ellis Horwood, 1980). 




Automatic control systems 



Automatic control systems, wherein a variable quantity is made to conform to a 
predetermined level, have been in use for several centuries; but, as with much 
technological development, the most rapid advances have taken place during recent 
years. Feedforward, or open loop systems, merely control the input, such as in the 
case of a machine tool cutter which simply follows a given guide or pattern. An early 
example of an open loop system is the Jacquard loom of 1801, in which a set of 
punched cards programmed the patterns to be woven by the loom, and no information 
from the process or results was used to correct the loom operation. Feedback, or 
closed loop, systems feed back information from the process to control the operation 
of the machine. One of the earliest closed loop systems was that used by the Romans 
to maintain water levels in their aqueducts by means of floating valves. Later, windmills 
were the spawning ground of several control systems, for example the sails were 
automatically kept into the wind by means of a fantail (1745), centrifugal governors 
were used to control the speed of the millstones (1783), and the speed of rotation of 
the sails was automatically controlled by roller reefing (1789). In the late eighteenth 
century centrifugal governors were also being used to control the speed of steam 
engines by regulating the steam supply. These devices provided much closer control 
than manual systems could, and they were cheaper to operate, so that the overall 
efficiency of the machine increased. This led to demands for even better control 
systems which operated within narrower margins so that the efficiency was further 
increased. 

In principle, many variables can be controlled by humans, but in practice this may 
be impossible, difficult, costly, or undesirable because of the need for continuous 
operation regardless of environment, large forces, and a fast response. The human 
reaction time of about 0.3 seconds is too slow for many applications. Further examples 
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of control systems are liquid level control by ball valve, temperature control by 
thermostat, and surface control such as a ship’s rudder or aircraft flaps by hydraulic 
servo. However, all types of control system can be modelled for analysis purposes, 
irrespective of the operating mechanism. The first theoretical analysis of a control 
system was published by Maxwell in the nineteenth century: this theory was soon 
generalized and followed by a large number of contributions to explain and improve 
the operation of control systems. In the early twentieth century the rapid development 
of automatic control systems took place because of the need to position guns and 
ships quickly and accurately, which led to the development of servomechanisms in 
the mid-1930s. In the 1950s the potential of multiple loop systems was investigated, 
and the introduction of computers opened the way for much greater complexity in 
control systems. Computer control is usually applied to industrial problems in one 
of three ways: supervisory control which continually adjusts the plant to optimum 
operation conditions, direct digital control, and hierarchy control which integrates 
the plants operation at every level from management decisions through to valve 
settings. For example, the new Waterloo rail terminal has a building management 
system which monitors and controls air conditioning, lighting, escalators and 
travelaters. 

Computer control is aided by the ability to measure and convert into electrical 
signals a wide range of system parameters such as temperature, pressure, speed, level, 
weight, flow, conductivity, and thickness. A fast rate of progress has been maintained 
to the present day, particularly since the introduction of microprocessors, so that 
only an introduction to automatic control systems can be attempted in a text of this 
length. 

The essential feature of an automatic control system is the existence of a feedback 
loop to give good performance. This is a closed loop system; if the measured output 
is not compared with the input the loop is open. Usually it is required to apply a 
specific input to a system and for some other part of the system to respond in the 
desired way. The error between the actual response and the ideal response is detected 
and fed back to the input to modify it so that the error is reduced, as shown in Fig. 
5.1. The output of a device represented by a block in a block diagram cannot affect 
the input to that device unless a specific feedback loop is provided. 




Fig. 5.1. Control system block diagram. 
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In general both input and output vary with time, and the control system can be 
mechanical, pneumatic, hydraulic and electrical in operation, or any combination of 
these or other power sources. The system should be absolutely stable so that if excited 
it will settle to some steady value, and it should be accurate in the steady state. 

This concept can be illustrated by considering a simple dynamic system of the type 
previously considered, as shown in Fig. 5.2. 



OUTPUT 
1 ► *o ,{| 




vVW 




Fig. 5.2. Dynamic system. 



INPUT 

► jr, m 



The system comprises a body of mass m connected to an input controlled platform 
by a spring and viscous damper. A specific input x s (r) is applied to this platform, and, 
as a control system, the response of the body or output x 0 {t) should be identical to 
the input. 

Considering the forces acting on the body, the equation of motion is 

mx 0 = k(x j - x 0 ) + c(x i - xj. (5.1) 

Equations of motion of this type have been solved for an harmonic input (equation 
2.18). For a general solution irrespective of input it is convenient to use the D-operator. 
Thus (5.1) becomes 

mD 2 x c = k(x { - x 0 ) + cD(xj - x D ) = (k + cD)(x; - x 0 ). 

It should be noted that although using the D-operator is a neat and compact form 
of writing the equation it does not help with the solution of the response problem. 
Now the force F on the body is mD 2 x 0 , so F = mD 2 x 0 , or 




The transfer function of a system is the function by which the input is multiplied to 
give the output, so that since F is the input to the body and x Q the output, (1/mD 2 ) 
is the transfer function (TF) of the body. 
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This is shown in block diagram form in Fig. 5.3. 



F 


1 


f ' \ 


*0 




Ud 2 / 





Fig. 5.3. Rigid body block diagram. 



For the spring/damper unit F = (k + cDXXj — x 0 ) as shown in Fig. 5.4. 

F 



x - x„ 



(k + cD) 



Fig. 5.4. Spring/damper block diagram. 



Because the input to the spring/damper unit is (x ; — x Q ) and the output is F, the TF 
is ( k + cD). These systems can be combined as shown in Fig 5.5. 




Fig. 5.5. Block diagram. 




Fig. 5.6. System block diagram. 

Fig. 5.6 is the conventional unity feedback loop form. Essentially, the spring/damper 
acts as an error-sensing device and generates a restoring force related to that error. 
Since 



(*i -*„)l 



k + cD\ 
mD 2 ) 






o’ 









Introduction 1 75 



/ cD + k \ 

X °~\mD 2 +cD + k) Xi ’ 

or 

x °_( cD + k \ 

Xj \mD 2 4- cD + kj 

This is the TF of the dynamic system with feedback, that is, it is the closed loop TF. 

The system response, or output for a sustained harmonic input motion 
X; = X { cos vt, has already been discussed, see section 2.3.2. If the output x 0 is assumed 
to be X D cos (vt — </>), then substitution into the closed loop transfer function (x 0 /Xj) 
gives 

— mv 2 X a cos (vt — </>) — cvX 0 sin (vr — (f>) + kX D cos (vt — 4>) 

= — cvX j sin vt + kX { cos vt. 

This equation can be solved by using the phasor technique (Fig. 5.7). 




Fig. 5.7. Phasor diagram. 



v rad/s 



From the phasor diagram, 

(kXf + (cvXf = (kX Q - mv 2 X 0 ) 2 + (cvXf. 

Hence 



X iy J(k 2 + (cv) 2 ) = X 0 J((k - mv 2 ) 2 + (cv) 2 ). 

This equation gives the amplitude of the sustained oscillation. However, there is 
also an initial transient vibration which is given by the complementary function; 
although this is a damped oscillation which dies away with time, it is often important 
in control systems, particularly if accurate positioning is required in a short time. 
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In the study of control system dynamics, the response of the system to a range of 
types of input must be considered; for example the impulse, step, and lamp sWwn 
in Fig. 5.8. Impulse excitation has been discussed in section 2.3.6. 







: Impulse 

dr 






▲ 


» r 




Step 

w t 


w *■ 

L 


Ramp 

c ► t 



Fig. 5.8. System inputs. 



For the system considered, if the input x ; is given a step change so that x- = X, 



( k + cD V 

\mD 2 + cD + k ) Xi 



= ( 1 + ^ D 2 

k k 



1 + 7 D )X 
k 




= X. 



That is, in the steady state x D = X, so that there is no error between the input and 
the output. However, the complementary function is given by (mD 2 + cD + k)x 0 = 0, 
so x„ = Ae~ l ’‘ ot sin (ft\/(l — C 2 )t + <j>). 

That is, x„ may also contain an initial damped oscillation, so that the system 
response to a step input would be as shown in Fig. 5.9. 
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Displacement 




The amount by which x D exceeds x t on the first oscillation is called the first overshoot. 
In a stable system, that is one in which the transient dies away, this will be the 
maximum overshoot. 

If the input is a ramp, x s = fit, 

so that 

<°=H d+ •■■X ,+ i d ) (w 

=(H d+ ')(^ + ^) 

= pt + -p --P = pt. 

That is, there is no steady-state error. However, the initial transient exists, so that 
the system output response is as shown in Fig. 5.10. 



Displacement 




Fig. 5.10. Ramp input response. 
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There are many needs and applications for systems that will move a remote device 
to some prescribed position in an acceptable and controlled way. The input for this 
movement may be a manual action such as turning a knob or moving a lever, or it 
may come from the motion or signal from some other device. Some systems have 
relatively low power at the input in comparison with that dissipated at the output; 
essentially, therefore, they are error-actuated power amplifiers, and they are usually 
referred to as servomechanisms, or just servos. 

5.1 THE SIMPLE HYDRAULIC SERVO 

Hydraulic remote position control (RPC) systems are used extensively; they generally 
rely on some form of supply control valve feeding actuators, pumps, and hydraulic 
motors. The type of control valve used can vary, but the most common is the spool 
valve which is shown diagrammatically in Fig. 5.11. 



Input 

position 




The input position of the spool may be controlled by direct mechanical linkage, 
or by electrical means. When the spool is centrally situated, as shown, both ports are 
closed and channels A and B are cut off. When the spool is moved to the right, B is 
connected to the high-pressure supply and when the spool is moved to the left, A is 
connected to the high-pressure supply. For small valve openings, it is usually assumed 
that the rate of flow of fluid is proportional to the valve opening, and that the fluid 
is incompressible. Some improvement in performance can be obtained by using the 
three-seal spool valve shown in Fig. 5.12. 



5.1.1 Open loop hydraulic servo 

Fig. 5.13 shows an open loop hydraulic servo. The input controls the position of the 
spool valve which directs fluid under pressure through a port to one end of the 
working cylinder. The working piston is connected to the controlled element so that 
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Input 

position 







Fig. 5.12. Three-seal spool valve. 



its position may be considered to be the output. 




If the area of the piston is A and the port coefficient is b, the flow equation gives 
Q = i4Dx 0 = bx i 



where Q is the rate of fluid flow. 
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Thus - “ ({Xs)*" 

that is, the system integrates the input. 
If a step input is applied, 

at t < 0, Xj = x 0 = 0, 

and at f ^ 0, x i = X. 

Thus when t ^ 0, 




so that the output increases with the time - as shown in Fig. 5.14. 



Displacement 




This type of response is usually undesirable in a control system, so the system is 
improved by applying feedback to close the control loop. 

5.1.2 Closed loop hydraulic servo 

Because of the limitations of the open loop servo performance, feedback of the 
output is usually added, as shown in Fig. 5.15. As the output reaches the desired 
value it acts to move the spool valve and thereby close the port. 

For the spool valve displacement y, from the figure, 

Xj + x 0 _ Xj-y 

l + m l ’ 

so that 

y = x { - (1/(1 + m))(x 0 + X;). 

The flow equation gives by = ADx„. Thus if y is eliminated, 

ADx a = — ( bl/(l + m))x 0 + (bm/(l + m))x ; . 



Control rod 
geometry 
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Output 



that is, 

[1 + ((/ + m)/bl)AW]x 0 = (m/f)x i? 
or 

[1 + TD]x 0 = (m/ l)x { . (5.2) 

where T = ((/ 4- m)/bl)A and has dimensions of time. 

T is known as the time constant of the system. Equation (5.2) is the equation of 
motion, and it relates the output to the input. It is sometimes convenient to consider 
the error e between input and output where £ = x ; — x G . 

Thus from (5.2), (1 + TD)x 0 = (m//)(£ + *<>)’ or 

X D = (l- ( m/I) + TU)-\m/fyL 

Hence a block diagram can be drawn to represent the system, as shown in Fig. 5.16. 



Inpur € 

^ 


m 

y 




(1 -SU 7D)- 1 


Output 


X. ^ 

1 


y — ^ 






x o 











Fig. 5.16. Servo block diagram. 
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It should be noted that in the idealized hydraulic servo considered, many factors 
have been neglected such as friction, slackness in connections, non-linearity in the 
spool valve, pressure losses, compressibility of the fluid, and flexibility of the metal 
parts. In a real system some of these factors may have to be taken into account if 
they have a significant effect on the system response. 

Response to a step input 

A step input is when x ; is suddenly increased to a value, say X, and when t < 0 
Xj = 0, and when t ^ 0, x ; = X. 

From equation (5.2), for the output x D the P.I. = ( m/I)X and C.F. = Be~' /T . Thus 
x Q = (m/l)X + 

Since x 0 = 0 when t = 0, B = — (m/l)X, hence 
x 0 = (m/l)X[\ - e“' /r ]. 

That is, x„ increases exponentially with time to a value (m/l)X; the system time 
constant Thas considerable effect on this response (see Fig. 5.17). 




After time T, the time constant of the servo, 

* 0 = (j)*d - e" *) = (jjx( 0.632), 

that is, the output reaches 63.2% of its final value at time T. It can also be shown 
that the output reaches 86.5% of its final after a time 2T, and 95% after a time 3 T. 
The factor ( m/l) amplifies or attenuates X to give the final x 0 value reached. This may 
be compared with the servo without feedback when x 0 = (b/A)Xt, and the output 
steadily increases with time. 
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Response to a ramp input 

A ramp input is an input which increases uniformly with time. Thus x { = fit is a ramp 
input, where /J is a constant. 

From equation (5.2), for the output x 0 , 

P.I. = (m/[)pt - (m/Q7fl and C.F. = Be“ t/T . 

Thus 

x D = Be~ tlT + ( m/l)pt — (m/l)Tp. 

If x 0 = 0 when t — 0, B = ( . Hence 
x„ = (m/l)Tp(e~ t/T - 1) + ( m/[)pt . 

For the case when m = I, 

x D = 7p(e~ ,/T - 1) + fit. 

Thus x 0 increases exponentially with time towards a value /?r ( = Xj). The error 
between the output and the input is 7 ]?(e _I/T — 1). 

Even when t = oo, this error is — Tp. This is known as a steady-state error because 
it persists after steady conditions have been attained (see Fig. 5.18). 




Response to a sinusoidal input 

For a sinusoidal input Xj = X { cos vt. 

For x Q the C.F. is an exponential decay, and this can usually be neglected. A P.I. can 
be assumed, x 0 = X 0 cos (vf — </>). 

Substituting this solution into equation (5.2) gives 



X Q cos (vr — 0) — vTX 0 sin (vt — cf>) = (m/[)X i cos vt. 



This equation can be solved by drawing a phasor diagram as shown in Fig. 5.19. 
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Fig. 5.19. Sinusoidal input response. 



Thus 

(m/[) 2 Xf = (14- v 2 T 2 )X 2 
or 

XJXt = (m/l)/y/(l + v 2 T 2 ) and (t) = tan~ 1 vT. (5.3) 

It is often desirable to show graphically how the steady-state response to a sinusoidal 
input varies with frequency. There are several methods of doing this; one is to draw 
a Bode diagram. The Bode diagram is a plot of gain (XJX J against frequency. The 
magnitude of (XJX ( ) is plotted on a log scale, usually dB. 




Fig. 5.20. Bode diagram. 




Sec 5.2] 



Modifications to the simple hydraulic servo 185 



The gain (or attenuation) in dB = 20 log 10 {XJX), so for equation (5.3) 

Gain (dB) = 20 log l0 (m/Z) - 101og lo [l + v 2 T 2 ]. 

Fig. 5.20 shows the diagram for the case m — l. 

Further discussion of the Bode analysis technique is given in section 5.7.2. 

An alternative diagram is the harmonic response locus which can be drawn by 
plotting gain and phase as a polar diagram. Since 

X, = (l/m)Xj( 1 + v 2 T 2 ) = [{IXJm) 2 + (lX 0 vT/m) 2 V /2 ’ 

the locus is as shown in Fig. 5.21. 




v = 0 




Fig. 5.21. Harmonic response locus. 



5.2 MODIFICATIONS TO THE SIMPLE HYDRAULIC SERVO 
5.2.1. Derivative control 

An intelligent human controller would take account not only of the instantaneous 
value of the input but also its rate of change, or the rate of change of error. This 
would help to avoid overshooting and would improve the system response. The simple 
servo can be modified to act in a similar way, as shown in Fig. 5.22. 

From the control rod geometry, 
the valve displacement 

y = (*, - v)/2. 

The flow equation gives 

by = ADx 0 

and force balance gives 

k(rx a — r) — cDv — 0. 

It is required to eliminate v and y and to give x„ as_/fXj). From (5.6), 
v = ( kr/(k + cD))x 0 . 

Substitute in (5.4): 

y = (x,/2) - (krj(k + cD)Xx 0 /2) 
and in (5.5) 
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High pressure 




ADx a = bxj 2 - (b/2)(kr/(k + cD))x 0 , 



that is 



(2AcD 2 + 2kAD + bkr)x Q — (bk 4- bcD)x ; 



or 



((2y4c/h/cr)D 2 + (2A/br)D + l)x Q = (1 + c!kD)x-Jr. 



(5.7) 



Comparing equation (5.7) with equation (5.2) which was derived for a similar system 
without a velocity element, it can be seen that both sides of the equation of motion 
have gone up by one order. Thus when x ; is changed the valve setting y is not reduced 
by the full value corresponding to the change in x 0 , because the faster x ( and x 0 
change, the less y is reduced. 

From equation (5.7), 



, , 2A 2Ac -1 

- |1+ 67 d+ m7 d 






so for a step input, Xj = X, 




The C.F. of x D is a damped oscillation or decay, so the response to the step input is 
as shown in Fig. 5.23. 
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For a ramp input, x i = fit, hence, 



x 



o 






pt cp 2AP 
r + kr br r 



so that a steady-state error exists. The initial oscillation is given by the C.F. as above. 

The steady-state response to a sinusoidal input can be found by drawing a phasor 
diagram. If x ; = X, cos vt, x„ = X 0 cos (vt — <f>), then from Fig. 5.24. 

x = V(1 + (cv/k) 2 ) 

>/((r — ( 2Ac/bk)v 2 ) 2 + (2 Av/b) 2 ) r 




Fig. 5.24. Phasor diagram. 
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At low frequencies v ~ 0, so X 0 -» X-Jr and <f) -» 0. 

At high frequencies v ~ oo, so X a -» 0 and 4> -* n/2. 

A detailed comparison of the simple system performance with that for a system 
with derivative control added depends on the system parameters; but, in general 
derivate control allows a faster transient response, with overshoot if too fast. 

5.2.2 Integra] control 

Both the simple system and that with derivative control will give steady-state errors 
when the input is a ramp, x s = fit. Real systems, in which fewer assumptions are made 
than here, can give steady-state errors with other kinds of input. To eliminate these 
errors a correction can be added to the system which is proportional to the time 
integral of the error: such a system is shown in Fig. 5.25. It can be seen that compared 
to the system with derivation action. Fig. 5.22, k and c have been interchanged. 

High pressure 




Equations (5.4) and (5.5) apply, but force balance gives cD(rx 0 — v) — kv = 0 and 
thus 



/ 2Ac \ f be bk 1 \ 

{ + 2 Ak + ber J X ° = \2Ak + her + 2 Ak + ber Dp ( ' ) 

Comparing equation (5.8) with equation (5.2) shows an extra term proportional to 
the integral of x { . The response to a step or sinusoidal input can be found from 
equation (5.8). 
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Example 43 

A servomechanism incorporating an hydraulic relay with displacement feedback 
through a dashpot and spring assembly is shown below. 




The velocity of the output ram ( Dx 0 ) is equal to K times the movement of the pilot 
valve from its neutral position. Given that the inertia of all parts and the effect of 
friction (other than at the dashpot c) may be neglected, determine the equation of 
motion of the system. 

For the control rod. 




where v is the spring displacement and y is the movement of the spool valve. 
Thus 



y — v X; — v 
L = L+m 



and 







L+ m 



L+ m 
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Force balance for the spring and dashpot gives 



kv — cD 




so that ( k + cD)u = 




The flow equation gives Dx 0 = Ky. 
Substituting for y and v gives 




which can be rearranged to give the equation of motion as 




Kmr \~| _ / KL V fe + cD \ 
(L+ m)(r + s))_ X ° \L+wJ\ c 



Example 44 

The figure opposite shows an hydraulic servo, in which it can be assumed that 

(a) all friction (except that intended to be present in the dashpot which has a viscous 
damping coefficient c) is negligible, 

(b) the oil is incompressible, 

(c) the inertia forces are negligible, and 

(d) the volume flow per unit time through either port in the valve is b times the 
displacement of the valve spool from its neutral position. 

Derive the transfer function relating the output displacement z to the input 
displacement x, this being measured from the position corresponding to zero oil flow. 
Find the value of z as a function of time following a sudden change in x of magnitude X. 

The flow equation gives by = ADz, where y is the displacement of the spool valve. 
If the spring extension is p, 

2(y + p) = (x + p) 



X 




kp = cD(nz — p) P 
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P I. is 



bk 



X = 



bk 



(cnb + 2Ak)D(l + \cnb + 2Ak 

cnb + 2Ak 



I Xt, 



C.F. is P + Q exp 



2 Ac 



so 



z = P + Qex p 



cnb + 2Ak 
2Ac 



t + 



bk 



cnb + 2Ak 



l Xt. 



Now when t = 0, z = 0, so P + Q = 0. And when t = 0, p = 0, so Dz = (6/2>l)A\ 
Thus 



fcX _ . _/ cnb + 2 Ak\ ( bk 

2A~ Dz = 0 + F \~2Ar~) + (^bTTAk 



SO 



1 , fc 

2^c 






cnb + 2 Ak 2 A 



or 



fe 2 c 2 nX 
(cnh + 2/1/c) 2 

Thus, 

/ bkX \ b 2 c 2 nX ( 

Z ~ cnb + lAkf + (cnb + 2Ak) 2 j 1 _ CXp 



/ cnb + 2 Ak 
\ 2 Ac 




Example 45 

In the hydraulic position control servo shown, the desired input position is converted 
to a voltage and compared with the output position voltage V a . Any error is applied 
to a power amplifier and position transducer to move the spool valve controlling the 
oil flow. The oil is delivered to the cylinder so that the piston acts on the load as 
required; a position transducer generates the signal V a as shown. 

With the parameters given in the figure, find the relationship between the output 
position and the input position, 

Xi 

If an input of 5 cm/s is applied, find the error between the input and output after 
2 s; neglect transient effects. Also calculate the required coefficient of the output 
displacement transducer for this error to be a minimum, and the size of this minimum. 
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High pressure 
oil 



Desired 

position 



Input position 
transducer 
1 v/cm 



Position 

Power transducer 
amplifier i mm/0.1 A 
400 mA/V 




Output position yA v Y/ 

transducer 

5 V/cm Piston 

area 20 cm^ 



Spool valve opening = — VJ 400 



100 



mm. 



Oil flow 



= 20 D* 0 cm 3 
= 5(Vj - VJ4 cm 3 . 

K = 5X 0 and V { = X it 

20 DX a = 5(*i - 5X 0 )4 



or 



X„ . 0-2 

X i 1 + 0.2 D 



Input is 5 cm/s, so X t = 5 1 cm. 
In the steady state 
0.2 



X„ = 



1 + 0.2D 



X, 



= 0.2[1 -0.2D + - ]5r 



- t - 0.2. 
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Hence steady-state error = — X 0 

= St-(t- 0.2) 

= 4r + 0.2. 

After 2s, ss error = + 8.2cm. 

For minimum error, let output position coefficient be a V/cm. 
Then 



V = 



and 



20 DX a = 5(X. t - a* 0 )4. 



Hence 



DX 0 = X, - a* 0 . 



Thus 



X n 1 



X., D + a’ 

1 ( D 

and steady-state output X a = - 1 (- • • 

a \ a 

5 5 

= - t ~ -j. 
a a z 

When t = 2, the steady-state error is 

£ = 5x2 — - x2h — 
a a 

for minimum e, de/da = 0. 

That is, 

J J ^ 10 10 

de/da = 0 = — , r, 

a a 



5t 



so a = 1. 

Therefore minimum steady-state error =10—10 + 5= + 5 cm. 



5.3 THE ELECTRIC POSITION SERVOMECHANISM 

Electric position servos are widely used in control systems; an example of turntable 
position control has already been considered in Chapter 1. 

Open loop systems are available, but in practice the operator often has to act as 
the feedback to ensure correct positioning. Usually, therefore, only the closed loop 
servo is considered. 
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5.3.1 The basic closed loop servo 

The elements of the basic closed loop position servo are shown in block diagram 
form in Fig. 5.26. The position of the load is compared to the desired position, thereby 
creating an error signal. The amplifier and motor respond to the error signal and 
act on the load position. 



Energy 

Amplifier 
and motor 



Inertia load 

and viscous damping 




Fig. 5.26. Closed loop position servo. 



The error signal 0 is the difference between the desired position (0;) and the actual 
position (0 O ) that is, 0 = 0; — 0 O . The amplifier and motor are such that the torque T 
produced by the motor is K9. For the load, 

T = (JD 2 + cD)0 o , 

since it is considered to comprise an inertia load with viscous damping. Thus 
m - 0 O ) = (JD 2 + cD)0 o 



or 



(JD 2 + cD + K)0 O = K0;. 

This is the same equation of motion as found for forced damped vibration, and 
therefore the solution is the same for a given input. 

Response of servo to a step input 0; = X. 

For 0 O , C.F. = v4e“ sin (o>V(l - ( 2 )t + </>), 
where £ = c/c c and co = ^J(K/J) rad/s, and 

P.I. = [1 + (c/K)D + (J/K) D 2 ] -1 * 

= [1 — (c/K )D + • • • -\X 

= X. 



Thus 



0 O = Ae Cto ‘ sin (cu^fl — £ 2 )t + <)>) + X. 



(5-9) 
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The response is shown in Fig. 5.27. 




Clearly the fastest response is when the damping is low. Unfortunately this is 
accompanied by a large first overshoot and a poor transient response, that is, the 
oscillations take a long time to die away. Increasing the damping slows the response 
but reduces the overshoot and improves the transient. Critical damping gives the 
quickest response without overshoot. 

To find the overshoot, dQJ&t is found and made zero and hence the times to 
overshoot are obtained. Substituting these times into (5.9) gives 9 a . See Example 46. 
There is no steady state error. 

Response of servo to a ramp input 9 i = fit 
For 9 a , C.F. as above (same system), but 

P.I. = (1 + (c/K) D + (J/K)D 2 r l fit = fit- ( c/K)fi 



= 9 i -(c/K)fi. 



Hence there is a steady-state error ( c/K)fi . 

Thus high damping to achieve a good transient response also results in a large 
steady-state error to a ramp input. 



Response of servo to an harmonic input 9 { — A sin vt 

The analysis is exactly as for the forced damped vibration system, so 



» = 

° V((K - Jv 2 ) 2 4- c 2 v 2 ) 



cv 



sin (vt — <fi), where (fi = tan 1 



K — Jv 
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Hence there are gain and phase errors which are functions of the system parameters 
and the frequency of excitation. 

Because of the shortcomings of the basic system, it is often modified to provide an 
improved response. 

Example 46 

In the closed loop system shown in Fig. 5.26 the input is given a step increase 9 V 
Find the output, 9 a , as a function of time, assuming that both 9 a and d0 o /dt are zero 
at t = 0. Hence find the time to the first overshoot and the corresponding magnitude 
of 9 a . 

From equation (5.9), 

0 o = 9 { + e~ i<0 'lA sin a\/(l — £ 2 )r + B cos a\/(l — £ 2 )t]. 

Now, when 9 0 = 0, t = 0, so 6 { = — B, and when d9Jdt = 0, t = 0, so 
0 = - + coV(l - C 2 )A 



A va - c 2 ) e " 



Hence 



6=9 - 9 e-^‘ 



c_ 

.Vd- 



sin co^I — £ 2 )t + 



COS £0^/(1 — £ 2 )r 



that is, 



e ° = 9 [ l _ v / (1 _ c 2 ) e ~ c " [sin K/a - e 2 )* + 0)] 



where 



tan = 



ya - c 2 ) 



Now at overshoot, time is given by d9Jdt = 0. 
That is, 



0 = + £ (oe ^ ^ sin (co^/fl - £ 2 )t + 0) 



^(1 W * 1 ~ C2) ^ cos (" V ( 1 - C 2 )* + 0 )]- 



Hence 



Va - c 2 ) 



= tan (co^C 1 - C 2 ) f + 0) 



= tan I 0)^(1 — £ 2 )r + tan 
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or — £ 2 )t = 0, n, 2n. . . 
so that the time to the first overshoot is 
n 

o»V(i - c 2 )' 

This is one half of the period of the damped oscillation. 
At this time. 



Oi 



V(i -{ 2 )“ p Wa-c J ) 



— £coti 



sin (n + (f)) 



= 0 - 



0 t 



exp 



- (,71 



l( - sin 4>) 



' V(i-c 2 ) K Wd-c 2 ) 



= $ + exp ( - ^ L/(l - £ 2 ) 



7i T Vd - C 2 ) F W(1 - c 2 >. 



Thus 



9=9; 



1 + expi 



Cn 



V(1 - C 2 )7 J 

Now if in a given servo, co = 5 rad/s, ( = 0.1 and 9 { = 10 deg. 
Time to first overshoot is 



5V(1 -0.1 2 ) 



= 0.63 s 



and 



= 10 



1 + expl 



0.1 7t 



V(1 -o.i 2 ) 



= 17.3° 



Increasing £ reduces 9 0 but increases the time to the first overshoot. 



Example 47 

The angular position 9 0 of a turntable is controlled by a closed loop servomechanism 
which is critically damped by viscous friction. The moment of inertia of the turntable 
is 400 kg m 2 , and the motor torque is 3.6 kNm/rad of misalignment between 0 O and 
the desired position 9 t . Given that ^ = 10 rev/min, find the steady-state position error. 

The equation of motion is 

(JD 2 + cD + K)9 a = K9 j 



so that 
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"• = [ 1 + ^ D + I D2 _ 

If 0; = /It, d a = Pt-{c/K)P, 



9 .. 



so the steady-state error is 



c 

K 



P- 



Now (cjK)P = 2 since damping is critical and c c = 2 y/{JK). 

Hence 



6 



ss 




= 40°. 



Response of servo to sudden loads 

Servomechanisms are sometimes subjected to sudden loads due to shocks and 
impulses, such as those arising from gust loads and impacts. Fig. 5.28 shows a servo 
with a sudden load T L applied. 



Inertia load 




Fig. 5.28. Servo with external load, r L . 

In this system, 

m - e o ) +T l = ( JD 2 + cD)0 o , 

so that 

(JD 2 + cD + K)6 0 =T l + K9 ;. 




Fig. 5.29. Servo response to sudden load. 
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If the input 6 i is zero, the response to the sudden load T L gives a steady-state error TJK. 

The complementary function gives the transient response so that the output position 
is as shown in Fig. 5.29. 



Example 48 

In a simple servo mechanism, J = 0.1 kg m 2 , K ~ 20 N m/rad error, and c = \c c . 
Given that a step input is applied, 

(i) find the frequency of the transient oscillation of the output, and the time to rise 

to the first maximum; 

(ii) find the position lag if an output of 10 rev/min is applied; and 

(iii) calculate the steady state error when a load of 1 Nm is applied. 



(i) 





= 14.1 rad/s, 



and 



C — 0.5 



so 

co v — - c 2 ) - 14.1 V[1 - (0.5) 2 ] = 12.3 rad/s, 

and frequency of oscillation,/ = 12.3/277 = 1.95 Hz. 

The time to the first maximum = (i)(l If) = 0.256 s. 



(ii) Steady-state error = — / 

/C 



y/(KJ) n = P 

K 1 a) 



2tt:10 



60 x 14.1 



= 0.074 rad = 4.2° 



F 1 

(iii) Steady-state error = — = — = 0.05 rad = 2.9° 

/C zO 



Example 49 

The block diagram of a satellite tracking system is shown 
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External 




The transfer functions of the amplifier plus motor, tracking device and sensor are 
as shown. The system is subjected to the desired position input, and also an external 
torque T L which is due to wind loads acting on the structure. 

(i) When T L = 0, find the closed loop transfer function, 0 O /0 S , 

(ii) Find K, and K 2 if 

(a) there is to be no steady-state error with a step position input, and 

(b) when a ramp position input of 2 rad/s is applied the steady-state error is to 
be 0.056 rad. 

(iii) A sudden steady wind can be approximated to a suddenly applied torque 
T l = 2 N m. Find the magnitude of the sustained displacement. 



(i) When T L = 0, the equation of motion is 

K 1 V 100 



(0i - 



1 + 0.2D/ \D + 8D + 100 



= 0 o . 



so 






and 



0 = = e, 



(D 2 + 8D + 100){1 + 0.2D) + 100 K X K : 

Took, 



Hence CLTF (closed loop transfer function), 

0c/ 100K^ 

6 i V(D 2 + 8D -I- 100X1 + 0.2D) + 100K 1 K 2 

( 100K, 



100 + 100K,K 2 + 28D + 2.6D 2 + 0.2D 3 
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(ii) With a step input 6 t = X, 

9 = _ — 122^1 (l- -D-)X 

° 100 + 100X^2 

so that the steady-state error 9 ^ is given by 
100K 1 _ 



X \ 1(X) + 100Ki x 2 



For Bgs to be zero, K x = 1 + K X K 2 . 
With a ramp input 9 l = pt, 

K j (. 28D 



0 = 



1 + k x k 2 



1 - 



Pt-(; 



100 + mK x K 2 
28 



D 2 



fit 



1 + K x K 2 \ r ' '100 + lOOK^)' 
Now K x = 1 + K x K 2 , so 



fib 



e 0 = pt- 



28 



P- 



100X/ 

When P = 2 rad/s, require 9$$ to be 0.056 rad. 

That is 

28 - * — ■ = 0.056, 

100 x 

so that 

K x = 10 and hence K 2 = 0.9 

(iii) With a sudden load T L applied the equation of motion is 



K2 °°{ 1 + 0.2D 



+ T l 



100 



D 2 + 8D + 100 



= B 0 - 



That is 



„ /D 2 + 8D + 100 K x K 2 
0A — + 



1 + 0.2D 






100 

Now K x K 2 - 9, so 

^(D 2 + 8D + 100X1 + 0.2D) + 900 



9, 



100(1 + 0.2D) 



or 
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0 „ - 



100(1 + 0.2D) 



1000 + 28D + 2.6D 2 + 0.2D 3 



100 28 

*1 + 0.2DX1- — D -)T l . 



1000 ' 1000 
When T l = 2, the sustained displacement is (100/1000) x 2 = 0.2 rad. 



5.3.2 Servo with negative output velocity feedback 

If viscous damping is not great enough, the damping effect can be simulated by 
providing a signal proportional to D0 O . This is shown in Fig. 5.30. 




Fig. 5.30. Servo with output velocity feedback. 



For this system 



m - 0 O - C,D0 o ) = (JD 2 + cD)0 o , 



that is 



(JD 2 + (c + KCijD + K)0 o = K0 

Thus the value of C, can control the total damping, and the step and transient 
responses can be improved at the expense of the steady-state error with a ramp input. 

However, if CjD0 o is added to 9 instead of subtracted, the opposite is true, but 
there is a danger of negative damping occurring with the associated unstable response 



Example 50 

A certain servomechanism is required to control the angular position 0 o of a rotatable 
load of moment of inertia J. The rotation of the load is subject to a viscous friction 
torque c per unit angular velocity. The mechanism has velocity feed back such that 
the motor torque is 





where £ is the error between the desired angular position input 0 i and the output 
position 9 a , and K and C 1 are constants. 
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(i) Draw a block diagram of the system and establish the equation of motion. 

(ii) Find the steady-state positional error when the input signal has constant velocity 
given by 




(iii) When a step input is applied the response overshoot must not exceed 10%. Find 
the minimum value of the damping ratio required, and express this ratio in terms 
of the given parameters. 

(i) The block diagram is as below: 




Thus 

10, -0 o - C.DdJK = [JD 2 + cD.]0 o , 
and the equation of motion is 
JD 2 9 0 + (c + KC l )D0 o + K0 o = K6 { . 



(ii) Input 0 { = (it. 



0 = 



1 + 



c + KC A d + J » 2 



fit 



K 

c + KC 1 

~k 



K 



P- 



steady-state positional error = 



I St 



C + KC r 

K 



\P- 



(iii) Damping ratio C = 



‘c’ 

C c 



c + KC 1 
2 J(KJ) ‘ 



At first overshoot, 
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Thus 



0.1 = exp[- C^/VU ~ C 2 )]> 

so that 



2.3 = 



Cn 

Vu - C 2 ) 



Hence 



c = 0.59 



C + KCy 

2 yJ{KJ) ' 



Example 51 

A linear remote position control system with negative output feedback consists of a 
potentiometer giving 8 V/rad error to an amplifier, and a motor which applies a 
torque of 3 N m/Y to the load. The load has an inertia of 6 kg m 2 and viscous friction 
of 12 N m s/rad. 

(i) Draw a block diagram for the system and derive its equation of motion. 

(ii) Calculate the maximum overshoot in the output response to a step input of 
2 rad, and 

(iii) Given that a tachogenerator is employed to provide negative output velocity 
feedback, derive the new equation of motion and calculate the velocity feedback 
coefficient needed to give critical damping. 

(i) The block diagram for the system is 




Now (0j - 6 0 )KG = (JD 2 + cD)0 o , 
So the equation of motion is 
(JD 2 + cD + KG)9 0 = KG6 { . 
where 

K = 8 V/rad error, 
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G = 3 Nm/V, 

J = 6 kg m 2 , and 
c = 12 N m s/rad, 

(ii) With a step input, overshoot is dp _s,I/v(1 
where 

c 12 1 

^ 2 J(KGJ) 2^(8 x 3 x 6) 2‘ 

Thus 

e — 1/2 ti/ v '[1 - 0/2) 2 ] _ g- 1813 _ QJ63 

and overshoot = 2 x 0.163 rad = 0.326 rad, for a 2 rad input. 

(iii) With output velocity feedback, the block diagram is 




So that 

[0 ; -9 0 - C.D6JKG = [JD 2 + cD]0 o 
and the equation of motion is 
[JD 2 + (c + c,KG)D + KG]0 o = KG6 
For critical damping 
c + c,KG = IjiKGJ), 
that is 

12 + (Cj x 8 x 3) = 27(8 x 3 x 6). 
Hence 

c, = 0.5 (N m s/rad)/(Nm/rad), 
or 



Cj = 0.5 s. 
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5.3.3 Servo with derivative of error control 

An alternative modification to the simple feedback servo is to add derivative of error 
control in which the signal is fed forward. 

This system is shown in Fig. 5.31. 




Fig. 5.31. Servo with derivative of error control. 



For this system 

m -0 O + c 2 m, - 0 O )) = (jv 2 + cD)d 0 , 

that is 

(JD 2 + (c + KC 2 ) D + K)0 o = K( 1 + C 2 D)0i. 

Because the damping term is increased the transient response dies away more quickly. 
For the response to a ramp input 0, = /3r, 

pI K(l+C 2 D)Pt 

K(l 4- (cjK + C 2 )D + ■ ■ ■ ) 

-(1 ~(c/K + C 2 )D + ---W + C 2 P), 

= pt- {C/K)P- 

that is, (c/K)p is the steady-state error due to the damping; this has not been affected 
by derivative control. Thus derivative control is a method for temporarily increasing 
the total damping term and hence the rate of convergence of the transient response, 
without increasing the physical damping or the steady-state error. 

5.3.4 Servo with integral of error control 

A further alternative is to feed forward an integral of the error signal as shown in 
Fig. 5.32. 

For this system 

Kid, ~0 O + (KJ D{6- dj) = (JD 2 + cD)d a 

that is, 

(JD 2 + cD + K + KKJ D)0 O = K(1 + KJT>)d, 
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Fig. 5.32. Servo with integral of error control. 



or 

(JD 3 + cD 2 + KD + KK l )0 o = K( D + KJ0,. 

For the response to a ramp input = [it, 

PI = K( k^ + + d K + - ) # = i = 

that is, there is no steady-state error. Although this is a very desirable quality, it is 
possible for this system to become unstable. This can be demonstrated by considering 
the response to an harmonic input. 

The left-hand-side of the equation of motion is third order, but the response to a 
sinusoidal input can be found by the phasor technique. 

If 0; = A cos vr and 0 o — B cos (vf — <£) in the steady state, then the equation of 
motion can be represented by Fig. 5.33. 




Fig. 5.33. Phasor diagram. 



AjdKv ) 2 + ( KK J 2 ) 
V((Kv- Jv 3 ) 2 + (KK 1 -cv 2 ) 2 ) 



cos(vf — <f>). 



If the denominator is zero, B is infinite and the system is unstable. For the 
denominator to be zero, 
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Kv - Jv 3 and KK X — cv 2 . 

Thus 

c = JK X and v 2 = K/J. 

This is a critical relationship for a stable response. 

The performance of the main forms of the electric position servo are compared in 
Table 5.1 overleaf. 

Example 52 

In the position control servo shown in block diagram form below, the amplifier 
voltage gain is 100 and the error detector produces 2 V per degree of error. 




The motor produces a torque of 0.1 N m per 100 V applied. The inertia of the 
load is 0.1 kg m 2 , and the viscous damping is 10 N m s/rad. Find the damping ratio 
of the system. 

Given that the input speed is 5 rad/s find the steady state error. 

Motor torque = 0.1/100 x 2 x 57.3 x 100 Nm/rad error 
= 11.46 Nm/rad error. 

Torque applied to load = 250 x 11.46 = 2865 Nm/rad error. 

For damping ratio, 

£ = c/c c = 10/2^/(2865 x 0.1) = 0.295. 

For = 5 rad/s, 

steady-state error = 10/2865 x 5 x 57.3° = 1°. 

Example 53 

An engine drives a load of constant inertia which is subjected to external moments 
which vary unpredictably. To limit the speed changes resulting from load changes, 
a governor is fitted. This measures the engine speed and changes the moment exerted 
by the engine to try to keep the speed constant. The changes are related by the equation 

(1 -I- TD)m = — ks, 

where T is a time constant, 

m is the change in engine moment from an initial steady value, 
k is a constant, and 

s is the change in speed from an initial steady value. 






Table 5.1 Comparison of main forms of electric position servo 



Input 




Basic system 


Basic system with output 
velocity feedback 


Basic system with 
derivative of error 


Basic system with 
integral of error 


Step 


Response 




Perhaps unstable 


Exponential decay oscillation, rate of decay increasing 


s.s. error 


MIT 




HIL 


Ramp 


Response 


\ 


Perhaps unstable 


7 

Exponential decay oscillation, rate of decay increasing 


s.s. error 


Proportional to 


Larger than D0j 


Proportional to Dfl; 


NIL 


Harmonic 


Response 


Gain and phase errors proportional 
to system parameters and frequency of 
excitation. 


Perhaps unstable 
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Given that the load moment increases by a step function of magnitude m L and then 
remains constant, derive an expression for the consequent speed change as a function 
of time. If the inertia of the engine and load is 4.0 kg m 2 , T is 0.35 s, and k is such 
that the damping of the closed loop is 0.6 of critical, find the final drop in speed 
following an increase of load, m L , of 1 1 N m. 



Inertia, J 




Increase 

m. 



Since s is the change in speed, 
m — m L = JDs. 

Now 

(1 + TD)m = — ks. 
so 

— ks 

m = — — = JDs 4- m, . 

1 + TD L 

For a step change in load m L , D m L = 0, hence 
JTD 2 s + JDs + ks = — m L . 

Solution comprises P.I. — mjk and C.F., so 



m, , 

s = — + e s "‘( A sin a>t + B cos cot). 

k 

Substituting initial conditions s = 0 and s = 0 at f = 0, gives 
£m L 



A - 



and 



R- mL 



thus 



s 



— - [— 1 + e C(0t (C sin cot + cos cot)]. 
k 



The steady-state error is — 

k 

From the equation of motion, c c = yJ(kJT), so 
c = 1.2 V(kJT) = J, 
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hence 



J 4 

1.447 ” 1.44 x 0.35 



7.936 kg m 2 /s, 



Thus steady-state error = 



11 

7.936 



— 1.386 rad/s. 



Example 54 

A linear servomechanism consists of a proportional controller which supplies a torque 
equal to K times the error between the input and output positions, and drives a 
rotational load of mass moment of inertia J and viscous damping coefficient (less 
than critical) c. 

Draw a block diagram for the mechanism and write down its equation of motion. 
During tests on a mechanism for which J = 5 kg m 2 it was found that 

(i) an external torque of 100 Nm applied to the output shaft gave a steady-state error 
of 0.2 rad, and 

(ii) a constant velocity input of 3 rad/s produced a steady state error of 0.3 rad. 
Given that a step displacement of 10° is applied to the mechanism when it is at 
rest, find the magnitude of the first overshoot. 




The equation of motion is 

JD 2 0 o + cD0 o + K0 O = Kdi. 

100 

(i) K = = 500 Nm/rad. 

(ii) 0.3 - x 3, so c = 50 N m s/rad. 

Since c c = 2j(kJ) = 2^(500 x 5) = 100 N m s/rad, 







